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CONTROL AND DYNAMICS 


Noise in Power Reactors—A Review 
of Experiment, Analysis, and Theory 


By Joseph A. Thie, Consultant* 


The dynamic behavior of reactors has been studied 
extensively for more than a decade utilizing the data 
of their randomly fluctuating variables. The term 
“reactor noise” is the accepted jargon for these time- 
varying random fluctuations in reactor variables, 
especially the ion-chamber current. Testimony to 
the widespread use of noise analysis has been given 
in numerous individual publications, from which 
collections and reviews have been made.!~ This re- 
view is based on a paper presented at the Japan— 
United States Seminar on Nuclear Reactor Noise 
Analysis held in Tokyo and Kyoto, September 1968. 


Characteristics of Noise Experiments 


Of the methods of obtaining dynamic information 
from reactors, stochastic techniques are among the 
more popular. Their attractiveness stems from the 
following features: 

¢ There is essentially no disturbance to the reac- 
tor, and this is quite favorable to operating policies. 

¢ Existing instrumentation and equipment are nor- 
mally used (although special data-recording equip- 
ment might be added). 

eComplex phenomena sometimes might require 
noise analysis as one of the modes of investigation. 

eAn experimental simplicity in noise tests has 
appeal. 

e Considerable precedent exists in monitoring for 
reactor stability by noise studies. 


Noise experiments have been performed in both 
zero-power reactors and power reactors. Table 1 
generalizes about these experiments. The multiplic- 
ity of available variables and the system complexity 
in power reactors are perhaps two underlying rea- 


*Address: P.O. Box 17, Barrington, Ill. 60010. 


sons for the research differences of Table 1. The 
future research potential is rated as substantial; a 
few specific suggestions are given below. 


Table 1 GENERAL COMPARISONS OF NOISE 
EXPERIMENTS ON POWER AND ZERO-POWER REACTORS 








Zero-power Power 
reactors reactors 

Motives Research Engineering 
Variables for analysis Few Many 
Environmental suitability Good Poor 
Relative number of experiments Extensive Average 
Data interpretation Easy Difficult 
Future research potential Average Substantial 





In contrast to zero-power reactor noise experi- 
ments in which only the neutron (or gamma) flux and 
control-rod fluctuations have been used, power- 
reactor experiments offer many more possibilities 
for variables to be studied. Table 2 lists those which 
have been used in noise experiments, along with sug- 
gestions for others. Investigations can have many 
facets when more than just neutron-flux noise is 
measured. Also, as one might expect, the amount of 
information obtained from a study is somewhat in 
proportion to the number of variables recorded 
besides the usual neutron flux. 


fable 2. VARIABLES FOR NOISE ANALYSIS 





Already used Other possibilities 





Neutron flux Mechanical displacement, velocity 
Temperature or acceleration 

Flow Local density 

Control-rod position Force and stress 

Pressure 

Valve position 

Pump speed 

Core gamma flux 
Coolant radioactivity 
Sound 
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Experimental Techniques 


To see the place of noise analysis among methods 
of investigating reactor dynamics, the latter are 
broadly classified in Table 3. Noise methods are the 
last two listed. It might also be mentioned that 
further subdivisions within this table are possible — 
such as categorizing according to data-gathering 
techniques or according to types of experimental 
equipment. 

The last method in Table 3 is found more frequently 
in the current literature than are the others because 


Table 3 
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more noise than the lower because of local steam- 
void fluctuations in its environment. A technique for 
expressing the extent of correlation at various fre- 
quencies is used here: the coherence is 0.94 at low 
frequencies but reduces to 0.7 at 2 cps. (Coherence® 
is a function of frequency defined such that it is unity 
for two signals having exactly the same fluctuations 
in a given frequency band but is zero if they are 
random with respect to each other.) 

An extensive listing of the power-reactor experi- 
ments to date is given in Table 4. Evidently the most 
popular (in number of reactors) has been the spectral 


EXPERIMENTAL METHODS OF OBTAINING DYNAMIC INFORMATION FROM POWER REACTORS 





Name 


External excitation 


Description 





Event analysis Any cause 
Rod (or other device) oscillation 
Pseudorandom excitation 


Noise analysis None 


Control rod (or other device) 


Control rod (or other device) 


Data analysis from a reactor transient, 
sometimes unexpected but usually planned 
Measurement of reactor response to 
oscillations 
Measurement of reactor response to 
almost random excitation 
Analysis of randomly fluctuating 
variables due to intrinsic causes 





of its simplicity and other reasons for popularity 
cited above. Often it consists of a Fourier analysis of 
the ion-chamber fluctuations, the result being the 
relative distribution of frequencies in the random 
signal. Recently, the trend among experiments is to 
analyze two or more signals simultaneously, obtain- 
ing mutual interactions as well as individual fre- 
quency structures. 

An example of the simultaneous analysis of two 
ion-chamber signals is given in Fig. 1. Individual 
noise spectra are given—one for the lower non- 
boiling part of the core, and one for the upper boiling 
zone of the same fuel channel. (The ordinate’s units 
are such that the integral of the curve over all fre- 
quencies is the square of the well-known root-mean- 
Square value of the power fluctuation.) At higher 
frequencies the upper ion chamber is seen to have 


Table 4 RESULTS ACHIEVED INA 




















1C 
2 Upper ion chamber 
= l }—__ 
= 
2 Lower ion chamber 
= Ol af 
we 
5 
ry 
@ = Coherence = 0.7 
3 _ Coherence = 0.94 
: | 
0.001 | 
0.1 1 10 
Frequency, cps 
Fig. 1 Results of a computer analysis of digitized ion- 


chamber recordings of the Pathfinder reactor at 190 Mw. 


T-POWER DYNAMICS STUDIES 





Method and excitation Reactors used 


Results 





Fluctuating control rod Kiwi-A3,’ MSRE,® Saxton! 
(pseudorandom) 
Noise analysis 


(no excitation) Calder Hall,'!? Chapelcross 


Dresden,'4 EBR-II,!5 EBWR,'§ 17 


BORAX-I, -II, -IV, and -V,9+10.1! 


Stability study, control-system study 


Stability study, control-system study, 
boiling detection, noise-source 
determinations 


13 
’ 


Elk River,? GTR,!® HBWR,? HFIR 


and ORR,‘+!9 HF BR,2° HRE- 


» 21,22 
“» 


HTR, JRR-1, JRR-2,?-4 3 Indian 
Point, Savannah, and Yankee,” 


Kyoto U.,?4 ML-1,”5 Pluto,” 


SPERT-IV,” SRE,2’ Senn,?8 
THOR,” Trino,?° VBWR and Big 


Rock Point,?}! LACBWR*? 
Pseudorandom plant Nerva,' Phoebus-1A,! SNAP 
control 
Intervariable cross- 
correlation and spec- 
tral analysis 


DFR,? DMTR,*? Enrico Ferm 


Pathfinder,® Saxton,? SNAP? 


Intervariable transfer functions 


i,*4 Stability study, intervariable transfer 


functions, noise-source determinations 
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analysis of just the neutron flux—presumably be- 
cause of its being somewhat simpler than other 
categories of experiments listed. The results tabu- 
lated indicate experiment incentives and utilization 
potentials. 

Noise experiments depend on the reactor power and 
other variables being affected by a randomly fluctu- 
ating excitation source. This can be an external 
device (such as a control-rod mechanism) or some- 
thing intrinsic in the reactor (such as water-tem- 
perature fluctuations stemming from inlet flow con- 
ditions). Whether external excitation or self-excitation 
is used to obtain information in a given frequency 
band depends first on the excitation energy available 
in that band and then on other factors, such as those 
in Table 5. 


Table 5 A GENERAL QUALITATIVE COMPARISON 
OF EXCITATION EXPERIMENTS WITH 
INTRINSIC-NOISE-ANALYSIS EXPERIMENTS 








Intrinsic 

Excitation noise 
Experimental complexity and cost More Less 
Interpretation of data Easy Difficult 
Disturbance to reactor system Some None 
Measures transfer function Always Sometimes 
Measures spectra No Yes 
Typical precision High Medium 





Among the external-excitation experiments, the 
pseudorandom has in recent times been preferred 
over the sinusoidal where it is not necessary to have 
a high concentration of spectral energy at discrete 
frequencies. Random excitation supplies a wide band 
of frequencies, all at once, throughout the duration 
of the experiment. Kerlin®* has discussed this energy 
content for various types of excitation. 


Theoretical Analysis Techniques 


It is usual to generate a number of statistical func- 
tions from the noise data to assist in its analysis. 
These range from the simplest—the mean value of 
variables —to rather sophisticated concepts such as 
the coherence in Fig. 1. Some noise experiments are 
engineered so as to obtain transfer functions (i.e., 
variation of system output amplitude with frequency 
for a fixed input amplitude to the system at all fre- 
quencies). In any event, from the magnitude and 
behavior of these quantities, it is possible to deduce 
information about the reactor’s dynamic character- 
istics. Categories of this information in order of 
increasing sophistication are: 

¢ Summary statements or qualities characterizing 
the reactor’s stability. 

* Values of significant parameters (frequencies, 
time constants, coefficients). 

¢ Qualitative and quantitative relations among the 
variables observed. 


Approaches for inferring this physically meaningful 
information are about as numerous as the experi- 
ments in Table 4. The success of a noise test often 
lies in the interpretation of the data, which is espe- 
cially difficult in the self-excited experiments. In 
these, investigators have found the following ap- 
proaches to have merit (in the order of increasing 
analytical effort required): (a) absolute magnitudes 
of variances of the variables, along with their changes 
with experimental conditions (power, flow, control- 
rod configuration, etc.); (6) examination of obvious 
characteristics of spectra and correlation functions 
(such as characteristics of maxima); (c) function 
behavior (such as the frequency, power, and flow 
dependence of spectra) compared to that of postulated 
theoretical models; and (d) simulations (such as 
combining a computer and a mock-up apparatus) of 
the reactor’s noise output for the purpose of match- 
ing measures of performance. 


One feature of noise tests of exceptional value to 
data interpretation should be singled out, for itis a 
rare experiment that has this feature: variety in the 
reactcr conditions used. For example, taking data 
under two significantly different control-rod patterns 
is a simple matter for most reactors and can give 
the control-rod contribution (if the rods happen to be 
shaking randomly) to the intrinsic noise source, but 
very few of the experiments in Table 4 have done this. 


Intrinsic Noise Sources 


In Table 4, types of experimental results obtained 
from noise tests are listed. These were inferred by 
the methods just described, with stability informa- 
tion being the most popular. However, experiments 
that determine the source of the intrinsic reactivity 
excitation noise are important to the fundamentals of 
nonexcited noise analysis as an experimental tool. 
Where this was an experimental result, it is given in 
Table 6. 

Although identification of a noise source is a first 
step, a quantitative understanding of its characteris- 
tics would eventually be desirable. To date, however, 
in virtually all the cases of noise-source identifica- 
tion, there is little quantitative information. One 
should ultimately expect at least empirical charac- 
terizations (much less difficult than a priovi under- 
standings) that would be useful in noise-analyzing 
reactors of a similar size and type. 


The status of understanding noise sources to date 
is nevertheless sufficient to allow five general state- 
ments to be made: 


1. The total reactivity-noise amplitude increases 
only in proportion to VN for N effectively equal 
and independent multiple sources (example: N ran- 
domly shaking control rods or separate inlet tem- 
perature zones).2, More generally, however, the 
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Table 6 IDENTIFIED INTRINSIC REACTIVITY-NOISE SOURCES IN POWER REACTORS 
Type Reactors Noise source 
BWR Dresden" Flow fluctuations outside the vessel 
Elk River Subcore temperature whose hetero- 


Pathfinder® 
VBWR?! 
DMTR®3 
JRR-2? 


DFR’ 
Enrico Fermi*4 


ML-125 


D,O research 
Fast 
Gas cooled 


Graphite CO, Calder Hall” 


H,O pool and research 


Homogeneous HRE-2 (normal flow)?! 
HRE-2 (reversed flow)? 
Liquid-metal compact SNAP-2? 
PWR Savannah 
Saxton? 
Water boiler JRR-V 


HFBR,”° HFIR and ORR," HTR! 


geneity is caused by feedwater 
mixing in the vessel 

Thermalization noise of high steam 
voids near in-core detector 

Flow affecting moderator density 

Inlet coolant temperature 

Coolant flow shaking core 

Inlet sodium temperature 

Inlet sodium temperature and flow 


Temperature of slowly moving water 
moderator 

Control rod and fuel movements 

Control rods and fuel shaken by 
coolant flow 

Fuel heterogeneity in solution 

Fluctuating flow paths of solution 
going through the core 

Coolant flow in turn affected by 
pump’s line-voltage fluctuations 

Control rods shaken by sea 

Nucleate boiling and also coolant 
flow 


Radiolytic gas bubbles 





Square root of the sum of the squares of the im- 
portance-weighted reactivities applies. 

2. The amplitudes of fuel and control-rod motion 
in relation to each other, when random, are dependent 
on mechanical dynamic relations and incident excita- 
tion forces’ frequency spectra. The resulting reac- 
tivity noise then depends on fluxes and geometrical 
relations in addition to these amplitudes. 

3. Imperfect thermal mixing in inlet coolant fluids 
results in temperature heterogeneities that depend 
directly on system temperature differences and in- 
versely on the mixing volume and turbulence. The 
spectral content of the resulting reactivity excitation, 
for frequencies above the reciprocals of mixing 
times and core transit times, is attenuated markedly. 

4. Inlet flow fluctuations become converted to tem- 
perature and reactivity fluctuations®’ by virtue of 
coolant transit dynamics through the core, and the 
noise is proportional to the temperature and steam- 
void coefficients. 

5. In-core detectors have a local thermalization 
noise of neutrons slowing down in a boiling-water 
reactor’s (BWR) fluctuating moderator density in the 
volume near the detector.* 


Useful Experiments 


To enhance the utility of noise applications in power 
reactors, some variety in the types of experiments 
and their accompanying theoretical interpretation is 
desirable —especially including types rarely or never 
performed to date. A few possibilities for experi- 
ments are suggested here, somewhat without regard 
to their relative merit: 
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«Commonly one uses binary pseudorandom se- 
quences as excitation functions, the excited variable 
spending random amounts of time in only two posi- 
tions or values. In ternary sequences, ‘ three values 
are used. Also, binary and/or ternary excitation of 
two or more different variables simultaneously can 
be done. 


¢ Spatial correlations can be made for one vari- 
able or among several variables. The BWR’s show 
diminishing in-core correlations with increasing de- 
tector separation and increasing local turbulences. 
Flow and temperature spatial correlations at the core 
entrance would describe the noise-excitation spec- 
trum in many reactors. 


e The use of coherence in measuring the amount of 
correlation between two random signals at various 
frequencies has been illustrated in Fig. 1. If two or 
more variables (such as flow and temperature) are 
influencing a third (such as the reactor power noise), 
one might take sufficient data to use partial co- 
herences’ to isolate the individual correlations. 


e By analogy to triple-coincidence counting in 
zero-power reactors,*® one might consider merits 
of experimental techniques using it in power reac- 
tors, i.e., on-line probability determinations, where, 
during time intervals A/; for three variables x, y, and 
Ze 


1X 5< MAB) < Xe as yt; + Ty) < Y%, 


and Z,:< 2lt; + 7%;) <Z, 


P,;=0 otherwise 
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Here Ty and 7, are delay times for coincidences of 
these variables with x. Along similar lines, multiple- 
order correlation functions have been defined’ and 
may be measured. 

e On-line measurement and display for the opera- 
tors’? of noise parameters or functions that would 
be most useful in power-plant monitoring are little- 
used applications at present. For example, one might 
monitor for control-rod difficulties if rod shaking is 
responsible for reactor noise. 

* Out-of-pile simulations,” scaled down for econ- 
omy, can be informative, such as in investigations of 
temperature and flow noise phenomena, 

e Monte Carlo “experiments” on computers might 
be valid for some test purposes. Here the computer 
would generate the random input for a dynamically 
simulated system. 

e One might consider applications of random- 
fluctuation analysis in dimensions other than time. 
For example, spatial fluctuations in a fuel rod’s 
gamma activity could be due to random material 
distributions. 


These suggestions by no means exhaust the possi- 
bilities of future noise investigations. One concludes 
that the research potential and benefits therefrom 
are significant. 
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Methods for Studying 
Space-Dependent Reactor Dynamics 


By Chun Hsu and Thomas P. Mulcahey 


To ensure the safety of the public, operating person- 
nel, and reactor power plant components, reactor 
dynamics (i.e., kinetics, stability, and control) has 
been the subject of considerable study in efforts to 
understand and control power reactors. Throughout 
the development process, the most recent tools and 
methods of modern physics and engineering have been 
adopted to analyze and explain observed phenomena 
and to design reactor and plant systems, control sys- 
tems, and safety systems that have maintained safe 
and satisfactory reactor operation. 

In the earlier stages of reactor development, when 
reactor size and power densities were more modest 
than they are today, it was satisfactory to analyze 
reactor dynamics with simpler calculations which in- 
volved only a few delayed-neutron groups and which 
assumed that neutron dynamics were independent of 
location within the reactor core. Although fuel burnup, 
breeding, and the presence and movement of control 
rods were known to introduce space-dependent ef- 
fects on reactor dynamics, for reactors of modest 
size and power their spatial effects usually could be 
calculated sufficiently accurately by methods that 
ignore the spatial effects. Occasionally, an important 
unexplained phenomenon, such as xenon oscillation in 
medium-size thermal reactors, would provide motiva- 
tion for considering the spatial effects. Still, the 
inability to compute spatial effects so that they could 
be compared with calculations that neglect spatial 
effects limited progress toward developing space- 
dependent models and methods. 

However, with plans being made for the develop- 
ment and construction of larger and more powerful ad- 
vanced reactors, itis becoming increasingly important 
that methods for calculating space-dependent effects 
be devised. Impetus for studying spatial effects has 
been gained from the development of methods for 
calculating the existence of potentially hazardous 
conditions in space-dependent reactors. For instance, 
large discrepancies have been reported in reactor 


REACTOR AND FUEL-PROCESSING TECHNOLOGY, Vol. 11, No. 4, Fall 1968 


power transients calculated with and without consider- 
ing spatial effects.! In the cases studied, the magni- 
tude of the power transient calculated by using 
space-dependent considerations can range from three 
or four times greater (in a 60-cm-diameter core) to 
more than several orders of magnitude greater (in 
a 240-cm-diameter core) than that calculated using 
standard nonspatial methods. From this we conclude 
that, for large reactors, the standard nonspatial meth- 
ods for calculating reactor kinetics are inadequate 
to predict some situations. To ensure the safe opera- 
tion of an advanced power reactor throughout its 
life, it is necessary to consider spatial dynamics in 
the design of reactor and plant control and safety 
systems. 


Although most effort to date has been devoted to the 
development and use of methods that exclude the 
space-dependent effects on reactor dynamics, more 
modern space-dependent models and methods are 
being developed to answer design and operational 
questions that arise. When completed, the methods 
will be capable of calculating heat transfer, fluid 
flow, and multiple-energy-group kinetics considera- 
tions. 


To give the reader a general understanding of the 
various methods and some insight into the current 
state of the art, this article discusses briefly the 
formulation of mathematical models and methods that 
are used to analyze reactor dynamics, then outlines 
and evaluates the concepts from the standpoints of 
their bases and application. For completeness, the 
basic mathematical theory of space-dependent reac- 
tor dynamics and the various approximation methods 
for obtaining solutions will be included, even though 
some of these concepts and methods were discussed 
in two earlier survey papers.’*® Also, because re- 
search in other engineering disciplines is related to 
the study of space-dependent reactor systems, some 
of the more important examples are mentioned. 
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Formulating Mathematical Models 


The formulation of a mathematical model of a re- 
actor system begins by describing the dynamic be- 
havior of the neutron distribution in the reactor. For 
power-reactor systems, it usually is acceptable to 
use the multienergy-group diffusion model‘ that has 
been used in many reactor designs and is generally 
considered quite accurate for power reactors. 

In a power reactor the large amount of heat en- 
ergy released by fission in the reactor coreis usually 
removed by coolant flowing through the core. The 
heat-transfer processes can be formulated separately 
by considering that the fission energy is the internal 
heat source and by constructing the equations of heat 
transfer and hydrodynamics. 

Besides the neutron distribution and the heat- 
transfer processes, other dynamic processes in power 
reactors affect the behavior of the reactor system. 
An example is the xenon process in a thermal reac- 
tor. Other processes, such as plutonium production in 
a fast breeder reactor, can be important in various 
reactor systems and can have a similar effect ona 
longer time scale (but may not be as important) as 
the xenon process in a thermal reactor. However, the 
mathematical description can be constructed easily 
from the radioactive-decay equations.°* 

These seemingly independent processes are im- 
plicitly coupled because the parameters inthe various 
processes usually depend on the states of all pro- 
cesses in the system. For example, the nuclear 
parameters are functions of the fuel and coolant 
temperatures due to the density change, Doppler ef- 
fect, etc. This coupling is inherently within the re- 
actor system; therefore control engineers usually 
refer to it as the “internal feedback.” 

The mathematical models thus derived are usually 
adequate for analyzing space—time reactor dynamics 
under normal operating conditions, and so they have 
been used to study space-time reactor systems. °® 
However, when adverse conditions are considered 
(e.g., accident analyses, loss of coolant, or stoppage 
of coolant flow), the model will be different and should 
be treated accordingly. 


Description of Methods 


Although an accurate mathematical description of a 
reactor system can be formulated, the model con- 
sists of one or more partial differential equations 
in space and time variables. Because these equations 
defy direct solution, other methods of obtaining in- 
formation about the dynamic behavior of the model 
are used. The six principal methods described are: 


. Point-reactor kinetics method. 
. Adiabatic method. 

. Modal method. 

. Nodal method. 


wD 


5. Spatial transfer-function method. 
6. Direct method. 


Five of them reduce the system description to a set 
of ordinary differential equations by various approxi- 
mation, linearizing, or averaging techniques. One re- 
tains the original equations and employs comparison 
techniques. 


POINT-MODEL REACTOR KINETICS 


The concept of point-model reactor kinetics can be 
stated simply as the study of reactor dynamic be- 
havior by means of spatially averaged or weighted- 
average quantities of the system variables. In this 
method, as its name implies, the reactor is con- 
sidered as a point rather than as a distributed sys- 
tem in space. It also can be interpreted by saying that 
the spatial distribution of the system variables does 
not change from some estimated shape. Whichever 
interpretation is used, the spatial variations of these 
variables are no longer considered in the analysis, 
and the only concern is the time variation of the 
averaged quantities. 

For deriving point-model reactor equations, the 
various approaches proposed are based on either 
physical or mathematical viewpoints. The basic as- 
sumption common to all these approaches is that the 
space-time solution of the system variables can be 
separated into a single product of a space-dependent 
function and a time-dependent function. A space- 
independent equation is obtained by performing spatial 
integration over the reactor region. The differences 
among the various approaches involve the choices of 
the space-dependent weighting functions. The most 
general discussion on the derivation of the point- 
model kinetics equations is given by Henry,’ who as- 
sumes the weighting function arbitrarily. The method 
reduces the system model to a set of ordinary dif- 
ferential equations. 

The most often used weighting function for the point 
model is the steady-state adjoint solution of the sys- 
tem equations. ® This choice is based on the fact that 
it causes the first-order part of the error in the 
computation of “reactivity” to vanish.’”° The physical 
interpretation is well understood as a result of the 
work by Hurwitz'! and Ussachoff” in the critical 
case, and Selengut'? and Kadomtsev'! in the sub- 
critical case. Other choices of weighting functions 
are compared by Gozani!® and Lewins. '® 


Pros and Cons. The development ofthe point model 
was very useful in the study of reactor dynamics in 
the early stages of reactor development. Its advantages 
were: 

¢The introduction of a concept of reactivity that 
enables reactor dynamics to be studied without know- 
ing details of the physical processes in the reactor. 

eThe simplification of the much more complex 
space—time reactor-system equations so that well- 
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developed methods for studying system stability and 
control can be applied directly to reactor systems. 


However, the basic assumption used in deriving 
the point model means that no knowledge ofthe spatial 
variation is obtainable with this model, so it is not 
suitable for studying large-size power reactors, fast 
or thermal. Thus the consequences of using this 
model are: 

¢ The stability analysis for the reactor might not 
be reliable because the parameters in the point- 
model equation are not reliable; furthermore, there 
is no way to account for the local peaking of the sys- 
tem responses (e.g., temperature) that might jeop- 
ardize the safe operation of the reactor. 

e The control-system design based on the point 
model might not be adequate; in particular, if control 
is to be optimal (e.g., to produce uniform burnup), the 
point model might not even be near true optimal 
without taking into account the spatial effect of the 
control method, 


Numerous analyses of reactor stability and control 
have used the point-model approximation; the major 
references are cited in Ref. 17. 


ADIABATIC METHOD 


The adiabatic method introduced by Henry and 
Curlee’®*”® is essentially based on the same principle 
as the point model. The difference is thatthe adiabatic 
method modifies the solution of the system variables 
as well as the spatial weighting function at various 
time intervals while retaining the point-model formu- 
lation. Lewins?’ and Marchuk and Orlov”! have dis- 
cussed using the time-dependent adjoint solution as 
weighting function; that approach, which also can be 
categorized as the adiabatic method, usually reduces 
the system either to a set of time-variant equations 
or to time-invariant equations between successive 
time intervals. 

This modified version of the point model implies 
that the flux shape responds instantaneously to the 
change of reactor properties. However, its inap- 
plicability for large power reactors means that it has 
the same drawbacks as the point model, although to 
a lesser extent. 


MODAL METHOD 


In concept, this method approximates the space — 
time solution of the system variables by a linear 
combination of some known space functions (modes) 
with undetermined time-dependent coefficients. With 
this method the space—time system equations can be 
reduced to a set of time-dependent equations of the 
modal coefficients by performing spatial integration 
over the reactor. The various types of modal meth- 
ods involve the selection or construction of the space 
functions. The most basic and common functions 
selected for this purpose are an infinite set of 
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orthogonal or biorthogonal functions. Although this 
method is much simpler for most one-dimensional 
problems,”"*3 the general approach is also suitable 
for systems of higher spatial dimensions“**> despite 
the complexity involved. 

Various types of space functions have been proposed 
for the modal method. Generally, they can be cate- 
gorized into two major areas: (a) the set of the 
space functions can .be the eigenfunctions obtained 
from solving the eigenvalue (steady-state) equations 
(sometimes the set of space functions is referred to 
as the natural modes of the system); examples are the 
 modes,”* w, modes, wz modes," 1 modes,”® Green’s 
Function modes, and the inhour modes;”’ or (b) the 
set of the space functions can also be arbitrarily 
chosen, but usually the functions in the selected set 
are mutually orthogonal (e.g., the Helmholtz modes”), 

The procedure of the method is well illustrated by 
Kaplan.*! Applications to the study of stability and 
control of reactors can be found in Refs. 32-37. 


Pros and Cons. The advantages of the modal 
method for analyzing stability and control can be 
summarized as follows: 

¢The method reduces the space-time reactor 
system to a set of time-dependent equations that can 
be analyzed using standard well-developed control 
theory. 

¢ The spatial variation of the system variables can 
be calculated by analyzing the coefficient of each 
spatial mode. 

¢The procedure is straightforward for linear sys- 
tems and, therefore, eliminates the complexity in 
the derivation of the model. 


However, despite these advantages, there are se- 
rious disadvantages when the reactor system is com- 
plex or nonlinear: 

e Except in extremely limited cases, the eigen- 
functions (natural modes) are very difficult to cal- 
culate. 

eIf the space functions are arbitrarily chosen, 
instead of calculating the natural modes, many modes 
will be needed and the resultant time-dependent 
equations will, in general, be coupled; therefore, in 
stability and control analysis, a large number of 
simultaneous equations must be treated — particularly 
when the system is nonlinear, the task may very well 
become insurmountable. 

¢ For a nonlinear system, whether the natural modes 
(for the linear part only) are obtainable or not, the 
resultant equations will generally be coupled; further- 
more, there is no definite knowledge of how many 
modes are required for a specified accuracy. 


In conclusion, the modal method is adequate for 
simple linear space-time reactor systems. For 
more complex or nonlinear systems, the method be- 
comes undesirable because of the amount of effort 
required in the analysis. 
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NODAL METHOD 


The concept of nodal method (coupled core) in- 
volves nothing more than considering a space—time 
reactor system as a finite number of points (nodes) 
or smaller reactor systems coupled by the interac- 
tion between adjacent nodes (e.g., due to flux leakage 
or heat transfer between interface boundaries). The 
procedure is to divide the entire reactor system into 
strategically located smaller regions or zones. The 
same basic assumption used in the point model (and 
adiabatic model) is imposed in each and every region, 
and the final node equations are derived by the same 
procedure as in the point-model formulation. The 
coupling coefficients are obtained by means of the 
boundary conditions given for each region. 

The finite-difference technique*® can be considered 
as a limiting case of the nodal method, This approach 
has been used for analog studies®*»*? and digital- 
computer programs‘! for obtaining the space—time 
solution of reactor systems. There are also cases in 
which some feedback effects were included in the 
analyses,‘°~* Despite the simple concept in this ap- 
proach, it is extremely unsuitable for stability and 
control analysis. The major drawback is due to the 
large number of node equations involved in the 
analyses (e.g., for a simple two-dimensional problem 
with four equations and a 10 by 10 mesh, 400 equa- 
tions must be handled). 

When the nodes are few and represent gross effects 
in each region, leakage parameters must be obtained 
from experiments or from time-independent calcula- 
tions. Two examples are the determination of the 
coupling coefficients (ij) in Avery’s theory of coupled 
reactors’® and the determination of the channel trans- 
fer coefficients in Wachspress’s multichannel flux 
synthesis.‘"»48 The reactivity concept used in the point 
model can also be modified for the nodal method by 
introducing a reactivity to each node. Examples of 
this type of analysis are those by Baldwin, * Danofsky 
and Uhrig®® on Argonaut reactors, and Cohn*! and 
Gamble” on reflected-reactor kinetics. Applications 
of this method to space-time reactor kinetics are 
described in Refs. 53 and 54; more references can 
be found in Ref. 55. 


Pyros and Cons. The advantage ofthe nodal method, 
in addition to those of the point model, for analyzing 
the stability and control of space—time reactor sys- 
tems is that it does take into account some spatial 
effects. Disadvantages of the nodal method can be 
summarized as follows: 

¢The reliability of the analyses depends on the 
complexity of the system as well as on how many 
nodes are taken in formulating the model; in most 
cases, when the nonlinear feedback is to be included 
in the analyses, the number of node equations re- 
quired for the desired accuracy can be very large, 
which is undesirable in the application of present-day 
control theory. 


¢ The determination of the coupling coefficients 
usually poses a problem and, also, it is rather dif- 
ficult to solve a set of spatial nonlinear equations 
whether it is time-dependent or not; to obtain the 
coefficients experimentally is no easy task because 
validly interpreting the data depends strongly on 
knowledge of the system solution. 


TRANSFER-FUNCTION METHOD 


In systems where the linear dynamic processes can 
be described by a set of ordinary differential equa- 
tions, the transfer-function representation of the 
processes is very useful in the sense thatthe stability 
and control analyses can be conducted by studying the 
system parameters without finding the solution of 
system equations. Intuitively, this concept can be 
carried over to the general space—time systems. In 
the study of space—time reactor kinetics, there has 
been much work done in defining the concept of 
space-dependent transfer functions. Examples are: 
the spatial transfer function derived for heat-transfer 
problems with boundary inputs and for boiling-water 
reactors by Sanathanan, °° the transfer function for 
reactor systems using modal expansion by Haugset,” 
the nodal method by Gage,*® and the spatial transfer 
function derived through age equations by Kylstra 
and Uhrig. The most general study on the concept 
of space-dependent transfer functions for reactors is 
by Nomura, who defines the transfer function by the 
concept of generalized Green’s transfer matrix. The 
resultant transfer function depends on the position of 
the detectors as well as the inputs. Other analyses 
including the study of modal interactions in measur- 
ing reactor transfer functions are those by Loewe*! 
and Naqi Zaidi.* Application of these methods to 
stability analyses for reactor systems has been 
studied by Guppy, Takeda,®4*® and Cohn.!% 

In general, because the transfer-function concepts 
only apply to linear systems, most of the analyses 
were derived either by treating linear space—time 
systems or by using perturbation theory to analyze 
system behavior subject only to small disturbances. 
Furthermore, the local disturbances induced by the 
presence of detectors and inputs are not included in 
the analyses. Therefore, if the method is to be useful 
as an experimental tool in the study of stability and 
control of space-time reactor systems, it is neces- 
sary to plan, implement, and interpret meaningful 
experiments. 


Pros and Cons. The advantages of the transfer- 
function method are: 

e Data can be obtained experimentally and, when 
the experiments are repeated at various locations, 
can be helpful in the basic understanding of space— 
time stability and control analysis. 

e The analyses performed on the transfer-function 
model can relate the system stability to the system 
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parameters without the necessity of finding the solu- 
tion of the system. 

e It is equally applicable to the various approaches 
in the study of space—time systems. 


The disadvantage of the method is that it cannot 
handle nonlinear systems directly, unless the analyses 
are carried out for various operating conditions ofthe 
reactor system. Therefore, it is not desirable to use 
this method in stability and control analyses in the 
design stage. However, the method should be suitable 
for predicting system stability from extrapolation of 
experiments made at various operating conditions, 
such as application to existing point-model reactors 
(e.g., EBWR and EBR-]). 


DIRECT METHOD 

This method treats the space-time reactor sys- 
tem directly from the governing partial differential 
equations, rather than reducing the system to a 
lumped-parameter model described by a set of 
ordinary differential equations. This direct approach 
has been studied by many (Wang,® Parks, Hsu,™ 
Kastenberg”) in the fields of power-reactor dynamics, 
aeroelasticity, and related fields. The basic theories 
used in the stability analysis of space-time systems 
are the stability theory of Lyapunov”? and the use of 
comparison functions,"' whereas for control applica- 
tions the theories used are the maximum principle 
of Pontryagin,” variational techniques, and dynamic 
programming.” The application of these theories to 
the general lumped-parameter system as well as 
reactor systems has been very successful; a few 
selected successful applications are given in Refs. 
75-79. The extension of these theories to general 
distributed-parameter systems was carried out by 
Zubov®® and Wang*! 8 for the areas involving stability 
analysis, and by Butkovskii and Lerner®*»* and 
Wang for areas involving optimal control. The ex- 
tended theories have been applied to space-time re- 
actor systems only relatively few times, the only 
publications to date being the study of simple linear 
homogeneous reactor system stability by Lundquist 
and Weissglas, ® the work in the stability analyses of 
general nonlinear reactor systems and optimal con- 
trol of reactor systems by Hsu,** 87 and stability 
analysis using the comparison function by Kasten- 
berg.£9,88 


Pros and Cons. The advantages of this direct 
method are that its use: 

¢Does not require solving the system equations. 

¢ Provides information about the effects of system- 
parameter variations on the system stability. 

eDoes not require engaging in any further ap- 
proximations beyond those in the derivation of the 
system equations. 


The one disadvantage of the direct method in 
stability analysis is that both Lyapunov’s method and 
the comparison function are approaches rather than 
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systematic methods, so there usually is no way to 
develop Lyapunov functionals or comparison functions 
for either the general linear or nonlinear systems. 
However, it is possible, with experience, to obtain 
practical design information for the specific systems 
of interest. 


Research in Other Disciplines 


The study of distributed-parameter dynamic sys- 
tems is not confined to nuclear power-reactor sys- 
tems. Related research is being conducted in other 
disciplines. Many of the developments in treating 
stability and control problems in other fields can be 
used directly or indirectly for analyzing power re- 
actor systems. Some of the more important develop- 
ments in other fields, other than those mentioned 
previously, are discussed briefly here. 

In space-time stability analysis, Rakhmatullina® 
applied the direct method to obtain regions of stability 
for nonlinear heat-conduction equations, Domshlak”’ 
obtained regions of stability for general parabolic 
nonlinear systems, Eckhaus”! discussed stability prob- 
lems using power-series expansions, and Wang” ob- 
tained the sufficient asymptotic stability conditions 
for stable nonlinear parabolic systems with feedback 
control and transport time delays. 

In space—time optimal control, work has increased 
in recent years. Lure” obtained the necessary opti- 
mality conditions for the general distributed system 
using the classical calculus of variations, and Egorov" 
formulated the optimal-problem quasilinear dis- 
tributed systems and obtained the optimal condition 
using Rozonoer’s method.* More recently, Wang and 
Tung” derived a maximum principle for diffusion 
systems in integral form, and Wismer”™ treated the 
optimal problem using multilevel control techniques; 
both papers used spatial discrete approximations. 
Kolb” developed optimal design concepts and tech- 
niques using finite integral transform, which is es- 
sentially one form of the modal method. Uzgiris”® 
obtained the necessary condition for optimality ofa 
linear parabolic equation with nonlinear boundary 
conditions using the gradient technique as applied to 
function space. Axelband’”’ and Yeh and Tou!” treated 
linear distributed systems with bounded controls. 
Axelband used the convex programming algorithm, 
whereas Yeh and Tou used conventional dynamic 
programming techniques. Brogan'™ gives results for 
the linear diffusion system with distributed control 
using dynamic programming and variational tech- 
niques; Sakawa!™ treated a similar problem with 
boundary controls. Briggs and Shen’ treated a 
simplified nuclear rocket system with constraints 
on the system states (temperature and thermal 
stresses). r 

Because most of the referenced developments have 
merits in their respective areas, their pros and cons 
will not be discussed here. 
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Conclusions 


Although the advantages and disadvantages of using 
the various approaches in the control and stability 
problems of spatially dependent power-reactor sys- 
tems have been discussed here in a general way, 
there is no clear indication that any one method is 
distinctly superior for the analysis of power reactors 
in general. The choice of approach for any specific 
reactor depends on many interacting factors. To make 
a reasonable choice, one must first carefully define 
the objectives of his investigation and formulate an 
accurate mathematical model of the system. Then he 
must define the geometrical, physical, and mathemati- 
cal characteristics of the system model and investi- 
gate the available analytical, computational, and ex- 
perimental tools. From this information he then can 
evaluate the various methods with respect to factors 
such as the type of approximation required, the com- 
putational effort involved, and the respective ac- 
curacy, as well as the likely consequences of using 
the various methods. By weighing every aspect of the 
evaluation, he can decide which method would be 
most suitable for the system, as well as the one most 
likely to yield reasonable results that will achieve 
his objective. 

Most of the references cited here are to work 
which is theoretical in nature and which constitutes 
extension of the well-known theories developed for 
lumped-parameter systems. In our opinion, the theo- 
retical investigations have reached a state of maturity. 
Although practical application of these theories is 
still in its infancy because of the computational dif- 
ficulties in soiving the practical problems, the in- 
tense interest among the researchers causes us to 
expect significant progress on practical problems in 
the near future. 
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REACTOR SAFETY 





Coolant Ejection Mechanisms 
for Use in Analyzing Fast Reactor Safety 


By Hans K. Fauske 


During steady-state operation, the coolant in a liquid- 
metal-cooled fast breeder reactor (LMFBR) is con- 
siderably below its boiling temperature throughout 
the core. But abnormal reactor behavior (reactivity 
perturbations and coolant-system failures) might 
greatly increase coolant temperature due to rapid 
heat insertion, a decrease in flow, or both, and might 
initiate boiling. When this occurs, reactivity might 
increase, but the rate of increase depends primarily 
on how rapidly sodium is ejected from fuel-coolant 
channels. This article reviews various models that 
have been proposed to describe liquid ejection from 
reactor coolant channels. A previous article! evalu- 
ated our knowledge of the superheat problem of liquid 
metals, a prerequisite for describing the ejection or 
expulsion of the liquid sodium from the reactor core. 


Rapid Channel Voiding 


Two basicaily different approaches have been 
adopted in theoretical studies of transient void forma- 
tion: (1) vapor growth according to thermodynamic 
equilibrium’ > — zero liquid superheat before as well 
as during the void formation—and (2) vapor growth 
based on nonequilibrium**— liquid superheating 
might occur in the single- and the two-phase regions. 
The various two-phase flow patterns assumed in these 
analyses are summarized in Fig. 1. 


Equilibrium Models 


Models pertaining to equilibrium vapor growth in- 
clude TRANSFUGUE-I (Ref. 2), BURP (Ref. 3), the 
MacFarlane analysis (Ref. 4), and TRANSFUGUE-II 
(Ref. 5). 


TRANSFUGUE-I 


The equations that describe this model’ are: 


Continuity: 
ap , 9 be 
ya + ie (pu) =0 (1) 
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Momentum: 
© a Oe 
Pa ta, ee * oe a (2) 
Energy: 
ah ah 
y(z,t) — p =F — pu aa (3) 


This set of equations applies to a homogeneous two- 
phase mixture where the density is defined by: 





PiPg 
= - 4 
° Pg (1 —x) +p,x (4) 
Because thermodynamic equilibrium is assumed, the 
vapor quality is: 


<a (5) 


EQUILIBRIUM MODELS 
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Fig. 1 Flow regimes 
used in expulsion models. 
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Because numerical instabilities occur in solving these 
equations, the following additional restrictions were 
imposed: (1) constant vapor density and (2) satura- 
tion temperature constant everywhere in the chan- 
nel. These restrictions imply that the calculated 
voiding rate becomes essentially independent of the 
momentum equation and is not modified as a result 
of pressure buildup (increased vapor density and 
saturation temperature). Thus, for the flow regime 
chosen, the TRANSFUGUE-I model will overestimate 
the void growth in the case of zero initial super- 
heat. The pressures calculated from TRANSFUGUE-I 
for a typical loss-of-flow transient are shown in 
Fig. 2 together with sodium vapor pressures cor- 
responding to calculated temperatures (Fig. 3). As 
pointed out in Ref. 2, a considerably larger pressure 
is necessary for this type of void growth (energy 
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Fig. 2 Inertial and vapor pressures resulting froma typical 
loss-of-flow transient, as calculated by TRANSFUGUE-I 
code.3 




















Boiling region —->=——"1 
1900 + _— Saturation 
Temperature profiles y, temperature 
that would exist 

1600 L without boiling / 0.88 sec 
uw Surface at ae 0.86 
© 1499 - (0.88 sec ancy 
2 4 0.84 sec 
s } 0.70 sec 
2 120 + 
= 
2 4 Surface > —_ ae 

Y. 0 sec 
Kaa 3 ite Bulk 
Pr git 
“a 
800 = 
it | | i | i 
1 3 4 6 7 
Axial Distance, ft 
| Blanket | Core | Blanket | Vent | 





Fig. 3 Sodium temperatures calculated to occur during 
typical loss-of-flow transient.3 


controlled) than is available from the only source of 
pressure, namely, the vapor pressure of sodium. 


THE BURP MODEL 


The occurrence of strong self-pressurization ef- 
fects from TRANSFUGUE-I calculations of power and 
loss-of-flow transients and numerical instability 
problems resulted in another model called BURP.’ 
The assumptions used to construct this model are: 
(1) the void spreads axially along the channel from a 
single starting point (Fig. 1) (note that the vapor is 
in contact with the heated surface), (2) the pressure 
in the void region is equal to the highest vapor pres- 
sure of the bounding liquid, (3) the self-pressuriza- 
tion effect suppresses boiling away from the growing 
void region, and (4) the void fraction within the vapor 
region is 100% (Ref. 3). 

Figure 4 illustrates core voiding and pressure 
buildup for the loss-of-flow transient (same condition 
as used in Fig, 2 for TRANSFUGUE-I). Being con- 
sistent with the vapor pressures, the maximum pres- 
sures shown in Fig. 4 are considerably less than 
those calculated from TRANSFUGUE-I. However, be- 
cause energy considerations are not included in the 
analysis, the method used for estimating the bubble 
pressure becomes questionable. Furthermore, if a 
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Fig. 4 Maximum pressure and location of void boundaries 
after boiling in a flow transient as calculated with the 
BURP model.3 


liquid layer is present on the heated wall (most 
likely the case and will be discussed later), the 
BURP model will predict too small a voiding rate, 
since evaporation of this layer will increase the vapor 
pressure over that predicted and will result in larger 
expulsion rates. 


THE MacFARLANE ANALYSIS 


This model* represents improvements over the 
TRANSFUGUE-I code in three areas: (1) the ref- 


REACTOR AND FUEL-PROCESSING TECHNOLOGY, Vol. 11, No. 4, Fall 1968 








182 REACTOR SAFETY 


erence pressure (or saturation temperature) along 
the channel varies according to the steady-state 
profile during the two-phase transient rather than 
being fixed at a constant value, (2) an annular flow 
model is assumed, allowing for slip rather than equal 
phase velocities (Fig. 1), and (3) a space-time 
dependent vapor density (feedback from inertia and 
friction considerations) is used rather than a con- 
stant value throughout the channel. 
The governing equations are: 


Continuity: 


a ra] 
=) [“p, +(1 — @) p,| + [A pgty 


+(l-a)p,u,]=9 (6) 
Momentum: 


) ) 
ot [apg + (1—a@)p,u,] + re fapguz + (1—a)p,u7] 


by or 
+ g lap, + (1- a)p) + — +: 


Energy: 
a 
yo(z,t) - i [apyhg + (1—a)p,h,] 


a ~ {G[xhg +(1—x)m]}\=9 (8) 


The well-known Lockhart—Martinelli analysis’ is 
used for estimating the void fraction for a given 
equilibrium quality (Eq. 5) and liquid—vapor density 
ratio, which represents the equation of state for the 
two-phase mixture. 

The improvements made in this analysis should re- 
sult in a somewhat smaller voiding rate than that 
calculated from the TRANSFUGUE-I model for the 
same initial conditions. However, numerical instabili- 
ties were also encountered in the MacFarlane analy- 
sis; therefore it was necessary to effectively decouple 
Eqs. 7 and 8 to eliminate inertia and friction feedback 
on the saturation temperature. As a result, pressure 
buildup larger than that available from the vapor 
pressure is calculated (see Figs. 5 and 6, which 
represent typical power excursions), again probably 
resulting in a conservative estimate of the core void- 
ing rate for the case of zero initial superheat and the 
flow regime chosen. This model has been incorporated 
in a more complicated calculation that includes the 
reactor kinetics equations and reactivity feedback. '° 

In the MacFarlane analysis, as well as in the 
TRANSFUGUE-I model, two-phase compressibility 
effects were neglected, and the momentum equation 
is used in the so-called “momentum integral” form."! 
This technique implies that pressure is effectively 
eliminated as an unknown in the momentum equation 
and is simply integrated across specified volume 
boundaries. In the case of high-pressure systems 
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Fig. 5 Pressure profiles at various times, as calculated by 
the MacFarlane analysis.!° 
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Fig. 6 Coolant-temperature profiles at various times as 
calculated by the MacFarlane analysis.!° 


(thermal reactors), this technique can be considered 
a valid approximation. 

In typical low-pressure systems (LMFBR’s) with 
large transients, however, the method becomes ques- 
tionable. Since large pressurization effects are likely 
to occur within the coolant channels with correspond- 
ing local momentum changes, these must be accounted 
for in the analysis. 


TRANSFUGUE-II 


This code® is basically the same as the TRANS- 
FUGUE-I code with the following improvements: (1) 
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the elimination of the numerical instability problems, 
making feedback from pressure changes to vapor 
density possible (it is not clear to this reviewer 
whether correction to the vapor quality is made), 
(2) the inclusion of a variable two-phase heat-transfer 
coefficient, and (3) the conversion from FORTRAN IV 
to FORTRAN H. A detailed description of this analy- 
sis, in particular the treatment of two-phase com- 
pressibility effects and numerical techniques used, 
would be welcome. 


Nonequilibrium Models 


As illustrated in Ref. 1, a considerable amount of 
experimental data have become available on liquid- 
metal superheat, clearly indicating that large super- 
heats can occur before boiling is initiated. Models 
taking into account incipient superheat include VOID, ® 
BLOW,’ and Judd’s analyses. ® 


THE VOID MODEL 


This model® represents the first attempt to take into 
account superheating of the liquid metal before boil- 
ing is initiated. For a specified superheat, vapor 
growth and pressure within the channel are computed 
by solving simultaneously the continuity, momentum, 
and energy equations describing a one-dimensional 
bubble-growth model that treats bubbles (assumed 
proportional to the active nucleation site density) as 
slabs of liquid and vapor (see Fig. 1). The heat trans- 
fer across each liquid—vapor interface is calculated 
using an empirical relation between the liquid and the 
vapor temperatures. The strong effect of superheat 
on the voiding rate calculated from this model is 
shown in Fig. 7. The amount of superheat present 
when boiling occurs, together with the rate of super- 
heat release, governs the driving pressure for liquid- 
sodium ejection. These factors also determine whether 
the reactor heat flux is an important parameter in 
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addition, as calculated by the void model.® 


subsequent vapor production leading to complete 
voidage of the coolant channels (Fig. 7). 

The choice of flow pattern used in this model is, 
however, questionable. For example, assuming heat 
transfer to take place only from liquid—vapor inter- 
faces perpendicular to the flow direction (disregard- 
ing the possibility of a liquid layer adjacent to the 
wall) might conceivably lead to underestimation of 
vapor growth and liquid expulsion velocities for a 
given superheat. 


THE BLOW CODE 


Experimental observations of potassium and sodium 
ejections from single channels, ” suggesting that 
the void growth occurs as a Single bubble that grows 
and eventually occupies the entire cross section except 
for a thin liquid film remaining on the heated wall 
(see Fig. 1), led to development of the BLOW code.’ 
It uses the well-known one-dimensional equations 
describing heat conduction in the fuel pin and cladding 
and Newton’s law describing the ejection of liquid. 
The assumed boiling process that supplies vapor to 
the expanding bubble from the liquid surface layer is 
verified by applying the kinetic theory of gases. 
Predicted values of the voiding rates from this model 
were compared with the potassium-ejection experi- 
ments.” Incipient superheat values were calculated 
by matching the ejection behavior. Figure 8 shows 
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Fig. 8 Naturai-convection kinetic behavior during boiling. 
Curves 1 and 2 are experimental data;!* curves 3 and 4 
are calculated with the BLOW code.’ 


ejection starting from natural convection with 65°C 
initial superheat as compared to 29°C measured, and 
Fig. 9 shows ejection initiated from forced convection 
with 25°C superheat as compared to 11°C measured. 
In other words, if measured superheats are used, 
the BLOW code will underestimate the experimental 
ejection rates. However, quoting Ref. 7, “Because of 
the difficulty of measuring the exact temperature at 
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Fig. 9 Forced-convection kinetic behavior during boiling. 
Curves 1 and 2 are experimental data;!* curves 3 and 4 
are calculated with the BLOW code.’ 


the point of bubble nucleation, the discrepancy be- 
tween the measured and the assumed superheat tem- 
perature is not too difficult to explain.” It is perhaps 
possible to explain parts of the deviations as being 
due to oversimplification of the flow regime. As an 
example, if the effects of initial bubble formation and 
liquid film thickness were included in the analysis, 
they might improve the predictions. 


JUDD’S ANALYSES 


Judd® was not sure whether one bubble or many 
bubbles control the expulsion (data on metallic and 
nonmetallic fluids'’-!5 show both flow regimes). He 
therefore examined both cases (Fig. 1). 

His single-bubble solution is quite similar to the 
BURP model,® except for: (1) consideration of initial 
superheating and (2) the bubble growth can be heat- 
transfer as well as inertia limited. (For a similar 
treatment of the bubble-growth problem, see Ref. 14.) 
Again, the shortcoming of this approach, as pointed 
out by Judd, is that of ignoring the possibility of a 
liquid layer in contact with the heated wall. 

In the analysis assuming many bubbles, it was 
found that the number of bubbles taking part in the 
expulsion has a small effect and depends mainly on 
the superheat. The increased surface area due to 
many bubbles leads to expulsion controlled primarily 
by inertia. Again, the latter analysis could lead to 
underestimation of the voiding rate, since the pre- 
dicted driving pressure for the liquid slugs is con- 
siderably less than the vapor pressure corresponding 
to the temperature of the superheated liquid. 


Recommendations for Improvements 


So that critical-heat-flux limits (vapor in con- 
tact with the heated surface), nonequilibrium, and 
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compressibility effects can be evaluated, the basic 
requirement for improving the proposed voiding 
models is reliable information about the flowregimes 
throughout the voiding period. The importance ofcon- 
sidering the two-phase compressibility effects was 
pointed out in Ref. 16, which reported sodium-ejection 
experiments that verified the single-bubble expulsion 
theory (liquid film on the heating surface) but also 
verified that condensation occurs as the bubble or 
vapor slug moves into colder regions of fluid, re- 
sulting in two-phase choking controlling the expul- 
sion rate at the lower liquid—vapor interface. Ref- 
erence 16 also discussed water experiments carried 
out to simulate sodium expulsion; the philosophy can 
be summarized by quoting the paper: “The obtained 
results have shown that there was no basic difference 
between water and sodium vaporization, if similar 
initial conditions are chosen. Particularly the dy- 
namics of ejection is similar and the importance of 
the role played by the critical flow phenomena was 
shown.” 

It appears, therefore, that much of the detailed 
information needed to refine existing voiding models 
can be obtained by using nonmetallic fluids that are 
considerably easier and less expensive to handle. 
This approach was verified by recent experiments ‘7-18 
in which extensive data on flashing flow of sodium and 
water were obtained, suggesting that no major dif- 
ferences existed and that the expulsion or flashing 
leading to critical flow is controlled hydrodynamically 
rather than thermally. However, as yet no satisfactory 
theoretical predictions are available for the two-phase 
choking phenomenon in the low-quality region. 

To evaluate accurately the magnitude of the com- 
pressibility effects requires a good knowledge of how 
pressure waves propagate in two-phase mixtures. Very 
few experimental data are available in this area, '® 
and recent data obtained using two-component mix- 
tures”’ illustrate the strong dependency of flow 
regime on the propagation velocity of pressure waves 
(Fig. 10). These findings emphasize the importance 
of knowing the flow regimes throughout the entire 
voiding period of a reactor coolant channel. 


Conclusions 


From this brief review of proposed models dealing 
with rapid channel voiding, it can be concluded that: 

*In the case of zero initial liquid superheat, 
TRANSFUGUE-I and the MacFarlane analyses will 
generally result in conservative estimates of the 
voiding rates (assuming the flow regimes used are 
correct), with the former predicting the largest 
values, whereas the BURP model will probably under- 
estimate the voiding rate if the heated surfaces re- 
main wetted. It should be noted that experimental 
verification for these models is still lacking. 

* In the case of sizable liquid superheat (the actual 
value of superheat present under reactor conditions 
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Fig. 10 Measured and calculated pressure-wave propaga- 
tion velocities in air—water mixtures. 


remains to be determined, see Ref. 1), most likely 
none of the equilibrium models will yield a dependable 
description of the voiding problem. 

¢ Of the models taking into account liquid super- 
heating, the BLOW code appears the most attractive 
because of its simplicity and extensive experimental 
verification. However, none of these models account 
for the two-phase expulsion that is likely to follow 
the piston-type ejection from a channel blocked at the 
inlet. 


¢ Additional refinements of the models require 
detailed information about the flow regimes through- 
out the voiding process such that heat-transfer 
limitations, nonequilibrium, and compressibility ef- 
fects can be accounted for. 


e At the present time, all the sodium boiling data 
have been obtained in single-channel experiments, and 
these have provided a basis for developing some of 
the one-dimensional expulsion models. Whether these 
models can adequately describe the coolant motion 
in a fuel subassembly depends on results of experi- 
ments with multipin clusters. 


Nomenclature 


ag = pressure-wave propagation velocity in gas 
arp = pressure-wave propagation velocity in two- 
phase fluid 
De = equivalent diameter 
friction factor 
g = acceleration of gravity 
G = total mass flow rate 
fluid enthalpy 
hg = vapor enthalpy 
hi = liquid enthalpy 
htg = latent heat of vaporization 
P = pressure 
t =time 
u = fluid velocity 
ug = vapor velocity 
ui = liquid velocity 
x = vapor quality 
Z = distance 
o(z,¢) = heat flux at channel wall 
a@ = volume fraction of gas or vapor 
p = fluid density 
Pg, = vapor density 
p, = liquid density 
T = frictional pressure 
o7/8z = frictional pressure gradient 
ratio of channel perimeter to cross-sectional 
area 
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NEW CONCEPTS 


Compact Power Concept Features a 
Fast Reactor, Heat Pipes, and Direct Converters 


By John J. Roberts and Edward J. Croke 


A simple, durable, and lightweight power generator 
might result if a compact fast reactor fueled by 
239Dy in a new high-temperature form were com- 
bined with heat pipes, the latest development in 
heat-transfer technology. Described here is a con- 
ceptual design for such a heat source that could take 
advantage of advancements in the technology of ther- 
moelectric or thermionic energy converters. The 
passive self-regulating nature of the system could 
permit unattended operation. A major objective of 
the 1-kw(e) design is the potential for straightfor- 
ward uprating of power to at least 10 kw(e). 

Although the design of the system has not been 
optimized or pointed toward a specific application, 
its possible features and performance are described 
here because they might be of interest in exploring 
the feasibility of various lightweight power sources. 
We believe that reactor systems of this type could 
compete with other power systems where light weight, 
reliability, and reasonable cost are important. 

Fuel Material. The design involves a new and un- 
proved fuel material, plutonium monophosphide (PuP), 
in which 7*°Ppu would be the fissionable constituent. 
Because ***pu is a by-product of the expanding num- 
ber of water-cooled reactors being built by utilities, it 
will be available in useful quantities.’ This fuel should 
cost from 20 to 70% of the equivalent in radioiso- 
topes for power in the 25- to 125-kw(t) range. A 
reasonably firm cost prediction’ of $10/g for **Pu 
implies a fuel inventory cost of about $380,000 per 
reactor, a figure that is essentially independent of 
the thermal-power level. Fuel burnup will be less 
than 0.05% for each 10,000 kw-hr of electricity 
generated. 

Heat Pipes. Heat pipes in both the reactor core 
and the space radiator would offer unique advantages 
for a 1- to 20-kw(e) system that bridges the gap be- 
tween a simple but low-power [0.3 to 0.5 kw(e)] 
conduction-cooled nuclear-thermoelectric genera- 
tor’ and the more complicated but higher-power 


[0.5- to 350-kw(e)] liquid-metal-cooled space power 
thermoelectric or turboelectric systems such as 
SNAP-10A or SNAP-50. Heat pipes used to extract 
heat from the core and reject it to the hot shell, to 
thermoelectric converters, and to space would avoid 
the compromises of reliability and operating life 
which occur with systems that have multiple liquid- 
metal coolant loops and high-temperature rotating 
machinery. Moreover, meteor damage would be less 
drastic in a heat-pipe radiator constructed in modules 
than in a conventional liquid-metal radiator or even 
a vapor-chamber fin radiator.° 

Converters. Recent improvements in thermoelec- 
tric and thermionic devices give these static con- 
verters special appeal for remotely sited electric 
generators. The demonstrated? feasibility of thermo- 
elements having hot-junction temperatures as high 
as 1200°C permits higher radiator temperatures, 
lower system weights, and greater radiation sta- 
bility (due to self-annealing) of thermoelectric ma- 
terials than previously were considered for space 
power generators. °® 

Although the minimum practical emitter tempera- 
ture for thermionic diodes has dropped’ to around 
1200°C, which is only a little higher than the maxi- 
mum temperatures of existing designs for nuclear 
reactors and radioisotopic energy supplies, in the 
1200°C range the state-of-the-art of thermoelectric 
converters is more advanced than it is for therm- 
ionic generators. But the availability of a reliable 
high-temperature source will spur the development 
of thermionic’ devices with probable efficiencies of 
15%. 

Although thermoelectric devices*® are only 4% ef- 
ficient, their availability led us to choose them for 
the reference design described here. 

Weight. In the 1- to 5-kw(e) range, the specific 
weight of radioisotopic thermoelectric power gen- 
erators is typically about 1000 lb/kw(e). As shown 
in Table 1, existing designs for unshielded SNAP -type 
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Table 1 WEIGHTS AND POWERS FOR COMPACT 
POWER GENERATORS 








a Power é Total 
Conversion system Specific 
Thermal, Electric, efficiency, weight, weight, 
kw(t) kw(e) % Ib lb/kw(e) 





Thermoelectric Version of Fast Reactor System 
25 1 4 525 525 
5) 5 4 1700 340 
Thermionic Version of Fast Reactor System 


46 5 15 500 100 
92 10 15 730 73 





thermoelectric generators that employ thermal re- 
actors weigh more than 800 lb/kw(e) in this power 
range. An unshielded SNAP-8 reactor—thermoelectric 
generator weighs 633 lb/kw(e) at 1 kw(e) and 293 
Ib/kw(e) at 20 kw(e).! But our reference fast reactor 
heat-pipe thermoelectric system would weigh only 
about 435 lb/kw(e) for a 1-kw(e) system or 290 
lb/kw(e) for a 5-kw(e) system. Because the reactor 
and heat pipes might supply the converter elements 
with up to 150 kw(t) at 1200°C, more efficient therm- 
ionic diodes could replace the thermoelectric ele- 
ments and produce up to 20 kw(e), at which power the 
total unshielded-system weight would be about 75 
lb/kw(e). (The lithium hydride moderator in SNAP-8 
has a vapor pressure’ of 5 atm at 1100°C, so its use 
with a thermionic converter does not seem feasible.) 


Potential Applications. The use of a ***Pu-based 
fuel in an unusually simple design that employs in- 
core heat pipes and a heat-pipe radiator could pro- 
duce a relatively lightweight and low-cost system 
that should be intrinsically reliable. Thus this fast 
reactor power supply could compare very favorably 
with radioisotopic generators and SNAP reactors in 
the 1- to 20-kw(e) range. 

As interest in communication and other satellites, 
manned space laboratories, and space probes rises, 
a flexible reactor power system that could be adapted 
for varying purposes would be highly desirable. The 
concept described here is intended to be such a 
device. No attempt has been made to optimize the 
design for a specific power or mission. Undersea or 
terrestrial versions of the system would be tailored 
to other means of heat rejection than thermal radia- 
tion. 


System Description 


Heat generated by fission of the ***Pu in each of 
36 unclad plutonium monophosphide® fuel elements 
in the reactor core is conducted (by helium) and 
radiated across a 2-mil gap to a heat pipe that ex- 
tends through the center of the fuel element (Figs. 1 
to 3). This reactor-system heat pipe is a niobium 
tube whose wicks consist of channels in the wall and 
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a niobium screen mesh in which lithium flows as 
the working fluid. Within the heat pipe, the hot 
vaporized lithium travels from the reactor core, 
through the reactor vessel, around two 90° bends, 
and down into a hot shell beneath the reactor. The 
hot shell is a cylinder formed by niobium webs 
welded between the 36 heat pipes. The 1200°C heat 
from the pipes is conducted from the hot shell to 
2222 thermocouples (thermoelectric converters) 
mounted on the outer surface of the shell (Figs. 1 
and 3). The 510°C outer surfaces of the thermoelec- 
tric converters are fastened to the radiator shell. The 
radiator, which consists of 100 rectangular-cross- 
section heat pipes containing potassium as the working 
fluid, rejects the heat at 500°C to space. 
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Fig. 1 Compact fast reacior 1-kw(e) power generator has 
heat pipes to convey heat from fuel elements to hot shell on 
which are mounted thermoelectric converters. Radiator 
has longitudinal heat pipes. 
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Fig. 2 Heat pipe from reactor system to hot shell is a 
niobium tube with wicks in the form of channels in the wall 
and a niobium screen mesh. 


When cooled by the thermocouples, the lithium 
condenses to a liquid and returns to the reactor core 
by capillary action in the wick. 


THE REACTOR 


As shown in Fig. 1, the reactor consists of fuel 
elements, reflectors, reactor vessel, shields, and 
reflector-type control elements. 


Fuel Elements. Each of the 36 unclad fuel ele- 
ments measures 30 cm long, is 2.6 cm across the 
hexagonal flats, and has a central hole 1.12 cm in 
diameter. The fuel material is plutonium monophos- 
phide (PuP), a novel compound® whose thermo- 
dynamic, chemical, and physical properties apparently 
are superior to those of plutonium dioxide (PuO,) 
or plutonium monocarbide (PuC), as is shown in 
Table 2. At the operating temperature the fuel density 
is 8.41 g/cm*®. The 84.9 lb of fuel makes up 67.5 
vol.% of the core. 

Plutonium monophosphide, a dark-gray metallic- 
appearing ceramic, has been prepared in high-purity 
40-g samples by reacting the metal hydride with 
the phosphide.* Because the monophosphide is the 
only stable high-temperature compound in the PuP 
system, the stoichiometric composition was ob- 
tained easily. With a fluid-bed technique, this same 
method can be used for producing the fuel elements. 

In the 1-kw(e) thermoelectric version of the sys- 
tem, the maximum temperature in the fuel is 1269°C, 


which is well below the melting or dissociation tem- 
perature. The PuP compound is extremely stable — 
up to 1500°C in vacuum and 1700°C in 1 atm argon; 
the evaporation and dissociation rates are negligible.° 
At 2600°C the compound melts and decomposes rap- 
idly.2 The linear expansion coefficient (11.9 x 
10-*/°C) and the thermal conductivity [0.09 watt/ 
(cm)(°C)] are nearly invariant with temperature.°*"° 
The PuP fuel elements do not have to be clad be- 
cause of (1) the low fuel burnup (<0.05% at 10,000- 
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Fig. 3 Arrangement of heat pipes from reactor, hot shell, 
thermoelectric converters, and radiator heat pipes. Ref- 
erence design has solid-type hot shell with single web be- 
tween adjacent heat pipes (at upper right). Advanced-design 
hot shell has a double web that forms a vapor chamber 
between the heat pipes (lower right); with wicks and lithium 
in the chambers, the heat and temperature are distributed 
more evenly than with single web. 
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Table 2. PROPERTIES OF CANDIDATE FUEL MATERIALS 
Bare sphere 
Theoretical Melting Thermal Critical Critical 
density, point, conductivity, radius, mass, 
Fuel g/cm? "C watts/(cm)(°C) cm kg Refs. 
Metallic 75U 18.8 1132 0.536 at 900°C 8.4 46 4 
Metallic 233U 18.8 1132 0.536 at 900°C 5.8 15 4 
Metallic 239Pu 19.7 680 0.209 at 600°C 5.0* 10 4 
235UC 12.3 2400 0.251 11.8 84 4 
235U0, 10.0 2750 0.025 at 900°C 14,2 120 4,16 
235 UO, —20 vol. Mo 10.1 0.24 at 870°C 16.27 168 14 
233 Puc 13.6 1654 0.1 at 400°C 7.6t 22 13 
239 PuQ, 10.0 2400 0.025 at 900°C 9.3 34 4 
238 PuC —20 vol.% ZrC 12:42 ~2000 8.9 33 4 
239 puO, — 20 vol.% Mo 9.69 <2000 10.8 51 14 
239puP 9.89 2600 0.09 10.3 45 8, 10 





*Enrichment 96%. 
+Theoretical density 94%. 
{Theoretical density 90%. 


hr design life), (2) the excellent stability of PuP, and 
(3) the fact that the coolant flows within the sealed 
heat pipes rather than in contact with the fuel ele- 
ments, which eliminates corrosion or erosion mecha- 
nisms in the fuel. In addition, because helium ther- 
mally couples the fuel element with the heat pipe, 
the material development problems usually involved 
with high-temperature bonding of cladding and fuel 
are avoided. These characteristics permit a rela- 
tively simple design and facilitate system fabrica- 
tion, assembly, and testing. The high-temperature 
compatibility of PuP and the niobium heat-pipe ma- 
terial remains to be investigated. 


Fuel Selection. The nuclear criteria for the com- 
pact reactor are similar to those proposed by Mac- 
Farlane’ and others for low-power conduction-cooled 
systems; i.e., a fast assembly reflected, if at all, 
for purposes of power flattening or neutron shielding. 
MacFarlane showed that, although some moderated 
cores have lower critical masses than cores that 
have pure fuel, the low melting points of most light- 
weight diluents are an overriding disadvantage. 

To study a single feasible reactor system, we 
selected an optimal reference fuel. Table 2 compares 
the characteristics of available candidate fuel mate- 
rials. Plutonium monophosphide, a novel fuel that is 
exceptionally well suited to high-temperature fast 
reactors, was chosen. 

Because the desired thermal-to-electric converters 
would operate at 1200°C, which means that the fuel 
must be stable up to at least 1400°C, the pure metals 
were unacceptable. Although *U systems probably 
would be cheaper than those with plutonium-based 
fuels, as Table 2 shows, they would be less compact. 
Thus all 7°U compounds were eliminated as candi- 
dates. Uranium-233 compounds were ruled out be- 
cause *y is unlikely to be available in sufficient 
quantities. 

Among the plutonium compounds, PuC is unde- 
sirable because of its relatively low melting point and 
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poor compatibility with most structural materials 
above 1300°C (Ref. 11). Although PuO, has low ther- 
mal conductivity, the maximum fuel temperature 
(~1500°C) at the outer surface of a PuO, fuel ele- 
ment would be considerably below the melting point 
of the fuel; however, at 1500°C, dissociation effects!” 
would compromise fuel-element performance sig- 
nificantly. Moreover, with PuO, fuel, the reactor 
design could not be uprated to 5 to 10 kw(e) because 
the fuel temperature would be excessive. 

The combination PuP—ZrC was considered because 
studies of uranium carbide have shown that the 
melting point of UC increases almost linearly when 
ZrC is added to form a solid solution; although a 
similar effect might be predicted for PuC—ZrC, this 
mixture has not been prepared, '® so it would not be 
an appropriate choice for this design. Similarly a 
PuO,—20 vol.% Mo cermet might have better thermal 
conductivity and stability than plutonium dioxide,'4 
but no data are available for this material. 


Axial Reflector Pieces. To the upper and lower 
ends of each fuel element are fastened 5-cm-high 
beryllium oxide reflector pieces (Fig. 1), whose 
cross-sectional dimensions are the same as those 
for the fuel elements. 


REACTOR VESSEL 

The core containment cylinder is the main struc- 
tural member of the system. Made of Nb—1% Zr, the 
32.8-cm-high vessel has a 30.3-cm outer diameter, 
0.25-cm-thick sidewalls, and 0.5-cm-thick upper and 
lower caps. 

A hexagonal niobium insert within the vessel fits 
around the core and supports and aligns the lower 
ends of the core heat pipes and the fuel elements. 
Also mounted on this insert is a Nb—1% Zr radial 
shield the same height as the fuel elements which 
reduces the radiative heat transfer from the core to 
the reflector and vessel wall. Six ribs protrude 
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radially from the vessel to support the radial reflec- 
tor and the six segments of the hot shell; the axial 
thermal shields and the radiator shell also attach to 
these ribs. 


RADIAL REFLECTOR 


A 5-cm-thick, 30-cm-high, 29.9-cm-OD radial re- 
flector of beryllium oxide surrounds the reactor 
vessel. By reflecting neutrons back into the core, 
this 65.5-lb cylinder reduces the fuel loading (from 
~60 kg for a bare core to 38 kg) and total mass. This 
combination of reflector and core (whose active 
length is 30 cm and critical diameter is 16.6 cm) is 
near the optimum from a weight standpoint. The 
reflector also flattens the core radial power profile, 
which tends to equalize the thermal load on the core 
heat pipes, which in turn ensures a more uniform 
temperature along the circumference of the hot shell. 

The reflector also serves as a radiation shield for 
the thermoelectric converters. 


Nuclear Calculations. The reactor criticality cal- 
culations were performed with the SNARG one- 
dimensional multigroup Sn code using the S4 option. 
Higher-order approximations only negligibly im- 
proved the estimated critical masses or multiplica- 
tion constants. Multigroup diffusion theory does not 
apply because the core radius is less than five mean 
free paths. Hansen and Roach!® 16-group cross sec- 
tions were obtained from ANL Set 201 augmented to 
include *'P, 

Although power peaking could occur near the outer 
edge of the core, this problem is avoided if the re- 
flector is less than 8 cm thick (Fig. 4). Within the 
constraints imposed by requirements for neutron 
shielding and the necessity for power flattening 
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Fig. 4 Power peaking at the core edge is negligible if re- 
Jflector thickness is less than 8 cm. 


without flux peaking, the reflector thickness was 
chosen to minimize reactor weight. The variations 
of reactor height, mass, and 39py inventory with 
reflector thickness are shown in Fig. 5 for the 
actual core configuration. Because this core plan re- 
quired a BeO reflector at least 1.5 cm thick to 
achieve criticality, the weight curve rises sharply 
for thicknesses below 4 cm. The minimum core and 
reflector weight corresponds to a BeO-reflector 
thickness of 5 cm, the choice for this preliminary 
design. 
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Fig. 5 How reflector thickness affects reactor height, 
mass, and 239 Dy inventory. 


The neutron dose rate at the surface of the hot 
shell would not necessarily be reduced by thickening 
the reflector because the power density in the core 
increases in proportion to the core volume reduction 
associated with the reflector saving (Fig. 6). 


Reflector Materials. Table 3 lists some of the 
reflector materials considered for the compact fast 
reactor. Beryllium oxide best satisfies the major 
requirements of a high melting point, a low dissocia- 
tion rate, a low density, and good neutron-scattering 
properties. Graphite is a close second choice. 

Graphite, with thermal conductivity'® as high as 
0.6 watt/(cm)(°C), might be preferred over BeO 
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[0.2 watt/(em)(°C)] for a purely conduction-cooled 
reactor design,’ but this advantage is unimportant for 
a system in which heat pipes provide the primary 
heat-transfer path. 

Beryllium carbide was eliminated because its 
dissociation rate!® is equivalent to 1 x io 
g/(sec)(cm’) at 1200°C, whereas that of BeO is only 
3 x 107' g/(sec)(cm?). The latter figure implies a 
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Fig. 6 Neutron dose to thermocouples during 2.5 * 10° kw- 
hr operation [10,000 hr at 25 kw(t)] as a function of reflec- 
tor thickness. 


loss of only 0.6 g during the 10,000-hr mission life 
of the reactor. 


AXIAL THERMAL-RADIATION SHIELDS 


As shown in Figs. 1 and 3, eachheat pipe undergoes 
two 90° bends between the core (evaporator section) 
and the hot shell (condenser section). In this transi- 
tion region a 4-cm-thick axial thermal-radiation 
shield of Nb—1% Zr laminations is installed to reduce 
the axial radiative heat loss to ~2% of the core 
thermal power. A similar 37.5-cm-diameter shield 
suppresses radiation heat losses at the opposite end 
of the core. The temperature at the outer radiating 
surface of these shields is 364°C. 


STRUCTURE AND RADIAL THERMAL SHIELDING 


The main structural member of the system is the 
core containment vessel. Six equidistant radial ribs 
that protrude from the vessel provide mounting and 
structural support for the reflector segments and for 
the six segments of the hot shell. These ribs also 
provide attachment points for the axial thermal 
shields and the radiator shell. Within the contain- 
ment vessel is a hexagonal niobium insert that fits 
around the core and provides support and alignment 
for the lower ends of the heat pipes. Also mounted 
on this insert is laminar radial shielding that is 
similar in structure to the axial shield. The radial 
shielding minimizes the thermal radiative heat trans- 
fer from the core to the reflector. As shown in 
Fig. 1, this type of thermal shielding is also em- 
ployed to minimize heat loss from the exterior 
surface of the radial reflector. 


NEUTRON SHIELD 


A shadow shield, if needed for the mission, would 
be mounted atop the reactor. Weight and effectiveness 
would be similar to SNAP systems. 


REACTOR CONTROL 


The reactor is designated to operate unattended. 
Both the temperature and power coefficients of re- 


Table 3 PROPERTIES OF CANDIDATE REFLECTOR MATERIALS 








Theoretical Reflector savings Mass of 
density, Melting for PuC sphere with 4-cm reflector, 

Material g/cm? point, °C 4-cm reflector, cm kg Refs. 
Cc 1.6 Sublimes* 1.03 6.0 16, 21 
Be 1.85 1284 1.83 ee | 4, 21 
Be,C 2.4 Dissociatest 2.01 V1 16 
BeO 3.025} 2550 2.02 8.9 16 
LiH 688 4 
ZrH, 1093 4 
Ni 8.9 1455 1,29 31.4 21 
Fe 7.86 1539 1.06 29.2 21 





*P.im = 107 at 1300°C. 
tPam = 4x10 at 1300°C. 


tBut 2.6 g/cm! at the 1200°C operating temperature of the radial reflector. 
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activity are estimated to be negative, which ensures 
stability at all times and enables passive control 
of the system during normal steady-state operation. 
Motion of the control elements for startup and shut- 
down can be preprogrammed to proceed automati- 
cally. 


Reactivity Coefficients. The axialtemperature co- 
efficient (—6 x 10~*/°C) was calculated by expanding 
a 15-cm_ half-slab of homogenized core material 
by a factor (1 + p,,, A7) and diluting the fuel number 
densities by their reciprocal. 

Two thermal-expansion effects of a temperature 
increase in the core influence the radial temperature 
coefficient: (1) the fuel pieces expand into the in- 
terelement voids and (2) the elements move outward 
because they are held by the heat pipes, which in 
turn are positioned by the lower end cap of the 
niobium reactor vessel. The first mechanism cannot 
be incorporated into the nuclear model because this 
effect does not change the homogenized fuel number 
densities of the core material. On the other hand, 
the expansion of the core container, which dis- 
places the elements radially, can be modeled by 
increasing the radii of the fuel annulus. This effect 
results in a radial temperature coefficient of —3.5 x 
10°). 


Startup. To initiate large reactivity changes for 
reactor startup or shutdown, three segments of the 
radial reflector are lowered or raised by small 
stepping motors (Figs. 1 and 7). The movable re- 
flector segments provide the best startup-control 
scheme in terms of reactivity worth, weight, and 
simplicity. Movable neutron-absorber rods or vanes 
were rejected because they would not provide enough 
reactivity control in a design that minimizes system 
weight. 

The number and size of reflector segments that 
would be movable for control purposes were decided 
on the bases of necessary reactivity variation and 
structural considerations. As Fig. 8 shows, the total 
worth of the 5-cm-thick radial reflector is 22.1% 
Ak. Although the curve of differential reflector worth 
shows that a slab of BeO immediately adjacent to the 
core is worth more than one at the outer edge of 
the reflector, practicality required that the control 
segments include the full 5-cm reflector thickness. 
The hexagonal array of the fuel elements and heat 
pipes and the location of internal structural supports 
favored a 60° span for each control segment. 

The reactivity decrement from criticality at zero 
power to full power can be estimated by multiplying 
the temperature coefficient (-4.1 x 10°°/°C) by the 
rise from the 30°C ambient to the 1220°C nominal 
operating temperature. With a shutdown margin of 
6% Ak assumed to be adequate, the required re- 
activity control is 11.1% Ak, which corresponds to 
three full-length (30 cm) 60° reflector segments. 


Reactor startup is accomplished by slowly insert- 
ing two of the three 60° reflector segments, each of 
which is connected to a small stepping motor by a 
screw drive (Fig. 7). When fully inserted, these two 
pieces bring the reactor critical at an equilibrium 
temperature of ~300°C. At this temperature, lithium 
is still frozen in the heat pipes, and heat transfer 
from the core is by pure conduction, so the reactor 
power is very low. Then the third 60° segment is 
inserted, under the feedback control of a temperature 
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Fig. 7 For reactor control during startup and shutdown, 
each of three 60° segments of the BeO radial reflector is 
lowered or raised by stepping motor (2°/step) connected to 
the reflector piece by a screw drive (10/in. pitch). 
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sensor, to bring the reactor up to its 1220°C op- 
erating temperature. 

The entire startup occurs slowly enough for uni- 
form heatup of the system, which passes through a 
series of quasi-equilibrium criticality-temperature 
states during the startup sequence. The desired 
heatup rate is attained by presetting the frequency 
of the pulses that activate the stepping motors. 
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Fig. 8 Reactivity effectiveness as a function of thickness 
of beryllium radial reflector. 


When the system stabilizes at its operating tem- 
perature and power, the startup controller is per- 
manently deactivated. Then steady-state operation is 
maintained by means of passive temperature- 
coefficient control. 


Steady Operation. During steady operation at 
power, the negative temperature coefficients of re- 
activity compensate for the small, gradual reactivity 
loss due to fuel burnup. Any reactivity reduction de- 
creases the power and fuel temperature, promptly 
causing the fuel to compress axially, and, after a 
delay for the lessened heat transmission to cool the 
niobium reactor vessel, thermal contraction causes 
the fuel elements to be displaced radially closer 
together. If the fuel temperature increases for some 
reason, the thermal expansion has the opposite re- 
activity effect and thus reduces the power. 

During the 10,000-hr mission life, 2.5 x 10° kw-hr 
of thermal power would be produced, and 13.2 g of 
*39y (2.19 cm? of the homogenized core material) 
would be burned up. This burnup, less than 0.05%, 
would not affect fuel properties. The reactivity dec- 
rement corresponding to this fuel loss is 2.67 x 
10~* ak, which can be compensated by a 6.5°C drop 
in the core temperature. If the maximum allowable 
temperature decrease during a 10,000-hr mission 
life were 50°C, the maximum reactor power at the 
end of the mission would be 192 kw(t), fuel tem- 
perature permitting. 
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HOT SHELL 


The hot shell, a 44-cm-long cylinder that has a 
33.5-cm outside diameter, is formed by solid niobium 
webs between the condenser sections of the 36 
lithium-bearing heat pipes (see Figs. 1 and 3), Like 
the heat pipes, the webs are Nb—1% Zr alloy. The 
shell is 1.3 cm thick and operates at 1200°C. 

The solid-web hot-shell configuration (Fig. 3) is 
the simplest to design, although not necessarily the 
optimal one. By using a double web to form a vapor 
chamber between adjacent pipes, and by putting 
lithium and wicks within the vapor chamber, we 
could distribute the heat and temperature more 
evenly. The more complex double-web shell also 
would be lighter, but the solid-web hot shell was 
selected for the reference design because of its 
Simplicity and because the solid webs provide mod- 
erately greater shielding for the thermoelements, 


REACTOR HEAT PIPES 


Lithium in 36 heat pipes transfers the heat from 
the reactor core to the hot shell. As shown in Figs. 1 
and 3, the 30-in.-long evaporator section of each 
heat pipe is in the center of a fuel element in the 
reactor. After making two 90° bends below the core, 
each heat pipe extends down into the hot shell. The 
transition section is 23 cm long; the heat-pipe con- 
denser section within the hot shell is 44 cm long. All 
36 heat pipes have 1.111-cm outside diameters 
and 0.076-cm walls and are made of Nb—1% Zr alloy. 
Each pipe is sized to accommodate the 0,88-kw(t) 
peak power generated in each of the six fuel elements 
in the inner ring of the reactor core, although the 
average power per fuel element in the core is 
0.764 kwi(t). 

At full power the lithium in the reactor heat pipes 
operates at 1200°C and 5.3 psia. Along the length of 
the pipe, the vapor pressure and temperature drop 
0.0235 psi and 0.655°C in passing from the evapo- 
rator to the condenser. 

By capillary action, wicks within the heat pipes 
return the liquid lithium from the condenser section 
in the hot shell to the evaporator section within the 
reactor core; the liquid pressure drops 0.1544 psi 
in the return. There are two types of wicks in each 
heat pipe: channels and screen mesh. The channels, 
of which there are 51 equally spaced around the 
inside surface of the 0.762-mm-thick pipe wall 
(Fig. 2), extend the length of the pipe and are 0.51 
mm deep and 0.204 mm wide. To prevent the high- 
velocity lithium vapor from interacting with the 
surface of the liquid lithium in the open channels 
during startup from low temperatures, which would 
tend to retard liquid flow and prevent startup,'" 
there is a single layer of Nb—1% Zr screen (Fig. 2) 
against the inner surface of the pipe. This 0.22-mm- 
thick screen has an equivalent pore radius of 0.08 
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mm, a permeability of 12, and a porosity (void frac- 
tion) of 0.5. 


Working Fluid. The reactor-system heat pipes 
operate at 1200°C to heat the thermoelement hot 
junctions to that temperature. Lithium was selected 
as the optimum heat-pipe fluid for 1200°C operation 
because: (1) its sufficiently high vapor pressure 
above 900°C, its high surface tension, and its low 
liquid density make it an excellent fluid for use in 
an efficient capillary pump; (2) its vapor pressure 
at 1200°C is low (~5 psi), so the pipe wall can be 
thin enough to give both low weight and low tempera- 
ture drop across the pipe walls; (3) the high latent 
heat of vaporization of lithium permits a small pipe 
diameter, which, among other benefits, reduces po- 
tential startup instability problems; and (4) the great 
heat-transfer capacity of a lithium pipe lends itself 
to a system design that can be uprated to higher 
powers with minimal modification of the heat pipes. 


Pipe Material. From the limited number of re- 
fractory metals compatible with lithium at 1200°C 
(niobium, tantalum, and molybdenum and their alloys), 
Nb—1% Zr was chosen for the heat pipe and wick. 
This selection was based primarily on heat-pipe 
tests at the Los Alamos Scientific Laboratory. '® 
Those tests verified the compatibility of the Li—Nb 
combination at temperatures in excess of 1100°C 
for operating times up to 3500 hr, including a num- 
ber of restart cycles. Lithium and niobium are com- 
patible at temperatures as high as 1150°C for up to 
6500 hr (Ref. 19). The LASL tests included runs at 
temperatures as high as 1300°C; however, 10,000-hr 
compatibility remains to be verified experimentally. 


Sizing the Heat Pipes. The excellent thermal 
properties of lithium make it easy for small-diameter 
heat pipes and wicks to carry the design power; 
they also keep the heat flux across the pipe surface 
low enough to preclude boiling within the wick. 

It is also important that the pipe diameter be large 
enough to: (1) keep the PuP fuel temperature below 
1700°C, (2) ensure a near-uniform temperature dis- 
tribution on the hot-shell surface that contacts the 
thermocouples, and (3) permit the system to be up- 
rated from its initial design power of 27.5 kw(t) to at 
least 125 kw(t) without significantly modifying the 
reactor design. 

The fuel-temperature limitation is not a problem 
with the 36-fuel-element configuration because the 
6 elements in the central ring do not get hotter than 
1269°C when the reactor operates at 25 kw(t). Fig- 
ure 9 shows how the peak fuel temperature at the 
core midplane would vary with reactor power for 
fuel elements of five sizes. Note that, for the 2.6- 
cm across-flats hexagonal fuel elements, the 1700°C 
limiting temperature for the PuP fuel would be 
reached at a thermoelectric power of 7.2 kw(e) 
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Fig. 9 Maximum fuel temperature reached in fuel ele- 
ments of five sizes as a function of system electricity gen- 
eration (for 4% efficient thermoelectric conversion). 


[180 kw(t)]. Figures 10 and 11 show the axial and 
radial distributions of temperature in the core. 

The governing criterion for the size of the heat 
pipes proved to be that associated with temperature 
uniformity on the surface of the hot shell on which 
the converters are mounted. A 20°C maximum allow- 
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Fig. 10 Axial distribution of temperature in outer edge of 
fuei element in center ring of core when the reactor is op- 
erating at 27.5 kwi(t). 
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Fig. 11 Radial distribution of temperature in fuel ele- 
ments along cove midplane when the reactor is operating 
at 27.5 kw(t). 


able circumferential variation of the 1200°C tem- 
perature was established (arbitrarily) for the solid- 
single-web shell. With the shell diameter the same 
as that of the 29.9-cm-OD reactor radial reflector, 
and with the number of heat pipes established as 36 
by the number of fuel elements, the pipe diameter 
determined the thickness and circumferential width 
of the web joining adjacent heat pipes. Onthe premise 
that the inner surface behaves adiabatically, the 
smallest heat-pipe diameter (outer) that would limit 
the circumferential variation to 20°C was deter- 
mined to be 1.111 cm (Ae in.). With this dimension 
established, the heat pipe easily met the other re- 
quirements, which were considerably less stringent. 

The theory of heat pipes in Ref. 20 considers heat- 
pipe equations in which the liquid and vapor in the 
heat pipe flow either adiabatically or under conditions 
of uniform heat input and extraction. In the compact- 
reactor core, however, with a “chopped cosine” axial 
power distribution, this assumption is not valid for 
the evaporator section of the lithium pipe, although 
it does apply to the condenser section. For design of 
the in-core heat pipe, the liquid and vapor pressure- 
drop equations were integrated over the evaporator 
length on the basis of a chopped-cosine power profile 
using the integrated average radial Reynolds num- 
ber, which approximates the actual case more closely 
than does the assumption of a uniform power input. 
The condenser section was sized on the assumption 
of uniform power extraction. 

The channel-type wicks in each heat pipe are sized 
to carry the full 0.88-kw heat load of a central fuel 
element. However, an additional ~148 watts of heat- 
transfer capacity is provided by the single-layer 
screen-type wick within the pipe. But the screen is 
counted on only to serve as a barrier between the 
liquid and the vapor to prevent instabilities during 
the low-temperature portion of startup, even though 
a screen wick is superior to a channel wick as a 
capillary pump. Because the radius of curvature of 
the liquid meniscus is effectively infinite in the axial 
direction for a channel wick, the capillary forces are 
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approximately half those in a screen-mesh wick. 
However, unlike the screen, under no circumstances 
can the channels milled in the pipe wall detach from 
the wall and cause hot spots. The channel wick also 
can be designed more readily than a screen wick be- 
cause the channel flow characteristics (flow area, 
hydraulic diameter, etc.) can be defined more ac- 
curately. 

The 0.102-mm effective capillary pore radius of 
the channel is sufficient to wick the liquid lithium 
to a height of 45 cm against earth gravity. This 
wicking height corresponds to the length of the con- 
denser leg of the heat pipe so the core-heat-pipe 
subassembly can function in earth gravity for ground- 
testing the space system. It also means that the 
system can be adapted readily for terrestrial or 
undersea use. 


THERMOELECTRIC CONVERTERS 


The 4/,-in, -diameter thermoelectric converters, of 
which there are 2222, are mounted on the outer 
surface of the 1200°C hot shell (Figs. 1 and 3). These 
1200°C thermocouples are of the type developed by 
Monsanto Research Corp.* With its 510°C cold junc- 
tion at the radiator surface, each of these 4% efficient 
converters produces 0.495 watt of 40-volt electricity 
when the reactor power is 25 kwi(t)., 


RADIATOR 


The radiator that rejects heat at 500°C to space 
consists of 100 rectangular cross-section heat pipes 
that form a 113.5-cm-long cylindrical shell around 
the thermoelectric converters and reactor (Figs. 1 
and 3), Although the full 113.5-cm length of the 
outer wall serves as the condenser section, only the 
44-cm-long portion of the inside wall in contact with 
the thermocouples functions as the evaporator section. 

The gross dimensions of the radiator were dictated 
by (1) the 39.5-cm outside diameter of the con- 
verter array, (2) the 510°C cold-junction tempera- 
ture of the reference thermocouples, and (3) the de- 
sire to reject 26.4 kw(t) at 500°C with an effective 
emissivity of 0.85 from a 4% efficient thermoelectric 
system. 


Radiator Heat Pipes. Each heat pipe has a 1.35- 
by 2.0-cm cross section, is made of Inconel tubing 
having a 0.051-cm (20-mil) wall, contains potassium 
as the working fluid, and has screen-mesh-type wicks 
of Inconel (Fig. 12). Inconel was selected mainly be- 
cause of its long-term compatibility with potassium, 
but also because of its reasonably good resistance to 
penetration by meteors (Fig. 13). The screen-mesh 
wicks are 0.314 cm thick along the inner and outer 
walls of the radiator rectangular pipes; the wicks 
on the lateral walls are 0.157 cm thick. Wick porosity 
is 0.5, pore radius is 0.025 cm, and permeability is 
12. With each of the 100 radiator heat pipes trans- 
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Fig. 12 Cutaway of radiator heat pipes, showing screen- 
mesh-type wicks. 


mitting 0.264 kw(t), the operating pressure is 0.5 
psia, the vapor pressure and temperature drop 
are 0.005 psia and 2.93°C, the liquid pressure drop is 
0.005 psia, and the total temperature drop across 
the radiator wall is 11.9°C. 


Fluid Flow. The configuration and operation ofthe 
radiator heat pipes are unusual in several respects. 
The heat pipes are double-ended, so heat is added in 
the central portion of the pipe, and the vapor flows 
axially in two directions as well as transversely 
(Figs. 1 and 3). The transverse flows of vapor and 
liquid within the rectangular pipe are considerably 
more complex and less symmetric than they would 
be in a circular pipe. The vapor produced in the 
evaporator section on the inner radius of the radia- 
tor shell condenses on the entire outer radius and 
flows back to the evaporator section through the 
upper, lower, and lateral wicks. No satisfactory 
analytical description has been developed for this 
complicated flow pattern. For this study the basic 
heat-pipe theory of Ref. 20 was adapted to determine 
the sizes of the wick and vapor passages. 

Unlike the primary heat pipes in the reactor and 
hot shell, in which the excellent wicking properties 
of lithium permit the system to be tested in an 
“upright” position, as shown in Fig. 1, ground tests 
of the radiator heat pipes (and thus of the overall 
integrated system) have to be conducted with the 
system in a horizontal attitude. For the radiator 
heat pipes to be tested vertically, the radiator-wick 
pore structure would have to lift the potassium 56,75 
cm (half the radiator height) to ensure that the lower 
half of the evaporator would be wetted. To do this, 
the wick pore radius would have to be 0.04 mm, 
which is extremely small. Even if this pore radius 
could be attained in practice, a thick wick in a large 
pipe with a large ratio of width to length would be 


required to keep the liquid pressure drop well 
below the potassium vapor pressure (~0.5 psi at 
500°C). However, with a pore radius of 0.025 cm, 
the wicking height is ~8 cm, which results ina 
design for the radiator heat pipe that has 1.34- by 
2-cm cross section; so ample margin exists for 
Capillary pumping in a horizontal attitude ground 
test of the integrated system. 


Radiator Fluid. The selection of an optimal fluid 
for the radiator heat pipes was less obvious than it 
was for the heat pipes in the reactor and hot shell. 
Lithium is not suitable for the 500°C radiator be- 
cause its vapor pressure is too low below 800 to 
900°C for an effective heat pipe. So the choice rested 
among: (1) the lower-temperature liquid metals, 
such as Na, K, Cs, and Rb (Hg was ruled out because 
of its questionable wetting properties); (2) the higher- 
temperature organic materials, such as the diphenyls 
(e.g., Dowtherm); and (3) the molten salts. 
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Fig. 13 Probable incidence of meteor damage to heat 

pipes of four materials as a function of wall thickness. 
Correlation of Ref. 22 was used to adapt data of Ref. 23 to 
materials of interest for this study. 


To aid in the selection of the best heat-pipe fluid, 
figures of merit were estimated on the basis of the 
relations in Ref. 20 for the limiting heat-transfer 
capacity of an optimized heat pipe operating with 
uniform power input and extraction in a 0-g environ- 
ment. This figure of merit for an operating mode 
with a high radial Reynolds number is 


Fy = L(p, py? /u,)" 
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and for a low radial Reynolds number is 
F, =Ly(p,p,)/(byt))” 


where L = latent heat of vaporization 
p = density 
U. = viscosity 

= surface tension 


~~ 


v, 1 = vapor and liquid 


As shown in Fig. 14, Na, K, Rb, and Cs are ap- 
proximately competitive in the ~500°F range of 
interest. Potassium was selected because of its 
better performance during startup from low tem- 
peratures, although above 500°F sodium clearly 
is superior to the other liquid metals. Cesium was 
eliminated because its relatively low surface ten- 
sion makes it less effective than potassium as a 
capillary pump. 

Heat-pipe experiments at LASL have demonstrated 
the performance of sodium and potassium. '" 
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Fig. 14 Selection of potassium as the reactor heat-pipe 

fluid was made partly on the basis of a figure of merit (Fy) 


that assumes uniform power input and extraction in a 0-g 
environment. 
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Meteor Invulnerability. Relative invulnerability to 
damage by meteors is a major advantage of a heat- 
pipe radiator constructed of an aggregate of inde- 
pendent cells. This concept involves 100 independent 
cells, so the puncture of any one heat pipe would re- 
duce the radiator capacity by less than 1%. 

Although meteor-penetration data are not yet fully 
reliable, Fig. 13 indicates that our system would 
encounter a meteor sufficiently large to puncture 
the 20-mil-wall Inconel heat-pipe radiator approxi- 
mately once in 2 years. Inasmuch as the radiator 
has been designed oversize by 10%, a reasonable 
attrition of heat pipes can be tolerated, so no meteor 
armor should be required. In contrast, although an 
advanced-concept vapor-chamber® finned Rankine- 
cycle radiator also would offer some invulnerability 
to meteor damage, it would require some meteor 
armor to protect the liquid-metal passages from a 
catastrophic puncture, 


Assembly and Testing 


The mechanical design of the reactor system en- 
ables straightforward assembly and handling. After 
the 36 reactor heat pipes are installed in the vessel 
(Figs. 1, 3, and 15) and welded to the flanges on the 
niobium heat-pipe bellows, the core is assembled in 
the vessel by sliding the fuel hexagons over the heat 
pipes. The central position, in what is basically a 
37 fuel-element array, is left vacant to simplify the 
routing of the heat pipes to the hot shell. 

After the core is assembled, the reactor vessel 
is purged and welded gastight under 1 atm helium. 
Then the three 60° movable reflector segments that 
constitute the startup control system are installed 
as three separate modules (Fig. 7), each of which 
includes a reflector segment, a stepping motor, a 
driving motor, and a covering shroud. To prevent 
unintentional startup of the reactor during shipment, 
prelaunch activities, and before the attainment of 
orbit, each reflector segment is held in the with- 
drawn position by a scram spring. During operation, 
the scram spring is restrained by a magnetic latch. 
With the three stationary reflector segments in 
place but the movable segments withdrawn, the 
shutdown reactivity margin is 6% If the three sta- 
tionary segments also are removed during shipment 
and handling, the bare core has a shutdown safety 
margin of 20%. 

For component development and system tests, the 
nuclear and nonnuclear parts of the system are 
separable. Electric heating can be used to test the 
assembled system without the fuel elements. 

Because the fuel is not bonded to the heat pipes 
in the reactor (heat transfer across the 2-mil gap 
between the fuel and pipe is by radiation and, mostly, 
conduction through the helium), no problem is posed 
by the differential thermal expansion between the PuP 
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fuel and niobium pipe during startup. With the nio- 
bium bellows welded on each heat pipe where it 
emerges from the reactor vessel, thermal stresses 
caused by differential expansion do not compromise 
the integrity of the pipes or vessel. 


Matrix hold-down 


Fuel position matrix 


Leaf spring with retainer 


Reactor vessel 


Fuel element 


Radial thermal shields and fuel support 


Axial reflector element 


Leaf spring and pressure plate 


Compression spring 


Bellows assembly 








oS Heat pipe 
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Fig. 15 Alignment and support system for fuel element. 


Optimization and Uprating 


As for any conceptual design, the design parame- 
ters and performance characteristics described for 
this 1-kw(e) compact fast reactor power generator 
are subject to considerable refinement. Nevertheless, 
the predicted weight (Table 1) of less than 450 lb 
for a 4% efficient 1-kw(e) thermoelectric generator 
indicates the performance potential of this concept. 
Doubtless, further design and development studies 
would modify the system performance and weights 
somewhat. 

If the conceptual design were carried to the 
preliminary-design stage, it would be worthwhile to 
optimize the complete system. An obvious improve- 
ment would result from a combined core, reflector, 
and reactor heat-pipe optimization. Desirable opti- 


mization steps such as power flattening in the core 
and improved mechanical design would have to be 
balanced against material temperature limits and 
the ability to increase the design power. Further- 
more, design changes such as the use of smaller 
heat pipes and a modified hot-shell structure would 
significantly reduce the specific weight, particularly 
for a high-power [5 kw(e)] version. 

The power of the 1-kw(e) thermoelectric system 
can be uprated to 5 kw(e) without significantly alter- 
ing the reactor design, although the hot shell and 
radiator would be enlarged. As shown in Table 1, 
the specific weight of the 5-kw(e) thermoelectric 
generator system would be reduced to ~290 lb/kw(e). 
The maximum fuel temperature associated with a 
5-kw(e) output would be 1590°C [vs. 1269°C for 
1-kw(e) operation], which is still within the range of 
stable performance for plutonium monophosphide 
under 1 atm helium.’ Other system temperatures 
would be essentially unchanged; passive tempera- 
ture-coefficient control would still be feasible. The 
reactor heat pipes in the 1-kw(e) system have been 
designed oversize to accommodate the 5-kw(e) power. 


THERMIONIC VERSION 


The same optimization steps would apply to a 
thermionic version of the system (the same hot shell 
would serve as the emitter heat source), where 
reactor-system performance improvements would 
be even more significant because the superior effi- 
ciency of the thermionic system (15% vs. 4% for the 
thermoelectric) would reduce system specific weight 
to about 100 lb/kw(e) for a 4.25-kw(e)-output gen- 
erator. If the reactor were used with a 1-kw/(e) 
15% efficient thermionic generator, the radiator 
length could be reduced from 113.5 cm to ~30 cm. 


IMPROVED HOT SHELL 


The solid-web hot-shell configuration (Fig. 3) 
selected for the reference design is the simplest but 
not necessarily the optimal design for the hot shell. 
As shown at lower right in Fig. 3, the vapor-chamber- 
web concept would be employed to achieve a virtually 
uniform temperature along the hot-shell surface. Then 
the hot-shell configuration would be somewhat more 
complex, in that the web structure between the heat 
pipes would contain lithium and would be lined with 
wicking material similar to that in the pipes from 
the reactor. This aggregate of vapor-chamber webs 
would weigh ~12 kg less than the solid shell. Be- 
cause approximately the same wicking height (44 cm) 
is involved, the wicks in the vapor-chamber webs 
would have the same pore size as the wicks in the 
reactor pipes, but the liquid and vapor pressure 
drops associated with the web would be considerably 
less because the distance over which heat is trans- 
ferred is much smaller (~0.5 cm). 
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We gratefully acknowledge the assistance of Ralph 
Carter, who solved many of the mechanical design 
problems; Jay Norco, who aided in the nuclear and 
thermal calculations; and Owen Kruger, who advised 
on reactor fuels technology. 

This article is an adaptation of a paper presented 
at the Intersociety Energy Conversion Engineering 
Conference, which was held in Miami, Fla., in 
August 1967. Details of a similar space power sys- 
tem concept are also provided in Argonne National 
Laboratory Report ANL-7422, A Heat-Pipe-Cooled 


Fast-Reactor Space Power Supply. 
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AQUEOUS PROCESSING 


Research and Development 
on Aqueous Processing 


By C. E. Stevenson and D. M. Paige 


Silicon in U—Al alloy has little effect on rate of dis- 
solution, but its presence in suspended solids may 
result in surface activity in solutions of such alloys. 
Tests are being made to determine conditions for the 
electrolytic dissolution of sodium-bonded EBR-II 
fuel elements, Semicontinuous dissolution of Zircaloy 
fuels in HF at the Idaho Chemical Processing Plant 
(ICPP) was not accompanied by excessive corrosion. 
Methods are discussed for the recovery of americium, 
curium, and neptunium from Shippingport reactor 
fuel. The radiation deterioration of extraction sol- 
vents is reviewed, and results of initial operation of 
the Eurochemic plant are presented. 

The performance of a centrifugal solvent extraction 
contactor, which has been operated with radioactive 
fuel solutions for 1 year, is discussed. A monitor for 
detecting traces of fissionable material left in the 
hulls from shear—leach processing of tubular-clad 
ceramic fuels is described. 


Process Methods 


EFFECT OF SILICON UPON THE REPROCESSING 
OF U—Al ALLOY FUELS 


In a continuation of earlier studies’ at the ICPP 
upon the effect of minor metallic constituents on the 
rate of dissolution of U—Al alloys in HNO;, Paige and 
Rohde? recently examined the metallurgical structure 
of both aluminum and 22 to 32% U—Al alloys that con- 
tain up to 11% silicon, the dissolution rates of such al- 
loys in boiling 6.6M HNO; containing 0.007M Hg(NOs3)., 
and the surface activity of the suspended siliceous 
solids resulting from dissolution of the alloy. Silicon 
is added in the fabrication of U—Al cores for plate- 
type fuel elements, such as those for the Engineering 
Test Reactor (ETR), to stabilize the UAl; phase and 
to prevent the formation of UAl,, which is present as 
hard, brittle intermetallic particles. 

Dissolution rates of the alloys were found to be re- 
lated to the metallurgical composition, i.e., the phases 


present. The observed rates decreased in the follow- 
ing order: Al >UAI; > U(Al,Si)3 > UAL, > Si. All al- 
loys tested, however, dissolved at rates that were 
practical for processing. Some tendency for the pref- 
erential dissolution of aluminum was noted, so that 
the possibility of enrichment of undissolved residues 
in uranium must be considered. Silicon present in the 
alloy as the intermetallic U(A1,Si); dissolves, but sili- 
con present in excess of one atom per atom of ura- 
nium apparently is not contained in the intermetallic 
and does not dissolve. This undissolved residue may 
require removal during reprocessing. Because the 
presence of even a small amount of silicon results in 
sufficient surface activity to require treatment (e.g., 
with gelatin) to obtain good extraction behavior, the 
concentration of silicon in the alloy is not important 
from this standpoint. Surface activity appears to be 
due to the presence either of polysilicic acid or of 
colloidal silica in the form of particles less than 35 py 
in diameter. 


AQUEOUS PROCESSING OF EBR-II FUEL 


The fuel used in Experimental Breeder Reactor II 
(EBR-II) at the National Reactor Testing Station in 
Idaho is highly enriched uranium alloy which is clad 
in stainless steel and which is bonded to the stainless 
steel with sodium. For the past several years, all fuel 
from this reactor has been reprocessed remotely by 
a pyrochemical process and refabricated into new fuel 
elements at the Fuel Cycle Facility immediately ad- 
jacent to the reactor.’ Small quantities of scrap from 
this process have been recovered at the ICPP.‘ To 
reduce the dependence of the reactor upon the contin- 
ued operation of the Fuel Cycle Facility as its fuel 
source, and to make additional hot-cell space avail- 
able for the handling of materials irradiated experi- 
mentally in EBR-II, it is now planned to provide for 
aqueous processing at the ICPP of all fuel discharged 
by EBR-II. The proposed process method includes 
electrolytic dissolution in HNO3, followed by TBP and 
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hexone extraction. Although extensive studies of the 
electrolytic dissolution of fuels of various types have 
been made at the ICPP,* new equipment will be re- 
quired for this purpose, and pilot-scale studies are 
now under way. ° 

An EBR-II fuel element consists of a uranium-alloy 
rod, 0,144 in. in diameter by 13.5in. long, contained in 
a cladding tube of 304L stainless steel with a 9-mil- 
thick wall; a 6-mil annulus containing a sodium bond 
for heat-transfer purposes separates the rod from 
the cladding. The uranium alloy is described as “fis- 
sium” and consists of 95 wt.% uranium at an enrich- 
ment of 52.2% with the balance being 2.5 wt.% Mo, 
1,9 wt.% Ru, 0.3 wt.% Rh, 0.2 wt.% Pd, 0.1 wt.% Zr, 
and 0.01 wt.% Nb. The latter elements represent me- 
tallic fission products not removed by the pyrochemi- 
cal process. Individual fuel elements are assembled 
into bundles of 91 with a hexagonal sheath and other 
structural components of stainless steel. Because of 
the sodium bond in the fuel element, the safety as- 
pects of the reaction of the dilute HNO; dissolution 
reagent with sodium were investigated. In the disso- 
lution of single fuel elements (less than 1 g of sodium 
in each), the rate of the sodium reaction was suffi- 
ciently low that no safety problem was apparent. Con- 
sideration has also been given to the possible corro- 
sion of the cladding and reaction of sodium during 
extended underwater storage. However, five fully ir- 
radiated fuel pins (~ 1 at.% burnup) have been retained 
for 10 months in the ICPP storage basin with no indi- 
cation of penetration of the cladding. 

In pilot-plant studies, two 10-day continuous dis- 
solving cycles have been carried out in which about 
1200 lb of stainless-steel tubing was dissolved (no 
sodium or fuel was present) using an insulated nio- 
bium basket in a 6- by 18- by 24-in. titanium vessel 
(the cathode), and a platinum-coated niobium anode, 
An average dissolution rate of 1.3 kg/hr was achieved 
with a current efficiency of 0.5 g/amp-hr at varying 
acid concentrations, The principal problem encoun- 
tered was the accumulation of a heavy sludge of par- 
ticles of iron, chromium, and nickel bound together 
with silicic acid or silica gel. The addition of gelatin 
to the dissolution reagent did not reduce the sludge 
accumulation. Further studies will include evaluation 
of the safety aspects of the sodium reaction and in- 
vestigation of methods for reducing or disposing of 
the sludge. 


CORROSION IN PROCESSING ZIRCONIUM-T YPE 
FUELS WITH HF 


The hydrofluoric acid dissolution process, used in 
the processing of Zircaloy-clad and -alloyed enriched 
uranium fuels at ICPP, requires the use of vessel 
materials of limited corrosion resistance to the com- 
bination of reagents that must be used at various steps 
in the process, For example, aqueous HF, to which 
relatively few alloys are resistant, is required for 
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dissolution of zirconium. Other reagents include 
small amounts of HNO3, which were used originally 
in the batch process to oxidize uranium and facilitate 
its dissolution, and Al(NO3)3 solution containing 0.02M 
chromic acid, the complexing agent that is added to 
the dissolver solution to enable solvent extraction of 
the uranium with TBP. Monel-400 is used for the 
dissolving vessel, and Carpenter-20, for the com- 
plexing tank. The solutions are handled thereafter in 
347, 348, 316, 304L, or 304 stainless-steel equipment. 

To increase fuel throughput, the process has re- 
cently been modified to permit semicontinuous disso- 
lution in 6.5M@ HF to which H;BO; is added to a con- 
centration of 4 g of boron per liter for nuclear safety. 
The 0.03M HNO; initially contained in the dissolver 
solution was found to be unnecessary in later opera- 
tions. Results of operating the semicontinuous pro- 
cess were reported previously. ’» Following that cam- 
paign, Zimmerman and Hoffman at ICPP conducted an 
evaluation of vessel corrosion which had occurred. 
The results of this evaluation have been reported 
recently.? 

The estimated rate of metal loss of the Monel-400 
dissolver, based on examination of a pipe specimen, 
was 7.5 mils/month of exposure over about a2-month 
period. This included losses incurred both in actual 
dissolution and in decontamination of the vessel. De- 
contamination was effected with hot (60 to 92°C) 7M — 
20M HF, 0.05M—0.1M chromic acid, and 6.0M H,SO, 
solutions over a time period of 244 hr. From labora- 
tory data the H,SO, decontaminant was most corro- 
sive (7 to 9 mils/month at 114°C). The total corrosion 
was not considered excessive since the dissolving 
vessel was constructed with a 188-mil corrosion al- 
lowance. Cyclic tests indicated that chromic acid cor- 
rosion of Monel-400 is increased by prior exposure 
to 20M HF. 

The 304L stainless-steel evaporator, used to con- 
centrate the product of solvent extraction, was also 
examined since there is some carry-over of fluo- 
ride to this unit. The evaporator also utilizes boron- 
containing 304 stainless-steel plates for nuclear 
safety. After 14,500 hr of exposure to product so- 
lutions from processing both zirconium-alloy and 
aluminum-alloy fuels (the latter solutions were 
fluoride-free), the estimated loss from the vessel 
wall was only 0.02 mil and that from the boron steel 
was 1.0 mil. A corrosion rate of 0.2 mil/month for 
the boron plates was calculated from the boron con- 
centration of the evaporator product. 


RECOVERY OF AMERICIUM AND CURIUM 
FROM SHIPPINGPORT REACTOR FUEL 


The Shippingport pressurized-water power reactor 
utilizes a seed-and-blanket type of core, in whicha 
ring of highly enriched uranium seed elements is sur- 
rounded by a blanket of natural uranium oxide. In the 
blanket, plutonium is formed and partially consumed. 
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Since the blanket is retained in this reactor for sev- 
eral years, appreciable amounts ofthe transplutonium 
isotopes americium and curium are also formed, A 
brief description was recently given by Rathvon, of 
the Atlantic Richfield Hanford Co., of the recovery at 
the Hanford Redox Plant of the americium and curium 
in connection with plutonium extraction from the first 
Shippingport blanket processing operation,” The blan- 
ket elements processed consisted of 13.5 tons of ura- 
nium (as UO,) in zirconium-alloy tubes and were 
estimated to contain 63 kg of plutonium, 1.2 kg of 
(44am + 43am), and 0.07 kg of “4Cm. The decladding 
and dissolution processes were previously described. 7 
After chemical decladding, the core material was 
dissolved in acid, and the uranium, plutonium, and 
neptunium were separated by solvent extraction with 
hexone, Americium and curium were then batch- 
extracted from the fission-product waste solution, 
with 30% TBP in a paraffinic diluent, stripped from 
the solvent with dilute HNO;, and concentrated by 
evaporation. Approximately 93% of the americium 
and curium were recovered by this means. No pro- 
cess details were given. 


NEPTUNIUM—PLUTONIUM RECOVERY 
FROM PUREX WASTE 


The accumulated, neutralized high-activity-level 
wastes from Purex processing operations at Hanford, 
which have been allowed to concentrate by self- 
boiling, contain significant amounts of neptunium and 
plutonium. These were accumulated prior to modifi- 
cation of the Purex process to provide for neptunium 
recovery. ” The transuranium elements are largely 
contained in a sludge in the waste-storage tanks, Be- 
cause plans have been made to sluice the sludge from 
the tanks with water to recover *’Sr from it, consid- 
eration has also been given to recovery of plutonium, 
as well as *°"Np, which is the starting material for the 
production of the isotopic power source 38pu (Ref. 
13). The conceptual recovery process that has been 
developed is based on coextraction of the transuranics 
from a nitric acid solution of the sludge with a diluent 
solution of di-2-ethylhexyl phosphoric acid (D2EHPA) 
and TBP. The utility of the D2EHPA reagent for 
heavy-element recovery has been discussed previ- 
ously.'4 Salt constituents of the sludge solution con- 
sist principally of the nitrates of iron, aluminum, and 
sodium at a total NO; concentration of about 3M (Ref. 
13), The transuranics are present at concentrations 
of 1 to 5x 10°°™M. 

The principal process operations consist of a 
sequence of batch contacts involving D2EHPA 
extraction—oxalic acid stripping (this permits use of 
available tankage for extraction), followed by evapo- 
rative concentration, This sequence of operations is 
repeated on the concentrated solution to provide fur- 
ther removal of uranium, iron, and other contami- 
nants. Because of its strong complexing action, oxa- 


late must be destroyed following the extraction 
sequences; this is accomplished by boiling the strip 
solution with HNO; and H,O). A simplified process 
diagram is shown in Fig. 1. 
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Fig. 1 Batch D2EHPA Np—Pu recovery process concept.'% 


For effective and simultaneous coextraction of 
neptunium and plutonium, it would be desirable to ad- 
just their valence states to the highly extractable 
quadrivalent state; however, all attempts to do this 
were unsuccessful. In the alternative procedure that 
was developed, Pu(IV), as found in the sludge solu- 
tion, is first extracted into the 0,1M D2EHPA—0.1M 
TBP solution in a paraffinic. diluent. After allowing 
the two phases to separate, but without removing the 
organic phase, the aqueous ‘phase is treated with 
0.03M ferrous sulfamate—0.01M hydrazine to reduce 
neptunium to Np(IV). Recontacting then extracts the 
neptunium without loss of plutonium. 


DETERIORATION OF EXTRACTION SOLVENTS 
UNDER IRRADIATION 

The limited irradiation stability of solvents used 
for the extraction and purification of fertile and fis- 
sionable materials has been indicated to be a factor 
of potential significance in defining the applicability 
of solvent extraction for fuel reprocessing.'° This is 
of particular concern in planning for the processing 
of fast reactor fuels because burnups of the order of 
100,000 Mwd/metric ton may be achieved in U—Pu 
ceramic fuels, and because economic considerations 
associated with the large fuel inventory in a fast re- 
actor may make it highly desirable to process fuels 
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after as short a cooling period as 1 month, Estimates 
of irradiation density in the aqueous feed to such 
processes are as high as 100 watts/liter. Blake'® has 
recently summarized the available information on 
methods of estimating the irradiation dose to sol- 
vents, as well as the literature on solvent deteriora- 
tion at expected levels, and has noted the specific 
effects of iodine and plutonium. It is observed that 
most of the reported data leads to erroneous and ex- 
aggerated predictions of solvent damage in process- 
ing; this is due in part to unavoidable confusion of 
chemical and radiation effects and in part to diffi- 
culties in estimating process doses realistically. A 
method of calculation is presented which is based on 
an equation developed for the radiation dose experi- 
enced in mixing. This is illustrated by the following 
equation applying to the operation of a mixer-—settler 
contactor: 


Em =Tm (Dg +GD,) P 


Em = organic-phase dose in mixing, watt-hr/ 
liter 
Tm = organic-phase residence time, hr 
D,,D, = beta and gamma power density in aque- 
ous feed, watts/liter 
G = fraction of gamma absorbed 
P = aqueous feed-volume fraction in mixer 


where 


An equation is also provided for exposure during 
phase disengagement. Literature data indicating ex- 
posures experienced in processing various thermal 
reactor fuels, where levels of up to 3 to 4 watt-hr/ 
liter per pass were noted, are evaluated on the basis 
of these equations. In comparison, estimated values 
for fast reactor fuel processing with short cooling 
are estimated to be only of the order of 0.1 watt-hr/ 
liter per pass, This value is thought to be quite tol- 
erable in consideration of the successful operation of 
extraction processes that provided adequate solvent- 
treatment equipment at accumulated solvent doses of 
tens of watt-hours per liter. 


In considering the effects of alpha radiation from 
high plutonium concentrations in fast reactor fuels, it 
was concluded that such radiation can contribute up 
to 10% of the total first-cycle dose in processing fast 
reactor fuels. In later extraction cycles the alpha 
effects are more significant. The prediction of the 
possible effects of the take-up of iodine by the sol- 
vent in processing short-cooled fuels is made diffi- 
cult by a lack of information of the behavior of iodine 
following dissolution. 


RESULTS OF EUROCHEMIC PLANT OPERATION 


The results of the first year’s experience in pro- 
cessing a variety of reactor fuels in the Eurochemic 
plant were described by R. Rometsch"" at the 1967 an- 
nual meeting of the American Institute of Chemical 
Engineers (AIChE) in New York in November 1967. 
This meeting was discussed briefly in the preceding 
issue of Reactor and Fuel-Processing Technology.® 
The Eurochemic plant structure and its process flow 
sheets were described in the Spring 1965 issue of 
Reactor Fuel Processing." In the first year of oper- 
ation, "” fuels from seven different reactors located in 
four countries were processed successfully (see Ta- 
ble 1). The core materials included both unalloyed 
uranium and U—Mo alloy varying in enrichment from 
natural uranium to 1.6% 7*°U. Fuels clad in both alu- 
minum and magnesium were handled satisfactorily, 
the former cladding being removed with NaOH and 
the latter with H,SO,. A total of 34.5 tons of uranium 
was processed in 208 on-stream days, with a pluto- 
nium recovery of 18.86 kg. 

Much of the material had undergone quite a long 
cooling period (3 to 9 years) so that high fission- 
product decontamination was not required. However, 
for the 3-month cooled material from the BR-1 re- 
actor, a first-cycle ruthenium decontamination factor 
of >10° and a beta decontamination factor of 5 x 10° 
were obtained for uranium. For Diorit, FR-2, and 
BR-1 fuels, minimum first-cycle zirconium decon- 
tamination factors were 70, 40, and 29, respectively. 


Table 1 FUELS PROCESSED DURING THE FIRST YEAR’S 
OPERATION OF THE EUROCHEMIC PLANT!’ 
(JUNE 1966—JULY 1967) 





Source of irradiated fuel (reactor and country) 











EL-1 EL-2A EL-2B Diorit FR-2 BR-1 Total 
France France France France Swiss W. Germany Belgium quantities 
Core U U U U U U 
Cladding Al Mg Mg Al Al Al 
Enrichment, % Nat. Nat. 0.72— 0.72— Nat. Nat Nat 
0.86 
Burnup, Mwd/ 25 200 300 1300 750 300 
metric ton 
Years cooled 9 2.5 3 1.6 1.6 0.25 
U, tons 2.0 3.8 5.6 6.1 6.7 1.9 34.5 
Pu, kg 0.03 0.65 1.06 6.9 4.1 0.61 18.86 
Process time, 17 31 37 84 208 
days 
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Plutonium was separated from uranium in the first 
cycle by a factor of 10*. Product specifications, fol- 
lowing two extraction cycles, a silica-gel treatment 
of the uranium, and final plutonium purification by 
oxalate precipitation, were met for both uranium and 
plutonium. 

It was felt that this initial operating period had 
adequately demonstrated the workability of the Euro- 
chemic multifuel plant concept, the feasibility of 
chemical decladding a variety of materials, and the 
operability of the natural-recirculation dissolver 
used. Somewhat more active waste was obtained than 
had been anticipated. This resulted principally from 
decladding operations, which amounted to 2 to 2.7 m® 
of U per ton for the Diorit, FR-2, and BR-1 process- 
ing runs. 


NEW AQUEOUS PROCESSING PLANT 


In a February 1968 news release,”’ Allied Chemical 
Corporation announced that it intended to file an ap- 
plication with the AEC for construction of a nuclear 
fuel-reprocessing plant designed to process spent 
fuel from U. S. power reactors. The plant will employ 
a modified aqueous process, will have a capacity of 5 
metric tons per day, and is expected to be operational 
in 1973. Primary consideration is being given to a 
location on the site of the AEC’s Savannah River plant 
near Aiken, S. C. No other details were released. 


Plant and Equipment 
Design and Operation 


PERFORMANCE OF CENTRIFUGAL MIXER-—SETTLER 


In spent-fuel processing, using solvent extraction, 
it has become more and more desirable to use a con- 
tactor that minimizes solvent exposure to fission- 
product radiation. This becomes even more desirable 
in processing breeder reactor fuel because econom- 
ics are improved by processing short-cooled fuel and 
thereby decreasing inventory charges. At a recent 
AIChE symposium on reprocessing irradiated fuels'® 
D. S. Webster presented a paper that related experi- 
ences using a centrifugal mixer—settler for sepa- 
rating plutonium and uranium from radioactive fission 
products.”! 

The unit tested is an 18-stage contactor made up of 
three 6-stage units, the design of which was previ- 
ously described in Reactor Fuel Processing.” 

The hydraulic capacity was tested prior to radio- 
active operations, and throughput rates up to the 
equivalent of 27 tons/day (70 gal/min) were achieved 
without excessive entrainment. Mass-transfer per- 
formance was also checked out in these cold runs, 
and less than 0.03% of the uranium remained in the 
waste stream. 

The performance of the unit, after a year’s opera- 
tion on radioactive fuels with fission-product gamma 


activity ranging from 50 to 200 curies/liter, was 
compared with that of conventional mixer—settler 
units: 

1. Mechanical reliability has been excellent (one 
defective unit was replaced remotely early in the 
run). 

2. Frequent shutdowns of the Purex process (over 
weekends) required only 15 to 30 min of solvent and 
scrub flow vs. 4 to 5 hr for pump-mix mixer— 
settlers. 

3. Steady-state operation on startup was reached 
in 20 min vs. 16 hr. 

4, Decontamination for the new unit was roughly 
equivalent; however, evaluation of decontamination 
performance is continuing. 

5. The short residence time lowered the gross 
gamma activity in the solvent to one-fifth of former 
values. (Higher-activity-level feeds will be tested.) 


This 1 year of operation, then, has demonstrated 
the advantages that were expected of the centrifugal 
mixer —settler; the demonstration of mechanical re- 
liability will require longer operating time. 


DELAYED-NEUTRON MONITOR FOR **U 
IN LEACHED HULLS 


One problem associated with the shear—leach pro- 
cessing” for spent fuels is making sure that no un- 
dissolved fuel is discarded with the hulls. When fuel 
material, such as UO:, is clad in stainless-steel or 
Zircaloy tubing, there is always the possibility that 
both ends of a short piece may be crushed closed by 
the shear, and hence not be leached by the acid. 

Strain et al.”4 have designed and evaluated a moni- 
tor to detect amounts of *°U as low as 10 mg in the 
presence of a gamma-radiation field of 10°r/hr. 

High levels of beta—gamma radiation can be ex- 
pected from both stainless steels (Co, tor, and 
58Co) and Zircaloy (Zr and *Nb). Table 2 gives an 


Table 2 DOSE RATES ESTIMATED FOR THREE SIZES OF FUEL 
LEACHING BASKETS THAT CONTAIN IRRADIATED 
STAINLESS-STEEL HULLS*4 


[Basis: Fuel Enrichment, 2% 2°5U; Burnup, 20,000 Mwd/ton; 
Irradiation Time, 6 x 10’ sec; Irradiation Flux, 3 x 10! neu- 
trons/(cm?)(sec); Basket Length, 150 cm; Depth of Hulls in 
Basket, 150 cm] 








Basket diameter, Decay time, Dose rate,* 
cm days r/hr 
13 14 12 nie 

30 9.05 x 103 
60 7.65 x 103 
120 5.91 x 103 
20 14 2.51 x 104 
30 1.91 x 104 
60 1.62 x 104 
120 1.26 x 104 
25 14 3.85 x 104 
30 2.95 x 104 
60 2.49 x 104 
120 1.93 x 104 





*At 8 cm from the outer surface of oasket. 
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indication of the various levels of radiation expected 
from various-sized batches of stainless-steel hulls. 

Four methods of detection were considered: 
(1) fission-product analysis by gamma spectrometry, 
(2) counting of neutrons produced by bombarding be- 
ryllium with the 2.18-Mev gamma radiation from 
Mce—iU4p, (3) differential gamma-ray absorption of 
fuel and cladding, and (4) delayed-neutron monitoring. 

The delayed-neutron monitoring method was se- 
lected for detailed evaluation and appeared to have 
the following advantages over the other three tech- 
niques: (1) variation in the type of cladding has no 
effect; (2) the analysis is independent of burnup or 
cooling time; (3) large quantities of hulls can be ana- 
lyzed rapidly; and (4) the analytical results can be 
handled by computers or nontechnical personnel. 

The theory of this method is that, when a fission- 
able material such as 7°U is bombarded by thermal 
neutrons, fission occurs, accompanied by the forma- 
tion of short-lived products that decay by emission of 
neutrons. Most of these radionuclides are isotopes of 
bromine and iodine, for example, Br, Br, 1371, etc., 
whose half-lives range from 0.23 to 54 sec. An aver- 
age of 0.02 delayed neutron is produced per ay fig- 
sion. This method has been demonstrated to detect 
7x10 °¢g of *°U per liter in a 0.5M nitric acid solu- 
tion after exposure in a reactor operating at 2.5 Mw 
and a neutron flux of 2.1 x 10’ neutrons/(cm’)(sec) 
(Ref. 25). 

The neutron generator used in the evaluation study 
was a commercially available Cockcroft—Walton ma- 
chine. The 14.7-Mev neutrons produced were mod- 
erated by paraffin or polyethylene to provide asource 
of thermal neutrons; however, the authors proposed 
water for a permanent installation. 

The delayed neutrons could be detected by measur- 
ing the products of neutron-capture reactions by 
means of such instruments as: (1) *°U fission cham- 
ber, (2) *He-filled proportional counter, (3) '°BF,- 
filled proportional counter, and (4) !°B-lined counters. 

A fission chamber was not considered practical for 
incorporation in the monitor because there is suffi- 
cient *°U in the fission chamber to produce an ap- 
preciable background of delayed neutrons when the 
fission chamber is subjected to neutron bombardment. 

Proportional counters containing *He or '°B make 
use of the reactions *He(n,p)°H or ‘°B(n,a)'Li to de- 
tect neutrons. The positively charged particles ionize 
the filler gas within the counter, thereby producing 
an electronic pulse that can be amplified and 
recorded. 

Various experiments and calculations by the au- 
thors led to the choice of a recently developed boron- 
lined chamber, which, combined with an electronic 
system having a nanosecond response, was used suc- 
cessfully to count neutrons in a gamma flux of 10° 
r/hr. 

Two prototype monitors were built and tested. 
One was designed with separate compartments for 
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irradiation and detection; the other had a single 
compartment. 

The double-compartment monitor was designed to 
assay small batches of leached hulls that had been 
removed from the leachant and the rinse tanks by 
means of a cup conveyor. In a proposed in-cell in- 
stallation of such a system (see Fig. 2), the hulls 
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Fig. 2 Proposed in-cell installation of double-compartment 
leached-hull monitor.*4 


from a single cup are deposited in the monitor for the 
period of time required for assaying for fissionable 
residues and are then dropped upon one of two con- 
tinuous conveyors. These conveyors either return the 
hulls to the leaching tank for further treatment or 
transfer them to a collector for subsequent disposal. 

The sensitivity of the double-compartment proto- 
type monitor for detecting *3517 is demonstrated in 
Table 3. 

In a variation of the shear—leach process, the hulls 
are carried through the leach and rinse steps ina 
porous, stainless-steel basket 20 cm in diameter and 
up to 150 cm long. The capacity of each basket is ap- 
proximately 50 kg (several thousand hulls). A single- 
compartment monitor was designed to assay the 
leached hulls without removing them from the basket. 
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Table 3 SENSITIVITY OF DOUBLE-COMPARTMENT MONITOR*4 





Repetitive analysis 





Sample (counts/20 sec)* Average* 
Background 

(20 hulls) 5, 9, 6, 4, 7, 6 6 
1 mg of 255 17, 15, 15 16 
5 mg of 2% 32, 47, 52, 54, 48 47 
6 mg of 235 38, 38, 60, 58 49 
100 mg of 235y 520, 583, 686, 506, 650 


910, 662, 687 





*Background not subtracted. 


A proposed in-cell installation of this monitoring unit 
is shown in Fig. 3; the prototype single-compartment 
monitor is shown in Fig. 4. 

The authors include detailed discussions of each 
component of the system and conclude that 1 mg of 
23517 could be detected with the two-compartment de- 
vice and 10 mg of 7“°U, with the single-compartment 
system. The double-compartment system has the ad- 
vantage of speeding up the process, since analysis 
and disposal of leached hulls could start at any 
point in the process. In comparison, the single- 
compartment system, which operates on larger units 
or baskets of hulls, may require recirculation of an 
entire basket if undissolved 7°°U is detected. 


The single-compartment model, however, would 
require less cell space, would be mechanically sim- 
pler, and lower in cost. The authors estimate the cost 
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Fig. 3 Proposed installation of the single-compartment 
leached-hull monitor.*4 


Fig.4 Prototype single-com- 
partment leached-hull moni- 
tor.4 The target of the neu- 
tron generator is inserted in 
the polyethylene moderator. 
The radiation shield for the 
Co gamma source can be 
seen beneath the moderator 
support. 
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to be less than $50,000 for the entire system. They 
also point out that the delayed-neutron monitor sys- 
tem might have other applications, such as analysis 
of uranium-containing process streams. 
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Volatility Processes 


By G. J. Vogel and J. W. Simmons 


Fluid-Bed Fluoride-Volatility Processes 
for Ceramic Oxide Fuels 


Fluid-bed fluoride-volatility processes are under de- 
velopment for the recovery of uranium and plutonium 
from irradiated UO, fuel. The processes reported 
here use approximately the same steps—chemical 
removal of cladding, oxidation of the ceramic oxide 
to a fine powder, fluorination to form volatile actinide 
hexafluorides, and purification of hexafluorides prior 
to their conversion to oxide for recycle tothe reactor. 
At different sites, different fluorination agents (Fo, 
BrF;, C1F3;, and C1F) are being evaluated. Develop- 
ment work at the Oak Ridge Gaseous Diffusion Plant 
(ORGDP), Oak Ridge National Laboratory (ORNL), 
Centre d’Etude de 1’ Energie Nucléaire (Mol, Belgium), 
Argonne National Laboratory (ANL), Battelle Memo- 
rial Institute, and Fontenay-aux-Roses is reviewed. 

Two general review articles were published re- 
cently. Fluidized-bed fluoride-volatility reprocessing 
methods being developed at ANL are described ina 
proceedings, ! and another article describes fluoride- 
volatility methods under development at AEC labora- 
tories for the reprocessing of plutonium-containing 
fuels.” 


DECLADDING AND OXIDATION OF OXIDE FUELS 


An experimental program is under way at ORGDP 
to determine the processing conditions and equipment 
requirements for optimum batch decladding and oxi- 
dation of UO, pellets.*** In one run, a 76-kg simulated 
fuel element (a 6-in.-square array of sixty 4-ft-long 
tubes loaded with 53 kg of annular UO, pellets) was 
processed through the decladding and oxidation steps 
in a semiworks-scale fluid bed (45.5 kg of 48- to 100- 
mesh alumina). In this run, decladding by treatment 
of the Zircaloy with 20 to 50 vol.% HCl followed by 
oxidation with 10.5 to 21 vol.% oxygen resulted in an 
overall HCl utilization of 29% during decladding and 
in elutriation from the reactor of 97% of the uranium 
charged. The elutriated material contained essentially 


no alumina, as is necessary if transfer by elutria- 
tion is to become practical. 

Four oxidation runs were performed in the same 
equipment with the charge consisting of ~55 kg of 
alumina or magnesium fluoride powder and 65 kg of 
UO, pellets (either annular or solid pellets).‘ It was 
concluded from the results of these runs that the re- 
actor can be operated successfully through the de- 
cladding and oxidation steps with either alumina or 
magnesium fluoride as the bed material, that a super- 
ficial gas velocity of ~2.5 ft/sec is needed for oxida- 
tion of a 2°/,-ft-deep bed of solid pellets, and that a 
gas velocity of ~1.5 ft/sec is needed for a2*/,-ft-deep 
bed of annular pellets. 

Reaction of stainless-steel cladding with HF—O, in 
a fluidized bed produces oxides of stainless-steel con- 
stituents. Reaction of Zircaloy cladding with HF-—O, 
may produce ZrO, shards.°> The decladding gas mix- 
ture may react with the fuel oxide also, producing 
particulate fluorides and oxides in a proportion de- 
pendent upon the composition of the gas. 

Two decladding runs were made at ORNL’ in which 
Zircaloy-clad fuel elements were declad with 20 to 40 
vol.% HF-—O,. A fuel assembly consisting of four 12- 
in.-long, 0.4-in.-OD Zircaloy tubes containing UO, 
pellets was used in each run. Runs were performed 
at 550°C in the presence and absence of a fluidized 
bed of alumina. Although it may be mandatory to re- 
move heat in this step by means of a fluidized bed, 
results indicated that exposure of the oxide to a 
fluidized bed should be minimal to minimize the 
amount of fines produced and to minimize uranium 
and plutonium losses. 

Bartlett and associates® patented a process for 
fluid-bed removal of zirconium cladding from fuel 
elements with a gas composed of 20 to 90 vol.% 
oxygen and 10 to 80 vol.% hydrogen fluoride at 500 
to 800°C. At Argonne it was found that a gas mix- 
ture containing at least 2 parts by volume of hydro- 
gen fluoride to 1 part by volume of oxygen does not 
react with the UO, pellets.’ A high-HF gas mixture 
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would thus react with and remove stainless-steel 
or Zircaloy cladding without affecting the UO). 

At Mol, one process cycle under investigation in- 
volves decladding of stainless-steel-clad UO, by re- 
action of the cladding with HF—oxygen mixtures at 
500 to 600°C. A series of nine 1-hr experiments was 
carried out® at 590°C to determine the effect of con- 
centration of the gaseous reagents on (1) the rate of 
the decladding reaction and (2) the degree of pulveri- 
zation of the uranium dioxide fuel pellets. There was 
no fluidized bed in the reaction vessel. Results 
showed that hydrogen fluoride concentrations of 40 to 
65% gave satisfactory decladding rates and that with 
65% HF in the gas mixture the rate of pulverization 
of the oxide pellets was at a maximum. 

Fuel assemblies, each consisting of nineteen 20- 
cm-long stainless-steel-clad fuel elements packed 
with natural UO,, were chemically declad at Mol.? In 
two tests the gas mixtures contained 15 to 40% HF, 
28 to 60% O,, and 0 to 53% N,. Heat transfer between 
the center of an assembly and the fluid bed was poor, 
and the temperature at fuel-element surfaces ranged 
from 500 to 800°C, indicating that the 2-mm distance 
between fuel elements in these assemblies was inade- 
quate. Study of chemical decladding was continued, 
using assemblies containing seven fuel elements 
5 mm apart.'° Each fuel element consisted of a 304L 
stainless-steel tube containing 20 depleted UO, pellets. 
Initial tests showed that prefluorination ofthe alumina 
bed with HF prevented agglomeration of alumina par- 
ticles during the decladding step. Four tests (three 
with prefluorinated alumina and one test with none) 
were performed at bed temperatures up to 650°C 
using 40% HF-60% O, as the reagent. Although heat 
transfer was very poor in the test in which alumina 
was absent, decladding was complete in this and the 
three other tests, and the uranium dioxide was 
pulverized and transformed into UO,F, and UF4. 
The authors concluded that it is probably not neces- 
sary to perform the decladding reaction in a bed of 
fluidized alumina. 

In a series of fluid-bed tests at Mol!’ in which 
4000-g charges of compacted uranium dioxide—2% 
fissium in 2000 to 4500 g of Al,O; were oxidized at 
500°C with 20 to 30% oxygen in a 10-cm-diameter 
reactor, only a small variation in temperature was 
observed along the bed height, and the uranium con- 
tent in the alumina increased uniformly as a function 
of time. At the end of the tests, no agglomerates were 
found and only a low percentage of unreacted UO, was 
present. The decrease in the proportion of alumina to 
one-third of the total charge did not pose new prob- 
lems. 


FLUORINATION WITH FLUORINE 


A series of experiments was performed at ANL 
in which 0.3- to 1.2-g charges of simulated spent fuel 
pellets containing UO,, PuO,, and nonradioactive 
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fission-product oxides were fluorinated with fluorine 
in a 2-in.-diameter fluidized-bed reactor. Each ex- 
periment was evaluated primarily on the amount of 
plutonium remaining in the alumina bed following 
recycle fluorination with 90 vol.% fluorine at 350 to 
550°C or 450 to 550°C. It was found that a 12-hr 
recycle-fluorination sequence gave as good plutonium 
removal as a 20-hr sequence. In other experiments 
in which various pellet-to-alumina weight ratios were 
used, it was found that the fraction of plutonium re- 
tained by the alumina decreased with an increase in 
the pellet-to-alumina charge ratio in the range 0.3 
to 1.2. In another experiment a single bed of alumina 
was used to process three batches of pellets; the plu- 
tonium retention on the final bed was 1.2% of that 
charged, which was the minimum value attained for 
this series of experiments. Another beneficial effect 
resulting from the reuse of alumina for several 
batches of fuel pellets is the reduction of the amount 
of solid waste. To minimize the amount of plutonium 
not reacted because of holdup on the upper cooler 
sections of the reactor, workers at ANL'® performed 
an experiment with a bed consisting entirely of —100 
mesh alumina instead of normally used 48 to 100 
mesh. After precoating of the off-gas filters with 
alumina, it was expected that the bulk ofthe plutonium 
fines would be returned to the reaction zone by blow- 
back of the filters. At the end of this run, excessive 
amounts of plutonium were found in the bed and on 
the upper sections of the reactor, and this fact indi- 
cates that the use of —100-mesh alumina as the bed 
does not decrease plutonium holdup in the upper sec- 
tions of the reactor. 


Three tests were performed at ANL'® to evaluate 
the removal of low concentrations of plutonium from 
alumina-bed samples by leaching with nitric acid, as 
an alternative to extensive fluorination of the bed. 
Samples containing 0.164 wt. % plutonium were leached 
for 2, 3, or 5 hr with 6N HNO, at 100°C. The per- 
centages of plutonium removed by leaching in the 
three tests were 95.1, 99.8, and 99.8%, 


Three campaigns of PuF, fluorination experiments 
were performed at ANL to demonstrate that PuF, can 
be produced and transported readily and collected 
quantitatively in a 3-in.-diameter fluid-bed reactor 
facility. The results of two of the three campaigns 
were reported earlier.'* The third campaign, like the 
other two, consisted of three separate fluorination 
runs (fluorination of about 260 g of PuF, with fluorine 
in each run), followed by a cleanup run in which plu- 
tonium deposited in the equipment was recovered by 
fluorination with fluorine.'> Three-hour fluorinations 
were used in this campaign instead of 5-hr fluorina- 
tions used earlier, and the mean production rate was 
4.1 lb of PuF,/(hr)(sq ft) instead of 2.4 lb of PuF¢/ 
(hr)(sq ft) obtained in the earlier campaigns. The 
fluorination rates, although substantially lower than 
those for uranium, are adequate for practical appli- 
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cations. A material balance for the third campaign 
shows that 98.9% of the plutonium charged was ac- 
counted for.'® This value is within the range (96 to 
104%) expected for a good material balance. 

Other engineering-scale studies were performed 
in the same ANL equipment to develop information 
on the basic behavior of PuF, in the fluorination pilot 
plant. !” Emphasis of the work was directed toward 
the step in which PuF, is fluorinated with fluorine 
for plutonium recovery. The material balance for 
six engineering-scale experiments and an additional 
engineering-scale plutonium-transfer experiment was 
98.7%, which is considered to be satisfactory.'* Intwo 
experiments, fluorination efficiency (defined as the 
amount of PuFs produced in a given period divided by 
the theoretical amount produced at equilibrium) was 
determined. Efficiency was near 50% at the start of 
each experiment and decreased rapidly as the pluto- 
nium content of the bed diminished. In the plutonium- 
transfer experiment, successful transfers of UF, and 
PuF, were carried out with 160.2 g of UF, and 89.4 g 
of PuF.. Recoveries were 99% and 98.4%. The results 
indicated that, with proper conditioning of equipment 
(prefluorination with fluorine or perhaps even with 
PuF,), the transport of PuF, in engineering-scale 
equipment is feasible. 

Instrumental methods for plutonium analysis of 
fluoride-volatility process streams are under devel- 
opment 13417 ot ANL as possible alternatives to chemi- 
cal analytical methods. A water-moderated neutron 
detector using *He proportional counters was con- 
structed and was tested with plutonium fluorides. Its 
specific neutron counting rate was determined to be 
20,800 + 400 counts/min per gram of PuF4, which is 
sufficiently sensitive and accurate for some process 
uses and indicates any changes in plutonium concen- 
tration at low levels with suitable uniformity to per- 
mit useful comparisons. A second type of neutron 
detector evaluated was a portable instrument with 
neutron shielding to minimize interference from any 
plutonium fluorides unavoidably near the process 
equipment being surveyed. Practical levels of direc- 
tional discrimination could not be achieved without 
sacrificing the portability feature of the latter detec- 
tor. 


FLUORINATION WITH ClF3 AND CLF 


Workers at Centre d’Etude de l’Energie Nucléaire 
(CEN) at Mol, Belgium, are developing a fluid-bed 
process for the recovery of uranium and plutonium 
from (U,Pu)O, fuel. A proposed flow sheet® provides 
for the reaction of compacts of (U,Pu)O, with a 
chlorine fluoride to selectively fluorinate uranium to 
volatile uranium hexafluoride. Plutonium is recov- 
ered in the following step as PuFs by reaction of the 
residue with fluorine. 

To be an effective agent for separating uranium 
from plutonium, a chlorine fluoride should react with 


uranium to form UFs but should not form PuF,. A 
thermodynamic study was carried out to allow estima- 
tion of the PuFs, concentrations that would be in 
equilibrium with various fluorinating agents. On the 
basis of Fig. 1 (Ref. 9), which shows that insignifi- 
cant amounts of PuF, are formed when CIF is in 
equilibrium with plutonium at temperatures of in- 
terest, CIF is being considered as a selective fluo- 
rinating agent. This reagent may be used as the 
selective fluorinating agent for uranium if the fuel 
is pulverized prior to the fluorination step. If the fuel 
is not pulverized, C1F; may be used as the selective 
fluorinating agent. '® 
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Fig. 1 Partial pressure of PuFg in equilibrium with ClF3;— 
CLF mixtures or F2 as a function of the temperature (atmo- 
spheric pressure).? 


In the investigations reported here, decladding of 
the oxide fuel at Mol was accomplished by reaction 
with HF—O,, which also pulverized the fuel to form 
mixed oxide and fluoride. After one stainless-steel- 
clad fuel assembly had been declad, the oxidized ura- 
nium and decladding products were reacted with 10 to 
25% CIF at 250 to 450°C to volatilize the uranium.”° 
Most gases from the reaction (but not UF,) were 
sorbed on a bed of sodium fluoride at 350°C. Analysis 
of samples from the bed indicated that volatile com- 
pounds of nickel, iron, and chromium had formed 
during the reaction with C1F. 

Laboratory-scale tests were performed at Mol? to 
determine the extent of reaction of plutonium tetra- 
fluoride with CIF, and ClF at selected temperatures. 
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Little plutonium was volatilized when PuF, was con- 
tacted with 30% CIF, in argon below 300°C. In four 
other tests, ° C1F was contacted with 50-mg samples 
of PuF, at 550°C. In two tests in which the gas mix- 
ture was 30% CIF in argon and in two other tests in 
which the mixture was 20% Cl1F —30% oxygen inargon, 
no volatilization of plutonium was observed. 

Other laboratory tests!’ were performed to deter- 
mine the effect of high plutonium concentration in 
(U, Pu)O, pellets on fluorination by C1F3. One investi- 
gation was carried out with an alumina bed using high- 
density fuel pellets. Upon fluorination at 300°C of 
pellets containing 0 to 20% PuO,, fluorination and 
volatilization of uranium from the alumina bed were 
complete. However, with plutonium contents of 20% or 
greater, the rate of uranium volatilization was sharply 
reduced. No volatilization of plutonium was observed 
during volatilization of the uranium. It was concluded 
that C1F, at 300°C can be used for direct fluorination 
of high-density mixed-oxide pellets. At the end of the 
runs, the residue was observed to be pulverized. 


In another investigation, UO,—15% PuO, micro- 
spheres were fluorinated with 30% CIF, in argon. As 
determined with a thermobalance, the rate of volatil- 
ization of the uranium was much lower than the 
volatilization rate for pure uranium dioxide in the 
same physical state. 

The influence of various ratios of CIF to CIF; (in 
the fluorination gas) on the rate at which PuF, is 
fluorinated at 350°C was also studied at Mol.'° At 
high ClF concentrations (70 to 90%) in the fluorina- 
tion gas, the rate of plutonium volatilization de- 
creased to one-tenth the rate observed when no C1F 
was present. 

Experiments were performed at Mol® to determine 
the effects of (1) temperature, (2) flow rate of the 
gases, and (3) concentration of the fluorinating agent 
(in nitrogen) on the quantity of uranium reacting with 
the fluorinating agent. Three series of experiments 
were performed in a 3.75-cm-diameter column; 
(1) the fluidized-bed fluorination of a UF,—UO,F,-— 
alumina mixture with C1F, (2) the fluidized-bed fluo- 
rination of a UF,—UO,F,—alumina mixture with C1F3, 
and (3) the fluorination of UO, with C1F;. In each run 
of a series, the same total amount of fluorinating 
agent was fed. With increased temperature, a sig- 
nificantly greater reaction occurred in each series of 
experiments. Increasing the gas flow rate decreased 
the total amount of reaction inthe runs in series 1 and 
3 above but increased the total amount of reaction in 
series 2. In series 2 and 3, an increased concentra- 
tion of the fluorinating agent gave a decrease in the 
total amount of reaction, whereas in series 1the con- 
centration of the fluorinating agent had no significant 
effect. Other tests at Mol® in which a UF,—UO,F,— 
alumina mixture was fluorinated with 15, 25, 35, and 
50% C1F showed that the reaction rate is a linear 
function of C1F concentration. 
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In a series of 16 tests,’ the quantity of UO, re- 
acted with 35% C1F; in a fluid bed was determined for 
two values of each of four variables: fluorination 
temperature (100 or 300°C), the presence or absence 
of fissium, the amount of alumina (0 or 20% of the UO, 
weight), and the ratio of height to diameter of the bed 
(3 or 6). The reaction time and amount of gaseous re- 
agent were the same in all tests. The maximum 
amount of UO, was reacted in the test carried out at 
300°C with a charge of 1000 g of UO, containing no 
fissium, in the presence of 200 g of alumina, and in a 
bed with a height-to-diameter ratio of 6. 

In other tests at Mol,’® UO, pellets were fluorinated 
with 26 or 31% CIF; in a 10-cm-diameter fluid-bed 
reactor, and the height of the reaction zone was mea- 
sured. With a gas flow of 6 cm/sec, the reaction zone 
had a height of about 20 cm. 


Behavior of Neptunium 


In some conceptual versions of the fluid-bed fluo- 
ride-volatility process, a gas stream containing UFg, 
PuF,, NpF,, and other volatile fluorides would be 
passed through a fixed bed of NaF (operated at 100 to 
400°C in various steps of the process).!° Neptunium 
hexafluoride would be removed from the gas stream 
by the reaction 


3NaF(s) + NpF,(g) = 3NaF - NpF;(s) + F(g) (1) 


The neptunium can be subsequently recovered from 
the 3NaF - NpF; complex as NpF, by fluorination with 
fluorine. 

Data on the extent to which NpF¢ would be taken up 
or released by a bed of NaF over a range of fluorine 
pressures is needed for a proposed method of sepa- 
rating NpF, from UF¢. The equilibrium constant 


Ky = (Prope) /(P¢,)” 


for the fluorination of 3NaF - NpF; by fluorine has 
been determined at a pressure of 700 torrs and tem- 
peratures’® of 250 to 400°C. It can be expressed as a 
function of temperature by the equation 


—3.147 x 10° 


= + 2.784 
T(°K) 


log K, (atm’?) = 


Several measurements, made at 350°C with fluorine 
pressures varying from 150 to 1500 torrs, showed 
that log Pypr, was a linear function of log Pr, with a 
proportionality coefficient of 0.49, which is compa- 
rable to the value of % expected for reaction 1. The 
value of the equilibrium constant, Kp, at 350°C was 
5.52 x 107% atm. 

From the values of the equilibrium constant, the 
following values for thermodynamic functions for the 
fluorination of 3NaF - NpF; have been calculated over 





NONAQUEOUS PROCESSING 213 


the range 250 to 400°C: mean enthalpy change, 14.4 
kcal/mole NpF,; mean entropy change, 12.7 cal/ 
(mole NpF,)(°K). The standard free-energy change at 
298°K was 10.6 kcal/mole NpFy. 


Oxidation and Hydrofluorination 
of Carbide Fuels 


At Mol the pulverization of uranium monocarbide 
by reaction with HF, O,, water vapor, CO,, and mix- 
tures of these reagents in a fluid bed, followed by 
fluorination of the pulverized material, is being in- 
vestigated,!! Compact UC in the form of high-density 
pellets was completely converted to U;O, with 44 vol. % 
oxygen at 500°C in one test and was completely con- 
verted to UO,F,—UF,—uranium oxide mixtures in 
two other tests in which HF—O, was used as the fluo- 
rinating agent at 450 or 500°C. 


Purification of Uranium Hexafluoride 
and Plutonium Hexafluoride 


In the fluid-bed fluoride-volatility processing of 
reactor fuel materials, volatile UF, and PuF, are 
produced by reaction of the fuel with gaseous fluo- 
rinating agents. Uranium may be separated from 
plutonium by selective fluorination with CIF, C1Fs, 
or BrF;; or, alternatively, uranium and plutonium 
may be simultaneously volatilized by reaction with 
fluorine and separated from each other in a subse- 
quent step. In these reactions, volatile-fluoride com- 
pounds of certain fission products are also produced 
and appear in the UF, and PuFs process streams. 
Some of the recent work on developing methods 
for removing the fission-product contaminants from 
process streams is described below. 


PURIFICATION OF URANIUM HEXAFLUORIDE 


An earlier issue of this review’’ discussed the use 
of NaF and MgF, sorbents to remove volatile impuri- 
ties from uranium hexafluoride prior to its return to 
the gaseous diffusion complex. A paper has been pub- 
lished describing the investigation.”! 


PURIFICATION OF PLUTONIUM HEXAFLUORIDE 


Cleveland” briefly described the methods, ad- 
vantages, and disadvantages of several techniques for 
separating plutonium from irradiated uranium in 
fluoride-volatility processes. 

The French are investigating the selective chemical 
reduction of plutonium hexafluoride with carbon diox- 
ide.”*> Tests showed that uranium hexafluoride and 
niobium pentafluoride apparently are stable in the 
presence of carbon dioxide and that ruthenium penta- 
fluoride reacts only slowly. Mixtures of UF, (30 to 60 
torrs), PuF, (30 to 70 torrs), and CO, (600 to 950 
torrs) were reacted at 200, 300, and 350°C and re- 


sulted in complete reduction of PuF, in15 min, 1 min, 
and 30 sec, respectively. 

Workers at ORNL are studying a purification 
procedure in which PuF, reacts with LiF to form a 
complex from which PuF, can later be conveniently re- 
covered by treatment with fluorine.”?""4" The tech- 
nique of purifying PuF, process-gas streams based 
on the sorption of PuF, on granulated LiF is also 
being studied at ANL.‘*'%?6 Two experiments were 
performed”® in which mixtures of ruthenium fluorides 
and plutonium fluoride were passed through a static 
LiF bed to determine whether stabilities of the com- 
plexes differ sufficiently so that (1) a large fraction 
of the ruthenium would be fixed as the solid complex 
and (2) a large fraction of the plutonium would remain 
in the vapor phase. In a 10-hr experiment with the 
bed at 300 to 400°C, the ruthenium decontamination 
factor was 2.5 x 107, and 1% of the plutonium was re- 
tained in the bed. In a 7-hr experiment at 300°C, the 
corresponding values were 9 x 10° and 26%. 

Two additional experiments were performed!® at 
ANL in which a gaseous mixture of PuF,, ruthenium 
fluoride, and fluorine was passed through LiF beds at 
400 and 500°C to determine the effect of bed temper- 
ature on the decontamination factor (DF) and on plu- 
tonium retention in the bed. Results of these and 
earlier experiments show that the best temperature 
for PuF,g—ruthenium fluoride separation is between 
200 and 500°C. 

In further work at ANL,’’ a PuF,—fluorine mixture 
was reacted with LiF at 300°C until the complex that 
formed contained 29.4 mole % PuFy. X-ray powder- 
diffraction analysis of the solid showed the presence 
of the LiF - PuF, complex. The sample was then 
heated to 400°C, fluorinated with fluorine, and weighed 
at intervals. A change in the rate of weight loss was 
observed during fluorination, suggesting that a second 
plutonium-containing complex was present in the final 
stage of fluorination. The presence of 4LiF- PuF, 
was confirmed by X-ray analysis of samples taken 
after the change in reaction rate. 

Several runs were performed at ANL in which the 
loading capacity of LiF for ruthenium fluoride was 
determined'® at 300, 350, and 400°C. The loading 
capacities expressed as milligrams of ruthenium per 
gram of LiF were 9 to 13 at 300°C, 17.3 at 350°C, 
and 9 at 400°C. 

A comparison was made of the behavior of ruthe- 
nium fluoride species generated by the fluorination 
with a 3:1 fluorine—oxygen mixture of (1) ruthenium 
metal and (2) ruthenium dixode.'* The volatile product 
was trapped in a cold trap at —78°C and then tran- 
spired by increasing the cold-trap temperature in 
50°C steps from 0°C to 150°C. The products of both 
starting materials behaved similarly —for example, 
both volatilized at about 100 to 150°C. 

To determine whether sublimation of PuF, from a 
cold trap at 0°C or 50°C, leaving ruthenium behind, 
might be a feasible technique of separating plutonium 
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from ruthenium, ruthenium fluoride and UF (astand- 
in for PuF,) were transpired from a cold trap at 0, 
50, 100, and 150°C (Ref. 19). Decontamination factors 
greater than 5.9 x 10° were calculated for four ex- 
periments. 


REACTION OF MIXED HEXAFLUORIDES 
WITH SILICON TETRACHLORIDE 


Speeckaert”” at Mol has proposed separating ura- 


nium from plutonium by converting them to chlorides. 
Uranium hexafluoride and plutonium hexafluoride are 
converted to the chlorides by reaction with SiCl,: 


UF, + 1.5SiCl, = UCl, <.c 1.5SiF, 
A G93 = —96 kcal/mole 


PuF i 1.5SiCl, = PuCl, a3 1.5SiF, + 1.5Cl, 
AG298 = —120 kcal/mole 


The separation of uranium and plutonium would be 
based on the difference of vapor pressures of the 
chloride species formed. The highest chlorides formed 
at UClg and PuCl,—PuCl;. The ratio of the vapor 
pressures of UClg and PuCl, at 300°C is about 1’, 
and the corresponding ratio for UCl, and PuCl, is 
many orders of magnitude larger. 


SEPARATION OF URANIUM FROM PLUTONIUM 
IN MOLTEN SALT 


The problem of separating uranium from plutonium 
and converting them to the dioxides for recycle to the 
reactor was examined at Mol.'' A process based on 
oxidizing uranium and plutonium in a medium of 
molten alkali chlorides was chosen. Plutonium and 
uranium hexafluorides may be dissolved in a molten- 
salt bath such as NaCl, KCl, or LiCl. Upon oxidation 
with HCl1—O,, HCl—air, Cl,—O,, or an alkali nitrate, 
uranium is converted to soluble uranyl fluoride and 
PuO, precipitates. Removal of the PuO, by filtration 
is followed by electrolysis to deposit UO, on a 
cathode. In tests carried out with a LiCl—KCl eutectic 
at 500°C, the quantity of plutonium dissolved in the 
salt after PuO, precipitation was generally found to 
be satisfactorily small. In tests in which uranium was 
charged in the proportion of 9 parts to 1 part of plu- 
tonium, no uranium was detected in the PuO, product. 


Filter Studies 


The experience at ANL with automatic blowback 
systems for primary filters in fluid-bed reactors has 
been reviewed in a topical report.” The effects of 
variables such as pulse duration, interval between 
pulses, pressure of blowback gas, and design of the 
nozzle-venturi assembly were investigated. It was ob- 
served that successful operation has been achieved 
with a wide variety of laboratory and pilot-plant units. 
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Battelle-—Columbus”’ studied the plugging of nickel 
off-gas filters that may occur in the decladding step 
in which Zircaloy is reacted with a HCl—nitrogen 
stream in a fluid bed to form volatile ZrCl,. The ob- 
jectives of the bench-scale, simulated decladding runs 
were to determine the causes of plugging and to find 
ways of alleviating this problem. It was determined 
that a major source of plugging was NiCl, formed in 
the pores of the filters by reaction of HCl with filter 
material. Although it was found that filter plugging 
from this reaction is reduced by lowering the operat- 
ing temperature of the filter to ~350°C, the minimum 
temperature is fixed by the need to volatilize all of 
the ZrCl, (sublimation temperature, 350°C), Another 
major source of plugging was deposition on the filter 
surface of fines formed when Zircaloy was hydro- 
chlorinated. These fines were largely ZrO, and were 
produced by reaction with oxygen present in the Zir- 
caloy or in trace amounts in the system. Deposition 
of alumina fines from the bed was a minor factor. 
Suggested as an alternative to nickel filters isa filter 
consisting of a packed bed of alumina, 


At ORGDP, several filter materials were evaluated 
for their ability to remain permeable during repeated 
cycles of hydrochlorination with HCl, oxidation with 
oxygen, and fluorination with fluorine.*4 Tempera- 
tures differed for different runs, and charge compo- 
sitions differed: either zirconium metal or Zircaloy, 
UO, present or absent, tabular alumina or prefluo- 
rinated magnesium fluoride. In one test, twoin-series 
sintered-nickel-fiber filters at 315°C were exposed 
to off-gases from the reaction of hydrogen chloride 
with zirconium metal in an alumina bed to determine 
whether plugging increased or decreased in the ab- 
sence of nonzirconium components of Zircaloy. Plug- 
ging developed more slowly in this 10-cycle test than 
in tests using a Zircaloy charge but reached essen- 
tially the same level of permeability by the end of the 
test series. In the next test a highly permeable filter 
fabricated of woven nickel wire anda sintered-nickel- 
fiber filter were installed in series, and filter tem- 
peratures were maintained at 315, 100, and 100°C dur- 
ing hydrochlorination, oxidation, and fluorination of the 
alumina—UO,-—Zircaloy charge. Excessive plugging 
within three cycles lent support to the possibility that 
plugging is of a chemical nature. Substitution of pre- 
fluorinated magnesium fluoride for alumina as the 
bed material, in atest similar to the preceding test, 
resulted in excessive plugging in the fifth cycle, 
indicating that this substitution would not solve 
the plugging problem. In another test an alumina— 
Zircaloy charge was exposed to decladding —oxidation — 
fluorination cycles, and the off-gases were passed 
through two filters fabricated of Inconel powder. Ex- 
cessive plugging of these filters occurred in the sixth 
cycle, indicating that they are not suitable for the 
three-step operation. 
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Heat Generation by Radioactive Materials 


Dunthorn at ORGDP*® studied the effect of heatup 
of an irradiated fuel element (a square array of fuel 
rods) and of the bed itself in fluid-bed processing. 
Information was needed to establish the following: 
design philosophy for equipment, the actions to be 
taken in the event of equipment malfunction, limits 
on the radiation levels of the fuel elements to be 
processed, and the number of bed recycles possible 
before heat generation becomes a serious problem. 
For each of three cases, he calculated the heatup 
expected for two burnups: first, 130-day cooled, 
10,000 Mwd/metric ton uranium fuel and, second, 30- 
day cooled, 25,000 Mwd/metric ton uranium fuel. For 
the first case, high-burnup fuel element in still air, 
the temperature is not expected to reach 500°C and is 
therefore not of serious concern. For either of two 
other cases, (1) a high-burnup fuel element immersed 
in a fluid bed and plugged with caked fluidizing medium 
and (2) ahigh-burnup fuel element covered with medium 
during a process shutdown, it is calculated that the 
maximum allowable temperature of 700°C could be 
exceeded in 40 min to 6 hr, depending on the tem- 
perature at which the decladding operation is per- 
formed and the level of fuel-element irradiation. A 
fourth case, a settled bed previously used to process 
irradiated fuel and generating 100 watts of heat per 
liter, would have to be cooled before being allowed to 
settle if it should be desired to reuse the bed for 
several fuel charges. 


Safety Aspects and Environmental 
Contamination Control 


Seven types of metallic micropore filters were 
evaluated at Mol!! for removal of <15-y U,O, and 3- 
to 6-y UF,—UO,F, dust from a gas stream. Such dusts 
will be elutriated in the off-gas from the Mol hot 
semiworks fluoride-volatility plant. On the basis of 
test results and economic factors, a 0.9-mm-thick 
Monel filter fabricated of rigid mesh with 2-u pores 
was chosen for this application. 

Also at Mol,'! caustic liquid scrubbers were eval- 
uated for washing the off-gases and ventilation air 
from the hot fluoride-volatility semiworks. Excellent 
removal of HF and fluorine from off-gases was 
achieved, and these gases and C1F3; were removed 
efficiently from ventilation air. 

Criticality calculations*! were performed at Mol for 
volatility processing of irradiated UO,—PuO, fuel 
from a fast 1000-Mw(e) reactor. The process is one 
of decladding with HF-—O,, fluorination of uranium 
with C1F,, and fluorination of plutonium with fluorine, 
all in fluidized beds. The criticality problems are 
most severe at the finish of the decladding step and 
during fluorination of plutonium. Six situations that 
could occur during processing were assumed, and the 
following factors were varied: the presence of a 


water jacket, weight of alumina bed, number of ir- 
radiated fuel assemblies, thickness of alumina in- 
sulation, and density of bed mixture. The applicable 
multiplication factor, K.,, was calculated for each 
situation. These data will be useful in determining 
the maximum number of irradiated fuel assemblies 
that can be processed simultaneously by fluorination. 

Servais®” has calculated doses that could be de- 
livered from process vessels containing an amount 
of plutonium tetrafluoride (23.24 kg) or plutonium 
oxide (20.00 kg) equal to that contained on one fast 
reactor fuel assembly. From these data he calcu- 
lated the biological shielding required to protect 
personnel working near such vessels. 


Evaluation of Fluoride-Volatility Process 


Allied Chemical plans to construct a fuel-reprocess- 
ing plant on the Savannah River near Aiken, S. C. It 
will have a capacity of 5 metric tons per day and will 
employ a modified aqueous process. *® 

A recent paper*4 by Allied Chemical mentions that a 
technical feasibility study of a conceptual fluoride- 
volatility process (employing aqueous processing for 
plutonium purification) for a Zircaloy-clad low- 
enrichment oxide fuel was made in 1967. The steps of 
the conceptual process were: (1) reaction of the 
Zircaloy cladding with HCl to form volatile ZrCl,, 
(2) oxidation of the fuel with air to pulverize it, (3) se- 
lective fluorination of U,0, with BrF; to form UFg, 
(4) fluorination of PuF, with fluorine to form PuFg, 
and (5) purification of plutonium by conversion to 
PuF, or Pu(SO,),, dissolution in nitric acid, and 
treatment by ion-exchange methods. Extensive use of 
fluid-bed reactors was proposed, and the conceptual 
plant design provided for remote operation and main- 
tenance and extensive shielding. 

The advantages and disadvantages of the conceptual 
design and the relation between plutonium content of 
fuel and plant capacity are discussed.*4 The limiting 
conditions for the conceptual process are a 30-kg 
limit on the mass of plutonium in fluidization reac- 
tor due to criticality considerations, the low rate of 
fluorination of PuF,, and the absence of simple, re- 
liable chemical decladding operations for certain 
fuel-element materials. The advantages of a volatility 
process, in comparison to an aqueous process, were 
stated to be that it can handle fuels after shorter 
cooling times, that krypton and tritium containment 
are easier, that radiation-stable process reagents are 
used, and that conversion of fission products directly 
to stable solid compounds can ease waste disposal. 


Fluorination of Uranium Compounds 
REACTION OF URANIUM COMPOUNDS WITH BrF; 


The kinetics of the gas—solid reactions between 
BrF; and U;0;, UO,, and UO; have been studied in the 
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temperature range 225 to 300°C and in the BrF; 
partial-pressure range 75 to 370 torrs.® All of the 
reactions probably proceed through the formation 
of UO,F, as an intermediate compound. The final 
products of all three reactions are thought to be UFg, 
O,, and Bry». In the reaction of U3;0, with BrF;, UF, is 
formed, probably by the interaction of product UF¢ 
and U;03. The following equations were derived to 
express the dependence of the rate constants on tem- 
perature and BrF; partial pressure for these reac- 
tions: 


U;,0,—BrF;: 
log k’ = 0.90 log P—2000/7 —0.220 
UO,-BrF;: 


log k’ = 0.84 log P —1630/T —0.270 
UO;—BrF;: 
log k’ = 1.05 log P —1680/T —0.767 


Values of the activation energy of 9.2, 7.5, and 7.7 
kcal/mole were calculated for the reactions with 
U;03, UO2, and UOs, respectively. 


REACTION OF UO, WITH FLUORINE 


Iwasaki*® has studied the kinetics and mechanisms 
of the fluorination of UO, pellets to UF, with fluorine. 
Earlier*’ he observed that the reaction proceeds in 
two steps, the conversion of the oxide to UO,F, and 
the conversion of UO,F, to UF.,, with the first step 
the more rapid of the two. In new work he postu- 
lates a three-step reaction: (1) the transformation 
of UO, pellets to partly fluorinated UO, powder 
(fluorine attacks the active sites at the grain bound- 
aries), (2) the formation of an intermediate (identi- 
fied as UO,F, in two cases), and (3) the formation of 
UF,. Below 400°C and with fluorine concentration up 
tc 40%, the highest studied, a thick layer of inter- 
mediate, yellow at the surface and green at the UO,— 
intermediate surface, is formed. Fluorination of 
this intermediate is the rate-limiting step of the 
three. Above 430°C (400 to 430°C is considered a 
boundary region) little intermediate is formed, and 
the conversion of UO, pellets to UO, powder is the 
rate-limiting step. With less than about 15% fluorine, 
the intermediate is black and has not been identified; 
with greater than 15% fluorine, the intermediate is 
white UO,F». 


REACTIONS OF ClF WITH URANIUM COMPOUNDS 


A series of experiments was performed at Mol® to 
determine the effect of reaction temperature (20 to 
300°C) on the nature of the gaseous products of the 
reactions between C1F and various uranium com- 
pounds. For the reaction between UO, and CIF at 
temperatures up to 100°C, the products are UF, Cly, 
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Cl0,F;, and O,; above 100°C, UF., Cl,, and O, are 
formed. When CIF is reacted with U3;03, the same 
products are formed as for UO», and ClO,F is formed 
up to 300°C. Reaction of C1IF with UO,F, gives the 
same products as with U;03, but ClO.F is formed up 
to 400°C and perhaps at higher temperatures. 


Reaction of Uranium Tetrafivoride 
with Water Vapor and Oxygen 


At Mol!! the reaction at 400 to 690°C of powdered 
UF, with water vapor and oxygen was studied by 
thermogravimetry and by chemical and X-ray diffrac- 
tion analyses. The UF, can react with water to pro- 
duce UO, which is oxidized to U;0, or with water plus 
oxygen to produce UO,F, which reacts with water to 
form UO;. A mixture of U;03, and UO3, formed by the 
competitive reactions, is the product. The quantity of 
UO; in the product increases with partial pressure 
of oxygen and decreases with temperature. 


Miscellaneous 


Bibliographies and articles of general interest in- 
clude the following: 

1. A bibliography, compiled by the International 
Atomic Energy Agency, on nonaqueous processing 
of irradiated fuel*® for the period 1955 to June 
1966. Among the subjects covered are the fused-salt 
process, the chemistry of fluorination processes, the 
Nitrofluor process, the chemistry and physical chem- 
istry of volatile fission-product fluorides, the chem- 
istry and physical chemistry of uranium hexafluoride, 
the chemistry and physical chemistry of plutonium 
hexafluoride, decontamination and separation of ura- 
nium and plutonium hexafluorides, and corrosion in 
the fluoride-volatility process. 

2. Cleveland*’ has presented the physical and chem- 
ical properties and method of preparation of many 
plutonium compounds, including the fluorides and 
oxides. 

3. A recently published French book‘? on uranium 
presents the method of preparation, physical proper- 
ties, and chemical properties of uranium compounds, 
including the fluorides and oxides. 

4. Gutmann has compiled a book on halogen chemis- 
try which includes the physical properties of the 
halogens, the inorganic chemistry of the halogens, 
and the physical and chemical properties of the halo- 
gen fluorides. *! Physical, physicochemical, and spec- 
tral studies of fundamental halogen reactions are pre- 
sented together with reviews of halogen-containing 
complexes, the electronic factors governing their 
formation, and the nature and properties of halogen 
bonds. 

5. Thirty-eight papers! on the physical chemis- 
try of protactinium were presented at an interna- 
tional conference at Orsay, July 2—8, 1965. Included 
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were papers concerning the complex halides of 


protactinium(V) 
properties of some alkali 


(D. Brown), the preparation and 


fluoride complexes of 


pentavalent protactinium (L. B. Asprey et al.), and 
the fluoride complexes of tetravalent protactinium 
(L. B. Asprey and R, A Penneman). 
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Compact Pyrochemical Processes 


By W. E. Miller and R. K. Steunenberg 


Pyrochemical processes are a Class of high-tem- 
perature procedures for the recovery and purifica- 
tion of nuclear reactor fuels and for the production of 
materials of interest to the nuclear industry. The 
processes are characterized by the use of molten 
salts and metals as media in which reactions or 
separations are effected. Typical liquid metals em- 
ployed in these processes are Mg, Zn, Cd, Bi, Cu, 
and alloys of these metals. The salts are usually 
mixtures of MgCl, with chlorides of alkali and 
alkaline-earth metals. 

Much of the pyrochemical-processing development 
work is concerned with the processing of fast breeder 
reactor fuel. Pyrochemical processes are particu- 
larly suited for this purpose for several reasons: 
(1) short-cooled fuels can be handled since the high 
radiation doses from these fuels do not adversely 
affect process chemistry; (2) process volumes are 
low, and consequently compact equipment is used; 
(3) solid waste streams of low volume are produced; 
(4) direct decladding of fuel by immersion in molten 
metals offers a promising solution to the decladding 
problem; and (5) the processes are adaptable to 
on-site plants. 


Salt-Transport Process 
for Fast Breeder Reactor Fuels 


Pyrochemical research and development work at 
Argonne National Laboratory is primarily concerned 
with the salt-transport process. Process flow sheets 
for a batch process! and for a process incorporating 
continuous steps” were presented earlier in Power 
Reactor Technology and Reactor Fuel Processing and 
Reactor and Fuel-Processing Technology. 

In the batch process, fuel elements are declad by 
dissolving the cladding in liquid zinc, which is 
subsequently removed from the UO,—PuO, mix- 
ture. The oxide mixture is reduced by magnesium 
in a molten salt—liquid Cu-—Mg system. The alkali, 
alkaline-earth, and rare-earth fission products are 
extracted into the salt phase, which is separated from 


the liquid metal. Uranium is essentially insoluble in 
the liquid metal, whereas plutonium is_ soluble. 
Plutonium is separated from the remaining fission 
products and the uranium by a salt-transport step in 
which the plutonium is removed from the Cu-Mg 
(donor) alloy and deposited in a Zn—Mg (acceptor) 
alloy. The metallic plutonium product is recovered 
by evaporating the Zn—Mg solvent metal. 

It may be advantageous to incorporate continuous 
or semicontinuous steps in the process to provide 
potentially higher capacities and more efficient sepa- 
rations. In the flow sheet that incorporates continuous 
steps, the decladding and reduction steps are essen- 
tially the same as those in the batch process. How- 
ever, the Cu-Mg-— Pu metal phase from the reduction 
step is processed continuously to extract the rare- 
earth fission products into a salt phase; and, in the 
subsequent salt-transport step, plutonium is extracted 
continuously from the Cu—Mg (donor) alloy into a 
salt phase and is simultaneously deposited in a Cd— 
Mg (acceptor) alloy. The plutonium is recovered from 
the Cd—Mg by retorting as in the batch process. The 
objective of the present work at Argonne is to produce 
the technology necessary for the development of a 
plant-scale process. Work is being done on process 
chemistry, process development, corrosion, kinetics, 
and equipment design. 


Zinc—magnesium alloys are used as plutonium 
acceptor alloys in the current pyrochemical flow 
sheet. More data on the solubility of plutonium in 
pure zinc and in Zn—Mg alloys were obtained.* The 
data are plotted in Fig. 1. The line for pure zinc is 
based on earlier experiments,’ but the data points on 
the line are from the recent work. 

As part of a program to investigate mass-transfer 
rates in liquid salt—metal systems, the distribution 
coefficient of cerium between Mg—Cd alloys and a 
55 mole % MgCl,—27 mole % NaCl—18 mole % KCl 
salt mixture was measured.*° The measurements 
were made at 441, 500, 527, 578, and 636°C, with 
magnesium concentrations ranging from 10 to 60 
at.%. The results are plotted in Fig. 2. 
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Fig. 1 Solubility of plutonium in zinc and Zn—Mg alloys. 


An alternative salt—metal system, namely, CaCl,— 
CaF,/Cu—Mg-—Ca alloy, is being considered for the 
reduction step. The solubility of CaO in a CaCl,—20 
mole % CaF, salt mixture was measured’ to deter- 
mine the limits of oxide loading of this salt mixture. 
Two experiments were performed in which CaO in 
excess of the estimated solubility was mixed with 
purified CaCl,—20 mole % CaF, in a baffled tantalum 
crucible at 900 to 950°C. Filtered samples were 
taken on both rising and falling temperature cycles 
after the salt had been equilibrated by stirring at 
each temperature for at least 2 hr. The solubility of 
CaO in the CaCl,-20 mole % CaF, mixture was 
determined by analyses of the samples and is plotted 
as a function of temperature in Fig. 3. A break occurs 
in the curve at about 735°C, and a change in the ratio 
of fluoride to chloride in samples taken above and 
below the break indicates that a solid phase contain- 
ing chloride precipitates below 735°C. 

Another experiment? was performed to study the 
reduction of UO,—PuO, solid solution in the CaCl,— 
20 mole % CaF,/Cu—Mg-—Ca reduction system and to 
determine if plutonium would coprecipitate with ura- 
nium in the reduction. After the reduction the salt 
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and supernatant Cu-—Mg were poured away from the 
uranium precipitate, and the precipitate was dis- 
solved in Zn—11 wt. % Mg alloy. Filtered samples 
were taken of both the Cu-Mg-—Ca alloy during the 
reduction and the Zn—Mg alloy after the dissolution 
of the uranium precipitate. Analyses of the samples 
showed that the plutonium-to-copper ratio in the 
Zn—Mg alloy was the same as that in the Cu-Mg—Ca 
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Fig. 2 Cerium distribution between Mg—Cd alloys and 
55 mole % MgCl,—27 mole % NaCl—18 mole % KCI salt.° 
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Fig. 3 Solubility of CaO in CaCl, —20 mole % CaF, (Ref. 5). 
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reduction alloy. This result indicates that none of the 
plutonium coprecipitated with the uranium, since any 
coprecipitation of plutonium would have resulted in a 
higher plutonium-to-copper ratio in the Zn—Mg alloy. 
The analyses also indicated complete reduction of the 
UO, —PuO,. 

Two decladding experiments were performed to 
study the removal of stainless-steel cladding from 
oxide fuels by dissolution of the cladding in molten 
zinc.® Simulated fuel subassemblies loaded with pel- 
lets and fines of UO, were immersed in molten zinc 
at 800°C. The subassemblies were constructed of 304 
stainless steel and contained 13 tubes, which were °/, 
in. in inside diameter and 12, in. long and had a wall 
thickness of 35 mils. The tubes were welded to three 
round plates 4 in. in diameter by '/, in. thick. The 
subassemblies contained 2.68 kg of stainless steel 
and were loaded with 2.49 kg of UO». 

Each subassembly was dissolved in 84.5 kg of 
molten zine contained in a 9°/,-in.-ID tungsten cru- 
cible. About 6 kg of CaCl,—15 wt. %CaF, salt was 
added to the melt to retard vaporization of the zinc 
and to wet and collect UO, fines, which are reduced 
in the subsequent reduction step. After ~30 min 
immersion at 800°C, the stainless-steel tubing was 
completely dissolved, but the round plates were 
relatively unaffected. The subassembly was im- 
mersed for an additional 7 hr to ensure complete 
dissolution of the stainless steel. Filtered samples 
were taken during the course of the experiments and 
analyzed for iron to follow the progress of the disso- 


lution. The results from one experiment are pre- 
sented in Fig. 4. The results show that a period of 
approximately 4 hr was required to dissolve all of 
the stainless steel in the subassembly. 

After the dissolution was complete, the molten 
zinc—stainless steel solution was transferred from 
the crucible. Subsequent inspection of the crucible 
showed that no undissolved stainless steel remained. 
The UO, pellets were distributed in a fairly uniform 
layer on the vessel’s bottom and were partially 
covered with a zinc heel. 

The zinc remaining in the crucible was removed by 
vacuum distillation. No major operational problems 
were experienced during the decladding and zinc- 
vaporization procedure in either experiment. The 
relatively zinc-free pellets remaining in the crucible 
after the decladding and zinc-removal procedures 
were suitable for the subsequent oxide-reduction 
step. 


Processes for Molten-Salt 
Breeder Reactor Fuel 


A number of recent reports describe the status of 
work directed toward the application of pyrochemical- 
processing techniques to molten-salt breeder reac- 
tor (MSBR) fuel. The molten-salt breeder reactor 
concept, which employs the Th-7y fuel cycle, is 
being developed at Oak Ridge National Laboratory. 
Oak Ridge has built and operated a molten-salt 
reactor experiment (MSRE) fueled with 2357, and 




















2.0 ING 
Oo > 2 
Sis ©) 
1.6 
ae 
= 
5 12 
3 
2 
3 
3 
— Of 
2 
0.4 
0 2 4 6 8 10 12 14 24 26 


Time, hr 


Fig. 4 Progress of decladding experiment. 
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preparations are now being made to operate the 
experimental reactor with °°U fuel. 

Continuous processing concepts for treating molten- 
salt core and blanket fuel were discussed in a pre- 
vious issue of Reactor and Fuel-Processing Tech- 
nology.? In one processing concept, core fuel is 
contacted with fluorine to oxidize uranium from UF, 
to UF., which volatilizes. The barren carrier Salt, 
"LiF -BeF,, is then vacuum distilled, and a highly 
concentrated solution of rare-earth salts in "LiF is 
left as a residue, which is discarded as a fission- 
product waste. The core fuel is reconstituted in a 
single operation in which UF, and H, are introduced 
into the purified salt and the UF, is reduced to UF. 
An alternative to the distillation process uses re- 
ductive extraction of the rare earths from the salt 
after the uranium has been removed from the salt by 
fluorination. One reductant that has been employed is 
lithium dissolved in bismuth. 

Efficient breeding of ayy requires that 233Da (bred 
from **Th) be removed from the neutron flux of the 
reactor to allow it to decay to *°U. The problem of 
separating protactinium from molten fluorides is of 
utmost importance to the thermal breeder concept, 
and the problem is being investigated at Oak Ridge. 
Methods of removing protactinium from blanket salt 
were reviewed previously. °® 

In addition to research and development directly 
related to processing MSBR or MSRE fuels, support- 
ing studies in molten-salt fuel chemistry are being 
carried out at ORNL. Some of these studies are 
concerned with solubility, vapor pressure, viscosity, 
thermodynamic, and electrical-conductivity measure- 
ments in molten-salt systems. 

The distribution of fission products within the fuel- 
containment vessel of a molten-salt reactor is of 
interest to the fuel processor. Not all of the fission 
products remain in the molten salt. Some deposit on 
metal surfaces, some on graphite surfaces, and some 
volatilize into the gas phase within the containment 
vessel. Although the processor is dealing only with 
the salt phase, a knowledge of the fission-product 
distribution helps to determine which fission-product 
elements will be present in the salt. This informa- 
tion is being accumulated from two sources. Opera- 
tion of the MSRE has yielded some information, and 
operation and postirradiation examination of in-pile 
thermal-convection loops have yielded additional data. 
The results of these studies are discussed in detail 
in two recent Oak Ridge reports.” !’ Table 1 gives 
the distribution of some of the fission products in the 
MSRE.* The noble metals, molybdenum, tellurium, 
niobium, and ruthenium, did not behave as expected. 
Appreciable fractions of these fission products de- 
posited on graphite and metal surfaces, and some 
entered the gas phase; only very small amounts 
circulated with the molten-salt stream. 

The relative ease with which an element can be re- 
ductively extracted from 2LiF~- BeF, into bismuth 
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Table 1 APPROXIMATE FISSION-PRODUCT DISTRIBUTION® 
IN THE MSRE AFTER 32,000 Mw-hr 








Inventory 
in MSRE Percent 
(dis/min)* Percent on Percent on Cover gas,f 
Isotope x 10!7 infuel graphite Hastelloy N %/day 
99Mo 7.91 0.94 10.9 40.5 17 
12Te 5.86 0.83 10.0 70.0 66 
108Ru 3.36 0.13 6.6 14.9 40 
95Nb 4.40 0.044 36.4 34.1 5.7 
%7r 6.00 96.1 0.03 0.06 0.14 
88Sr 5.02 77.0 0.26 33 
1317 4.00 64.0 1.0 16 





*This total inventory was calculated from the power history of 
the MSRE, assuming no losses to the cover-gas stream. 

*These values represent the percentages of daily production 
rate lost to the cover gas per day. 


metal has been expressed in terms of extraction 
potentials. **° The potentials that have been deter- 
mined for certain fission products and other fuel 
constituents indicate the following ease of extraction:? 
Li< Be< Ba< Eu < Th< Nd< Ce<Sm< Zr< La< U. 

In a separate experiment it was shown that ura- 
nium could be quantitatively separated from cerium 
by the reductive extraction technique (see Fig. 5). 
It is possible that this procedure could be used as an 
alternative to fluorination for the removal of uranium 
from the salt prior to other fuel-reprocessing steps. 
A similar experiment with uranium and zirconium 
showed a somewhat less effective separation, but one 
that would still be adequate if two or more stages 
were used, 
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Fig. 5 Separation of uranium from cerium by reductive 
extraction with lithium from LiF-—BeF, (67-33 mole %) 
into bismuth at 600°C (Ref. 9). 


The status of the work on the extraction of prot- 
actinium from molten fluorides into molten metals 
is discussed in two recent Oak Ridge National Labo- 
ratory reports.*!° Studies of the removal of prot- 
actinium from solution in LiF—BeF,—ThFy, (73-—2— 
25 mole %) salt by addition of thorium metal were 
reviewed earlier in Reactor and Fuel-Processing 
Technology. ® Earlier experiments were performed 
with static systems in which simulated blanket salt 
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was in contact with bismuth and was contained in 
mild steel. These experiments substantiated the 
reduction of protactinium from the salt mixture but 
failed to demonstrate its quantitative dissolution in 
the metal phase. Since the extraction method utilizes 
the molten metal only as a transport medium, the 
apparently very low solubility of 33a or its carrier 
in the molten metal is not necessarily detrimental to 
the process. A dynamic system was constructed and 
was operated by recirculating bismuth containing 
dissolved thorium through the simulated blanket salt 
mixture. Recovery of protactinium on a steel-wool 
column in the bismuth circuit accounted for about 
43% of the "pa initially in the salt phase. However, 
the recovery of 7°°Pa was effected mainly by filtra- 
tion rather than by adsorption. The filtered solids 
showed that the *°Pa was associated with both iron 
and thorium. 


More recent experiments’? have been designed to 
demonstrate the reductive extraction of *°Pa from 
the blanket salt mixture into a molten metal and its 
back extraction into a second salt mixture by hydro- 
fluorination. The apparatus (shown schematically in 
Fig. 6) contained all of the liquid phases in graphite 
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Fig. 6 Extraction vessel for *%Pa removal from molten 
fluorides.!® (All material in contact with salt and metal 
phases is graphite.) 


to preclude the possible adsorption of *°*Pa on solid- 
metal surfaces. The initial experiment in this series 
demonstrated essentially complete removal of 233. pg 
from the blanket salt, but only 65% of the *°Pa was 
found in the recovery salt. Losses of 7*Pa that 
occurred in this experiment and in a later experi- 
ment were attributed to adsorption of the 33a on 
metal particles formed by reduction of iron impuri- 
ties in the blanket mixture. 

Another experiment was conducted!® in which tin 
was added to the metal phase to increase the iron 
solubility in the bismuth. In this experiment approxi- 
mately 1.78 kg of the simulated LiF—BeF,—ThF, 
(73-2-25 mole %) blanket salt and 1 me of **Pa 
were in contact with 2.98 kg of bismuth that contained 
0.41 kg of tin. Approximately 0.533 kg of LiF—BeF, 
(60—40 mole %) was also placed in contact with the 
molten-metal mixture to serve as a recovery salt. 
Figure 7 shows the transfer of 233 during reduction 
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Fig. 7 Extraction of *%Pa from LiF—BeF,—ThF, (73-2—- 
25 mole %) by reduction with thorium metal into LiF -BeF, 
(60-40 mole%) by hydrofluorination via molten Bi-—Sn 
mixture at 650°C (Ref. 10). 


by thorium-metal additions to the blanket salt and 
hydrofluorination of the recovery salt. Approximately 
90% of the *%Pa activity was found in the recovery- 
salt mixture, 2% was present in the molten-metal 
phase, and 8% remained in the blanket salt. 

One of the proposed methods for processing the 
fuel stream of the MSBR includes a distillation step. 
In this step the LiF—BeF, carrier salt is distilled 
away from the less-volatile fission-product fluorides 
that must be removed from the reactor continuously. 
Current work is producing the information required 
for the design of distillation systems for MSBR 
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processing. Studies of the effect of a concentration 
gradient in a semicontinuous still, of factors limiting 
the vaporization rate, and of errors in measured 
relative volatilities are discussed in a recent ORNL 
report. !! 


Miscellaneous 


Mullins and coworkers’® have reported the dis- 
tribution of “Am between liquid plutonium and a 
fused salt (NaCl—KCl) containing PuCl,. Equilibrium 
measurements were made at 698, 730, and 775°C 
over a PuCl, mole-fraction range of 6 x 10~° to 
2x 1072. 

A recent report discusses the removal of ameri- 
cium from plutonium by electrorefining. ' In 31 
production runs, the plutonium metal containing an 
average of 267 ppm *4!Am was produced from a feed 
metal containing 960 ppm 241Am, An average of 1% of 
the feed plutonium was retained in the salt as the 
chloride and an additional 11.2% was retained in the 
salt as finely divided plutonium metal. 

Applications of molten-salt processing to plutonium 
recovery at Dow Chemical Company, Rocky Flats 
Division, were reported recently by Long.’ Three 
areas of molten-salt processing were investigated, 
namely, electrorefining of plutonium, recovery of 
plutonium from electrorefining salts, and extraction 
of americium from plutonium. 

Neptunium production at Rocky Flats was de- 
scribed in another Dow report.!® Conversion of NpO, 
to NpF, with HF was followed by reduction of the 
NpF, to metal using calcium metal as the reductant 
in a bomb-type reduction. Iodine was added to the 
reduction charge to yield a more fluid slag. The 
average yield for eight reductions was 97.3%. The 
size of the NpO, batches charged was between 100 
and 400 g of neptunium, and the impurities in the 
metal produced ranged from 179 to 1000 ppm. 

A Japanese report discusses the reduction of ThO, 
by calcium.'® The work was performed to determine 
the most favorable conditions for commercial prepa- 
ration of the metal. Unlined stainless steel was the 
most suitable material for a reactor vessel; the 
yield deteriorated when a graphite reactor vessel 
was used. The optimum particle size of ThO, was 
about 150 yp. A high yield of metal could be obtained 
by using an amount of CaCl, flux equivalent to 20 to 
40% of the weight of ThO, and about 25% more cal- 
cium metal than the amount theoretically required 
for the reduction. The reduction was most efficient 
when the system was heated for 4 hr at 950°C or for 
3 hr at 1050°C. Under the above-mentioned condi- 
tions, the yield of thorium metal was about 90%, and 
the amount of residue that was insoluble in HCl was 
about 3 to 7%. 

Another Japanese report!” describes experiments 
in which UO; was reduced by calcium vapor. The 
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resulting uranium powder was melted to determine 
the optimum conditions for preparing a uranium- 
metal powder that had good fusibility and high purity. 
Coarse uranium powder (gray-white) had better 
fusibility than fine uranium powder (black or brown- 
black). The optimum conditions for preparing coarse 
uranium powder were a reaction temperature of 
1100°C, a Ca/UO,; mole ratio of 9, and a reaction 
time of 2 to 3 hr (for this experimental scale). The 
yield percentage of coarse uranium powder was more 
than 90% under these conditions. Further experi- 
ments were conducted on the preparation of a small 
uranium-metal ingot from the uranium-powder prod- 
uct. The best procedure for making small uranium 
ingots was to place the metallic calcium outside of 
the crucible, but within the reaction vessel, and to 
heat the powder in 1 atm argon at a temperature 
greater than 1200°C. The purity of uranium ingots 
obtained under these conditions was greater than 
99.9%. 

A French report!® describes a study of the elec- 
trorefining of beryllium in a molten KCl—LiCl-—BeCl, 
bath at an average temperature of 475°C. The cur- 
rent production rate is 100 g of beryllium per 24 hr 
for a charge of 25 kg of electrolyte. The feasibility 
of obtaining a high-purity metal was demonstrated. 
The sum of the principal metallic impurities of 
commercial beryllium was lowered from a maximum 
of 4000 ppm to 35 ppm after electrorefining. How- 
ever, there is a possibility of contaminating the fused 
salts with the materials constituting the cell, e.g., 
pyrex of silica. Prolonged contact of the baths with 
the cell materials caused contamination, principally 
by oxygen and sodium. Ion exchange between the 
lithium in the bath and the sodium in the pyrex oc- 
curred in the fused eutectic KC1—LiCl at 500°C. It 
was confirmed that the impurities concentrated in 
the immediate vicinity of the cathode. This phe- 
nomenon appeared to occur only with impurities that 
were more reducible than beryllium. 

Corrosion tests from the Molten Salt Reactor 
Experiment are discussed in an ORNL report. The 
first group of Hastelloy N surveillance specimens 
from the reactor core was removed after 7832 Mw-hr 
of reactor operation. During this period the specimens 
were held at 645+ 10°C for 4800 hr, and they accu- 
mulated a thermal dose of 1.3 x 107° neutrons/cm?, 
The high-temperature ductility of these specimens, 
which were irradiated in a salt environment, was 
reduced, but the decrease was similar to that ob- 
served for Hastelloy N specimens irradiated in the 
Oak Ridge Research Reactor in a helium environ- 
ment. The low-temperature ductility was also re- 
duced, and this effect was attributed to the growth of 
intergranular precipitates. The specimens showed no 
evidence of corrosion. However, a carbon-rich layer, 
1 to 2 mils in depth, was noted in areas where 
Hastelloy N and graphite were in contact. The me- 
chanical properties of the Hastelloy N appear to be 
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adequate for the continued satisfactory operation of 
the MSRE. Test results are presented for the effects 
of several variables on the tensile ductility of irra- 
diated and unirradiated Hastelloy N. These variables 
included test temperature, strain rate, and pre- 
straining. 

Pyrochemical preparation of U—3 wt.% Mo alloy 
powders was reported at a recent American Chemi- 
cal Society meeting.”’ The alloy powder, which had a 
particle size of —200 mesh, was prepared for powder- 
metallurgy applications by bomb reduction of mixed 
uranium and molybdenum oxides with ground cal- 
cium. The oxygen content of the mixed oxides ranged 
from 12.1 to 15.3 wt.%. In some reductions, 0.1 to 
1.0 mole of CaCl, per mole of uranium was used to 
reduce the average particle size of the alloy powder 
and to aid in dissolution of the slag. The particle- 
size distribution, molybdenum content, and purity of 
the alloy powders were compared with the oxygen 
content and purity of the oxide feeds. A higher oxygen 
content in the mixed oxides appeared to lead toa 
coarser metal alloy powder, other things being equal. 
In some cases the molybdenum content was higher in 
the larger-sized fractions but was lower when the 
larger particles consisted of filigree agglomerates 
of small particles. The metal alloy powders con- 
tained significantly less Na, Mg, Si, Mn, Cu, and Pb 
than the mixed-oxide feeds. The best yields of —200- 
mesh metal alloy powder, and the best-controlled 
molybdenum content, were obtained using a mixed- 
oxide feed with an oxygen content of 12 wt.%. 


A French report?! discusses an electrolytic method 
for the separation of uranium and plutonium fluorides 
in LiCl—KCl eutectic. The purpose of the work was 
to perfect the separation of the uranium and plutonium 
products of the fluoride-volatility reprocessing cycle 
and to obtain UO, of high purity which could be used 
directly in the refabrication of nuclear fuels. The 
first step of the procedure consists in dissolving a 
mixture of the two fluorides in the LiCl—KC1 eutectic 
at 450°C. The solution is then oxidized to obtain the 
soluble uranium oxyhalogenide and to precipitate 
plutonium in the form of insoluble PuO,. Uranium 
dioxide is deposited at the cathode by electrolysis. 
The residual PuO, is treated with hydrogen chloride 
(mixed with hydrogen if necessary), which renders it 
soluble in the form of PuCl;. The plutonium metal can 
be recovered at another cathode by electrolysis at a 
higher potential. This new method has the advantage 
of using known experimental techniques. The separa- 
tion tests carried out on uranium and cerium (a 
stand-in for plutonium) gave separation factors of 
100. 


Another French publication” reports the prepara- 
tion of UO, crystals by electrolysis of a bath con- 
sisting of 20 wt.% UO,F, dissolved in LiCl—KCl 
eutectic at 450°C. Of the UO, powder obtained, 80 to 
90% was in good crystalline form. 
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WASTE DISPOSAL 


Progress in Research and Development 
on Waste Disposal 


By Phillip Fineman 


Conversion of High-Level-Activity 
Wastes to Solids 


FLUID-BED CALCINATION 


In March 1968 the Idaho Chemical Processing Plant 
(ICPP) terminated on schedule its second waste- 
processing campaign! in which the aqueous fuel- 
reprocessing wastes that are produced there were 
treated by fluid-bed calcination in the Waste Cal- 
cining Facility (WCF) at the ICPP.’ 

During this campaign, which had been in contin- 
uous operation for nearly 2 years, 880,000 gal of 
aluminum nitrate-bearing wastes and 103,000 gal of 
fluoride-bearing zirconium wastes were processed 
through the WCF and converted to 12,200 cu ft of 
granular solids for an average overall volume reduc- 
tion of 11. The calcined product is stored in under- 
ground bins. The process problems encountered in 
the campaign were minimal, and these were corrected 
by adjusting operating parameters. °® 

Since the completion of the campaign, the WCF has 
been placed on standby status, and preventative 
maintenance and minor repairs of WCF equipment 
are under way.” 


POT CALCINATION PROCESS 


A recent Pacific Northwest Laboratory (PNL) 
progress report* describes continuing developmental 
work on the pot calcination process® for the solidifi- 
cation of high-level-activity wastes. This process, 
which was originally developed at Oak Ridge National 
Laboratory (ORNL), is one of three solidification 
processes that PNL is planning to demonstrate on 
an engineering scale with plant-level-activity wastes, 
Treatment of wastes by the ORNL process involves 
carrying out the calcination and solidification steps 
in the same vessel or pot. The pot also serves as 
the permanent storage container for the solidified 
waste product. 


Two additional demonstration runs‘ (PC-4 and 
PC-5) were completed using simulated PW-1 wastes* 
prepared from Purex-process plant wastes. Arrange- 
ment of the equipment (designated as Mode C form 
of operation) for the two runs is shown in Fig. 1. In 
this arrangement the calciner pot is fed directly from 
the waste evaporator, and all of the pot condensate 
is continuously recycled directly to the evaporator. 
Thus a constituent that is volatilized from the pot 
will continue to increase in concentration in the feed 
stream until a steady-state concentration is attained. 
Waste solution is continuously supplied to the evap- 
orator to replenish that fed to the calciner pot. In 
run PC-4, sodium hydroxide solution (4M) was fed 
to the top plate of the acid fractionator to determine 
if scrubbing of the fractionator off-gas would remove 
additional ruthenium. 

In the two runs the recycling of the pot condensate 
resulted in decreased volatilization of ruthenium. 
Less than 5% and 3%, respectively, of the total ru- 
thenium in the feed to the pot was volatilized in runs 
PC-4 and PC-5. In prior runs conducted with PW-2 
waste,t 30% of the ruthenium volatilized. Although 
the reason for the low ruthenium volatility with PW-1 
waste is not known, it is thought that the high iron 
content (~0.8M) of PW-1 wastes may be a key factor. 
The use of sodium hydroxide solution in run PC-4 
to scrub the fractionator off-gas resulted in a ru- 
thenium decontamination factor (across the fraction- 


*PW-1 waste is a first-cycle Purex-process waste. It 
is high in nitric acid (~5M), free of sulfate ion, and con- 
tains added iron (~0.8M) from the dissolvable container 
used to transport chopped fuel. Because the plant wastes 
available for preparation of the PW-1 waste feeds for 
these two runs contained excessive amounts of sodium, 
sulfate was added to prevent the formation of volatile 
sodium and cesium oxides. 

*PW-2 waste is a combined first- and second-cycle 
Purex-process waste. It is high in nitric acid (~5M) and 
contains iron (~0.4M) and sulfate (~0.8M) ions from the 
ferrous sulfamate reductant. 
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ator) that was 30 times greater than that for run 
PC-5, in which no sodium hydroxide solution was 
used. 

Although operating difficulties (foaming and plug- 
ging) were encountered in the two runs, the runs were 
successfully completed. Early in run PC-4, foaming 
took place in the calciner pot and caused the expul- 
sion of some of the feed solution to the safety-relief 
pot. To prevent further expulsion of feed, the solution 
level in the pot was lowered to allow more freeboard 
for foaming; however, no further foaming was ob- 
served during the remainder of the run. In run PC-5, 
sufficient freeboard for foaming was maintained. 


Run PC-5. In this run an 8-in.-OD pot, also fab- 
ricated of 304L stainless steel, was used for the 
calcination of 520 liters of feed containing 1,200,000 
curies of mixed-fission-product activity (including 
18,000 curies of ruthenium). A calcine volume of 55 
liters was produced for a volume reduction of about 
9. This waste would be the equivalent of processing 
about 1.3 metric tons of 20,000-Mwd power-reactor 
fuel that had been allowed to cool for 1.5 years. 

Internal heat generated (total of 4300 watts or 78 
watts/liter) by the radioactivity in the pot was about 
85% of the maximum thermal amount permitted for 
a pot containing solidified waste; at the maximum 
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Although foaming was evident, no feed was expelled. 
Plugging of the sparger and the specific-gravity 
dip tube in the waste evaporator occurred in both 
runs; the plugging condition was not corrected during 
the course of the runs. Additional details for each of 
the runs are presented below. 


Run PC-4. A total of 1080 liters of feed contain- 
ing 650,000 curies of mixed-fission-product activity 
(of which 47,000 curies was ruthenium and 30,000 
curies was strontium) was processed in a 12-in.-OD 
pot fabricated of 304L stainless steel. The volume of 
calcine produced was 105 liters, which corresponds 
to a volume reduction (feed to calcine) of about 10. 
This waste would be the equivalent of processing 
about 3 metric tons of 20,000-Mwd power-reactor 
fuel that had been allowed to cool for 4.7 years. 

The self-generating heat density of the calcine 
product was 25 watts/liter. At equilibrium (ambient 
air temperature of 38°C), the maximum centerline 
temperature of the pot was 477°C, and the tempera- 
ture drop from the center line of the pot to its wall 
was 312°C, 
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Equipment arrangement (Mode C operation) for PNL pot solidification for runs PC-4 and 5 (Ref. 4). 


value, upon which the design criteria for pot storage 
are based, the internal pot temperature of 900°C or 
the pot-wall temperature of 425°C would not be ex- 
ceeded. At equilibrium the pot centerline temperature 
(ambient air at 38°C) was 725°C, and the temperature 
difference between pot center line and its wall was 
495°C. 


CONVERSION TO GLASS 


The results of a design and cost study of a glass- 
fixation process developed at Harwell were presented 
at the 1967 American Institute of Chemical Engineers 
Symposium on Recent Advances in Reprocessing of 
Irradiated Fuel.’ The process known as FINGAL 
(Fixation in Glass of Active Liquor)® incorporates 
highly radioactive nitric acid wastes for permanent 
storage in a chemically inert and highly insoluble 
glass. These wastes arise from the reprocessing of 
fuel from the U. K. nuclear power program. The 
FINGAL process, which has been under development 
at Harwell for some 10 years, is being considered as 
an alternative to the presently used system of liquid 
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storage.’ A simplified flow sheet of the process is 
shown in Fig. 2. Processing of the nitric acid wastes, 
which have been previously concentrated by evapora- 
tion, is carried out on a batch-production basis and 
involves using three cylindrical vessels* of type 
25/20 stainless steel connected in series. The con- 
centrated wastes, to which are added silica and 
boraxf (the glass formers), are fed to the first 
vessel maintained at 1050°C. Here the slurry mix- 
ture is dried and the solid residue fused into a glass, 
Process off-gas from the glassmaking vessel is 





























stream.{ For successive new process cycles, the 
other two vessels are progressively moved into the 
glassmaking position; the filter in this vessel is 
then dropped and incorporated in the glass melt. A 
new vessel with a new filter is installed in the pro- 
cessing train. (Design drawings showing details ofthe 
proposed glass vitrification plant and the storage 
vault are given in the original paper,’ to which the 
reader is referred.) 

The authors compared the costs of building the 
FINGAL process plant and operating it for a 20-year 
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passed successively through the other two vessels, 
each of which is equipped with a filter to remove 
volatilized ruthenium and other fission products, and 
then the off-gas is scrubbed with nitric acid and 
sodium hydroxide before discharge tothe atmosphere. 
When the process cycle is complete, the vessel con- 
taining the solidified glass product is removed from 
the furnace and transferred to dry storage in a vault 
where the glass-filled vessel is cooled by an air 





*The vessels are 1 ft in inside diameter by 8 ft high and 
will be filled to a glass depth of 6 ft. Each vessel will hold 
the waste arising from the processing of approximately 
17 metric tons of Magnox fuel. A vessel of this size rep- 
resents a compromise between obtaining the maximum 
amount of surface area available for heat transfer from 
the glass and reducing the number of cylinders which have 
to be produced from a large reprocessing program. 

*Tests at Harwell have shown that up to 40 wt.% of waste 
oxide can be incorporated in a B-—Si glass without ad- 
versely affecting its properties as a containment medium.’ 


2 Simplified flow sheet of Harwell FINGAL process. 


7 


period with that of the existing method of storing 
concentrated acid wastes in stainless-steel tanks 
equipped with cooling coils. The capital and operating 
costs for the two schemes are based on the following 
assumptions: 


1. The reprocessing of 1500 metric tons of Magnox 
fuel per year. 


tThe following thermal values were established for the 
storage of glass-filled vessels in dry vaults: (1) maximum 
centerline temperature of the filled vessel, 500°C; (2) max- 
imum temperature at the surface of the vessel, 200°C; 
and (3) maximum cooling-air outlet temperature, 200°C. 
Sufficient cooling air (natural draft supplemented with fans) 
must be passed across the stored vessels of glass for a 
long period to maintain the temperatures below these 
values. The natural draft is furnished by a chimney 200 ft 
high by 10 ft in inside diameter, which can serve five 
storage vaults. Any one storage vault is capable of con- 
taining a total of 488 vessels stored in rows of 8 abreast 
with a maximum of 61 rows. 
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2. The tank storage of concentrated nitric acid 
waste (concentrate volume of 18 liters per metric 
ton of Magnox fuel). 

3. The availability of one standby storage tank. 

4, Present-day costs. 

5. An operating period of 20 years. 





Cost, millions of pounds (British) 


FINGAL process 
and glass storage 





Concentrated 
liquid storage 








Capital value of facilities 3.78 2.41 
Operating costs 1.74 2.62 
5.52* 6.03f 





*Equivalent to $12.5 million (U. S.). 
+Equivalent to $16.9 million (U. S.). 


The authors point out that, of the two schemes, 
liquid storage has a cost advantage, an advantage that 
would be very much greater if the liquid wastes were 
evaporated to higher concentration factors. However, 
some of the cost advantage would eventually be lost, 
once it became necessary to start replacing liquid- 
storage tanks. About 20% of the cost of operating the 
FINGAL process is attributed to the cost of the 
stainless-steel vessels; and, although no method of 
reducing this cost is apparent at present, this may be 
possible in the future, They also point out that this 
process has now been sufficiently developed in the 
United Kingdom so that there is enough information 
to utilize it directly on a substantial scale, and there- 
fore no further extensive development is planned, 
However, the cost of introducing this process is not 
considered justifiable on safety grounds alone, and it 
is believed necessary to demonstrate that there are 
economic advantages to the system before it should 
be chosen to replace the already established liquid- 
storage system. 


Selected Current Literature 


Studies related to the storage of high-level-activity 
wastes in underground storage tanks are described 
in two reports, one issued by PNL’ and the other by 
SRL (Savannah River Laboratory).'’ In the PNL 
report,’ Willingham describes the use of a computer 
program developed for calculating steady-state tem- 
peratures in heat-generating and non-heat-generating 
systems. This program offers a reasonably accurate 
method of estimating temperatures expected in and 
around underground storage tanks containing radio- 
active wastes. These temperatures are important to 
know so that safety limits set by past experience are 
not exceeded. Programs with constant or tempera- 
ture-dependent thermal conductivities, internal heat 
generation, and constant film or contact coefficients 
can be accommodated. The results are calculated by 
an iterative process, which the author describes in 
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detail. In the SRL report,'? Fenimore and Horton 
have examined the possible consequences of a release 
of waste due to tank leakage into the surrounding soil 
at SRL. (Tank design and operating controls minimize 
the possibility of this occurring at SRL.) In their lab- 
oratory experiments, they passed simulated wastes* 
containing radioactive cesium and strontium through 
small columns (2 in. in diameter by 6 in. long) of 
SRL soil, a sandy clay soil that surrounds the storage 
tanks. The rate of movement of cesium and strontium 
through the soil and the feasibility of recovering these 
radionuclides from the soil by leaching were studied. 
Their results indicated that, if radioactive waste 
leaked from a storage tank, the movement of stron- 
tium and cesium through the surrounding soil would 
be retarded because of reduced permeability? of the 
soil. The movement of strontium would accelerate as 
the waste became diluted with groundwater, but 
cesium would continue to move at a low rate. Stron- 
tium was readily leached from the waste-soil mixture 
by dilute solutions{ of acid or salts, but cesium was 
not readily leached by any solution tested. 

A series of articles devoted to the solvent extrac- 
tion of useful fission products generated in the re- 
processing of spent nuclear fuel is being presented 
in Isotopes and Radiation Technology, a Technical 
Progress Review published by the AEC.'! The sec- 
ond? of these articles reviews the liquid — liquid 
extraction of cesium. Details are presented of spe- 
cific extraction schemes as well as safety considera- 
tion of the organic solvents used. 

Blomeke and Harrington,'® ORNL, review the waste- 
management experience of several nuclear power- 
plant prototype stations (Dresden 1, Big Rock Point, 
Humboldt Bay, Elk River, Indian Point 1, and Yankee). 
The report discusses the sources and characteristics 
of wastes of these power stations and their waste- 
management systems. Descriptions of the systems 
and operational and cost experiences are presented. 

Proceedings" have been published of the 1967 
International Atomic Energy Agency symposium’® on 
the disposal of radioactive wastes into the ground. 
English translations'® of eight Russian papers pre- 
sented at the symposium have become available. 

ORNL has developed a computer (FORTRAN II) 
code for handling the vast amount of data obtained 
in its program!" for the fixation of highly radioactive 
wastes into stable, solid media. This code was de- 





*High-level-activity liquid wastes at SRL are essentially 
an alkaline solution of sodium salts, and consist, typically, 
of the following constituents: NaNO, (3.5M), NayCO; (0.3M), 
Na,SO, (0.2M), NaOH (0.3M), and NaAlO, (0.6M), The 
wastes also contain such fission products as: Sr, ®zr— 
SNb, Ru, Ru, 13’Cs, 41Ce, and Ce. 

‘Soil permeability was reduced due to dispersal of clay 
aggregates by sodium ions from the waste solution. 

iThe solutions used for the leaching tests were untreated 
SRL groundwater and 0.1N solutions of KNO3, HNO;, and 
Ca(NOs)o. 








WASTE DISPOSAL 231 


signed to provide data storage, retrieval, and inter- 
pretation. Details of the computer code were presented 
in a previously issued ORNL report.'® ORNL has 
recently issued a supplement!® to this report in which 
is described an option for use with the computer 
code. The option was developed to permit arranging 
the data in order and comparing them with data col- 
lected either in the earlier part of the program or 
any data that may subsequently be stored in the pro- 
gram. The option can be used to select and compare 
solid compositions on the basis of their properties, 
the additives used in their formulation, and chemical 
equivalents (molecular weight of an element divided 
by the oxidation state of the element) of various ions 
in the solids. 

The International Atomic Energy Agency has issued 
an annotated bibliography”” on the recovery of fission 
products. The bibliography contains 720 references 
to the literature published during the period 1960 to 
1966. 
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